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2.6    2GCF  ROUTINES  (Continued) 


2.6.75    ROUTINE  U 

2.6.75.1  Purpose 

The  purpose  of  this  routine  is  to  compute  the  difference  between  the  present 
controller  offset  and  the  desired  controller  offset  and  correct  for  this 
difference.     This  correction  is  limited  to  between  +6  seconds  and  -3  seconds 
per  cycle,    so  the  necessary  compensation  may  extend  over  several  cycles. 

2.6.75.2  Relationship 

This  routine  is  entered,    if  required,    from  Routine  T,    the  Interval  Advance 
Routine  which  is  executed  once  per  second  for  each  local  controller  in  the 
system.     Routine  U  is  called  by  Routine  T  at  the  start  of  the  programmed 
A-Phase  Green  provided  that  no  transition  is  requested  or  is  in  progress. 

2.6.75.3  Assumption 


None, 


2.  6.  75.4    Approach 

This  routine  is  entered  at  the  start  of  the  programmed  A- Phase  Green, 
provided  that  a  transition  has  not  been  initiated  or  is  not  in  progress  at 
this  controller.     The  routine  will  perform  the  following  series  of  compu- 
tations. 


1.  The  difference,  EPSOF,   between  the  present  controller 
offset,  TCK1,   and  the  desired  offset  value,   WRKOF,  will 
be  computed. 

2.  EPSOF  will  be  converted  so  that  the  shortest  time  between 
the  two  offset  values  will  be  chosen.     The  converted  cor- 
rection will  then  be  imposed  over  a  series  of  cycles,  with 
the  constraint  that  a  given  cycle's  correction  will  always 
be  greater  than  -3  seconds  and  less  than  +6  seconds  to 
minimize  the  transitional  effect. 

3.  EPSOF  will  be  added  to  the  first  interval  duration. 


2.6.75.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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C      EXIT      ^ 

^    TO  ROUTINE  T   J 

Routine  U  (Concluded) 


2.6.  76    ROUTINE  V 

2.6.76.1  Purpose 

Routine  V  apportions  the  split  of  the  critical  phases  of  an  operational 
critical  intersection  according  to  the  measure  (detected)  vehicular  demand 
on  these  phases. 

2.6.76.2  Re  lationship 

Routine  V  is  entered  once  per  second  for  each  controller  being  processed. 
The  routine  is  part  of  the  Level  2  controller  Do- Loop  Routines.     It  follows 
Routine  T  and  exits  to  Routine  Al.     For  every  controller  satisfying  the 
critical  intersection  criteria,    it  calls  Subroutine  UPKV  which  unpacks  K3PH 
and  LSI  for  that  controller.     It  also  unpacks  the  four  values  of  GAABC, 
GLABC,    GMABC,    SABCT,    and  GUABC  corresponding  to  that  controller. 
Subroutines  GDSUB  and  PSUB2  are  also  called. 

2.6.76.3  Assumption 
None. 

2.6. 76.4  Approach 

Routine  V  is  one  of  the  routines  comprising  the  controller  Do- Loop  (defined 
in  Routine  Al,    Do- Loop  Routine)  and  is  scheduled  by  the  Do- Loop  Routine 
at  the  beginning  of  the  first  programmed  controller  interval.     Routine  Al 
passes  unpacked  bit  arrays  for  32  controllers.     Routine  V  checks  each  of 
these  32  controllers  to  see  if  the  Interval  One  Start  Flag  and  semi- actuated 
controller  switches  are  set  (INT1S  and  SEMSK).     The  routine  will  perform 


the  critical  intersection  phase  split  computations  for  these  controllers  pro- 
vided all  the  following  conditions  are  satisfied. 

•  The  intersection  has  been  designated  by  the  operator 
as  a  critical  intersection  ("CICSF"  FLAG). 

•  The  intersection  has  the  required  detector  complement 
for  CIC  operation  ("NQL"  number  of  critical  links 
array). 

•  The  intersection  may  operate  as  a  CIC  at  the  present 
time  of  day  ("CTODF"  FLAG). 

•  The  detectors  on  all  critical  links  feeding  the  inter- 
section are  functioning  properly. 

•  The  intersection  is  not  starting  either  a  transitional 
cycle  or  the  first  cycle  following  a  transition  ("CICBP" 
FLAG). 

•  The  demand  on  at  least  one  of  the  critical  links  has 
exceeded  a  threshold  constant  stored  in  arrays  GUABC 
and  GLABC  as  described  in  the  following  paragraphs. 

•  A  normal  operator  offset  change  has  not  been  introduced 
at  the  start  of  this  cycle  ("FSTCF"  FLAG). 

If  the  above  conditions  are  satisfied  Routine  V  performs  the  following  com- 
putations : 

•  The  green  demand  time  (GD)  is  computed  for  each  selected 
critical  link  feeding  the  critical  intersection  by  calling  the 
Green  Demand  Subroutine. 


•  If  a  controller  phase  has  two  critical  links  assigned  to  it 
(indicated  by  contents  of  the  "LDSET"  array)  the  green 
demand  time  will  be  set  equal  to  the  larger  of  the  two 
link  demand  times. 

•  The  green  demand  time  on  all  of  the  critical  links  are 
checked  to  determine  if  at  least  one  exceeds  an  upper 
threshold  ("GUABC"  array)  before  initiation  of  CIC 
control  begins.     To  maintain  CIC  Control,   once  ini- 
tiated, at  least  one  of  the  critical  links  must  exceed  a 
lower  threshold  ("GLABC"  array).     These  two  thres- 
holds provide  a  hysteresis  correction  to  prevent  switching 
back  and  forth  between  area  and  critical  intersection 
control. 

•  The  adjustable  intervals  associated  with  each  critical 
phases  will  be  corrected  by  the  amount  of  change  in  that 
phase. 


2.6.76.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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13 


IR9  =  GAABC(I'R4> 
IR12  =  GMABC(IR4) 


IR7  = IR4+1 

IR13  =  GMABC0R7) 


©■ 


PSUB2      SUBROUTINE 
FOR  GA,  GB 


IR9  =  SABCTUR4) 


IR12  =  IR12-IR9 


IR7  =  VRINKD+IR6 


INTDUUR7)  =  IR12 


INDX5  =  VRINT(2)+IR6 


INTD4(INDX5)  =  GB-SUMB 


CNT0R(K12)  =  INIDU0R6) 


Routine  V  (Continued) 


14 


3440 


3560 


IR7  =  IR4+1 


SUMC  = 
SABCTIIR7) 


GA  =  GAMIN 

GB  =  GBMIN 

GC  =  GABC-GAMIN- 

GBMIN 


IR7  =  VRINTC3) 

+IR6 


INTDU(IR7) 
GC-SUMC 


GDI  = 

GD2  = 

=  GDB 
=  GDA 

1 

PSUB2 
SUBROUTINE 
FOR  GC,  GA2 

1 

G12  =  GAABCUR4) 
GAB  =  G12 

i 

GA  =  MINIMUM 

OF 

GA1,  GA2 

1 

IR7  =  IR4+1 

1 

GB  =  GAB-GA 

J 

1  1770 

G1MIN  = 
GMABC  =  IR7 

<2> 

1 

G2MIN  = 
GMABC(IR4> 

1 

PSUB2 
SUBROUTINE 
FOR  GB,  GA1 

J 

GDI  =  GDC 
GD2  = GDA 

i 

G12  = 

GAABC0R7) 

1 

IR7  =  IR7+1 

i 

G1MIN 
SABCKIR7) 

i 

G2MIN  = 
GMABC0R4) 

Routine  V  (Continued) 


15 
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Routine  V  (Concluded) 
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2.6.  77    SUBROUTINE  GDSUB 

2.6.77.1  Purpose 

This  subroutine  computes  green  demand. 

2.6.77.2  Relationship 

This  subroutine  is  called  from  Routine  V  by: 

CALL  GDSUB  (Index,    GA) 

where: 

Index  specifies  the  link 

GA  returns  the  computed  demand  in  seconds. 

This  subroutine  calls  UPHW. 

2.6.77.3  Assumption 
None. 

2.6.77.4  Approach 

Subroutine  GDSUB  computes  vehicle  demands. 

2.6.77.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 


18 


C      ENTER  J 


CALL  UPHW 

(AHCCI,  CNTLN,  AHCC) 


I 


CALL  UPHW 

(VOLI,  CNTLN,  VOL) 


I 


AHCC  = AHCC/128 
IA  =  AHCC*  Cl/32 
IA  =  IA/256 


I 


ICC2  =  VOL*  KCL(CYCLC) 

7200 
IB  =  ICC2*C2/2 
ICROS  =ICC2*  AHCC 


I 


IC  =  ICROS*  C3/16 
IC  =  IC/512 

GA  =  IA  +  IB  +IC 
GA  =  GA*IRR/32 


I 


C      RETURN      J 


Subroutine  GDSUB 
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2.6.78    SUBROUTINE  PSUB2' 

2.6. 78. 1  Purpose 

This  subroutine  manipulates  the  arguments  to  compute  a  split  with  respect 
to  the  ratio  of  two  green  demands  and  checks  if  within  bounds. 

2.6.78.2  Relationship 

This  routine  is  called  from  Routine  V  by: 

CALL  PSUB2  (IR9,  IR10,  IR11,  IR12,  IR13) 

where: 

IR9  =  Green  time  available 

IR10  =  A-phase  green  demand 

IR11  =  B-phase  green  demand 

IR12  =  Minimum  A-phase  green  duration 

IR13  =  Minimum  B-phase  green  duration 

2.6.78.3  Assumption 
None. 

2.6.78.4  Approach 

Subroutine  PSUB2  allocates  green  time  proportional  to  the  green  demands, 
so  that  both  the  A-phase  and  B-phase  durations  are  greater  than  their 
respective  minimum. 

2.6.78.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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IR12  =  IR9  -  IR13 


c 


I 


RETURN 


J> 


C       ENTER   J 


ROUND  IR10 
AND  IR11  DOWN 
TO  A  EVEN 
NUMBER 


I 


IR11  = 
(IR11XIR9) 

I  RIO 


I 


ROUND  IR11 
DOWN  TO  A 
EVEN  NUMBER 


YES 


< 


RETURN 


J 


IRIO  =  IR9  -  IR11 
IR12  =  IR11 


YES 


IR13=  IRIO 

1 

c 

RETURN 

) 

Subroutine  PSUB2 
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2.6.79    SUBROUTINE  UPKV 

2.6.79.1  Purpose 

Subroutine  UPKV  unpacks  bytes  in  a  manner  similar  to  Subroutine  UN  PS. 

2.6.79.2  Re  lationship 

This  subroutine  is  called  from  Routine  V  by: 

CALL  UPKV 
This  subroutine  calls  subroutines  UPBYT  and  UPHW. 

2 .  6 .  79 .  3    Assumption 
None. 

2.6. 79.4    Approach 

Subroutine  UPKV  makes  repeated  calls  to  Subroutine  UPBYT  to  unpack  data 
used  by  Routine  V. 

2.6.  79.5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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& 


c 


ENTER 


I 


") 


CLEAR  LINK 
ARRAY 


I 


CALL  UPBYT 
CK3PHI,  KL2,  K3PH) 


I 


LSI  =  6  *  KL2 
1=  LSI 
LSO=  LSI +  6 


CALL  UPHW 
CCNTLI,  I,  LL) 


YES       /£& 


Subroutine  UPKV 
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CALL  UPBYT 
(LNKPT,  LL+1,  IP) 


Subroutine  UPKV  (Continued) 
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CALL  UPBYT 
(CYCLI,  KL2,  CYCLC) 


I 


CYCLC  =  CYCLC +1 
jzz  4  *|<L2  -3 
K  =  l 


i 


CALL  SUBROUTINE  UPBYT 
TO  TRANSFER  DATA 
FROM-.  TO 

GAABI  (J)  -  GAABC  (K) 

GLABI  (J)-,  GLABC  (K) 

GMABI  (J)  -  GMABC  (K) 

SABCI  (J)  -  SABCT  (K) 

GUABI  (J)  -  GUABC  (K) 

VRINI  (J)  -.  VRINT  (K) 


Subroutine  UPKV  (Concluded) 
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2.6.80    ROUTINE  W 

2.6.80.  1    Purpose 

Routine  W  provides  logic  for  selecting  one  of  the  following  area  control 
techniques  for  each  section: 

•  Time  of  day  selection  of  pattern 

•  Manual  operator  selection  of  pattern 

•  Automatic  selection  of  pattern  by  traffic  responsive 
algorithm 

•  Inhibits  CIC  operation  during  certain  periods  in  the  day 
at  those  intersections  where  prevailing  conditions  require 
such  action. 

The  routine  re-formats  controller  parameters  from  disk  into  arrays  for  use 
in  other  routines.     It  re-formats  the  link  parameters  used  by  Routine  V 
and  it  sets  up  flags  and  arrays  required  by  the  controller  and  data  proc- 
essing routines. 

2.6.80.2    Relationship 

Routine  W  is  entered  once  every  15  minutes  or  whenever  the  operator  calls 
for  a  new  Manual  or  Time-of-Day  Pattern(s).     The  routine  is  scheduled 
by  Routines  A,   H,   or  M  and  entered  from  Routine  A.     It  exits  to  Routine  A. 
It  calls  Subroutines  WSUB  and  TRSP. 
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2.6.80.3  Assumption 

Data  is  stored  on  the  RAD  in  appropriate  formats. 

2. 6.80.4  Approach 

The  primary  change  to  Routine  W  was  to  delete  the  first  generation  traffic 
responsive  code  which  selected  traffic  signal  timing  patterns  through 
matching,    link  by  link,    surveillance  volume  plus  weighted  occupancy  with 
volume  plus  weighted  occupancy  associated  with  each  first  generation, 
traffic  responsive,    and  signal  timing  pattern.     Code  to  schedule  the  second 
generation  traffic  signal  optimization  program  (LSTSQS)  was  inserted  in 
place  of  the  deleted  code. 

The  second  generation  optimization  scheduler  executes  if  either  the  desired 
mode  or  current  mode  is  traffic  responsive.     This  scheduler  merely  reads 
in  the  time-of-day  pattern  from  the  RAD  and  sets  the  Level  5  Start  Flags 
(LEV5SF)  for  those  sections  which  have  been  selected  for  the  traffic 
responsive  mode  of  control.     These  flags  in  turn  are  examined  by  the 
MACRO  loop  executive  when  it  is  executed  at  the  5MINCT  mark.     If  any  are 
set,  the  Level  Five  executive  calls  the  optimization  routines.     The  time-of- 
day  pattern  is  required  to  obtain  the  durations  of  the  fixed  intervals  in  each 
controller  cycle. 

If  the  operator  selects  the  manual  or  time-of-day  mode,    Routine  W  resets 
the  LEV5SF  flags  for  the  appropriate  sections.     This  reset  of  the  flags 
inhibits  the  Level  5  optimization  routine  for  the  appropriate  sections. 

It  should  also  be  stated  that  Routine  W  no  longer   sets  the  Transition  Flags 
(TRANSF)  in  the  traffic  responsive  mode.     This  function  is  performed  when 
necessary  in  the  Priority  Level  5  software. 
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In  addition,  Routine  W  no  longer  waits  to  execute  the  traffic  responsive 
mode  until  the  next  15 -Minute  Mark  (FMMARK).     Instead,  the  LEV5SF 
flags  are  set  immediately  upon  request  of  the  traffic  responsive  mode.    The 
Level  5  optimization  algorithm  is  not  executed  immediately,  however.     In- 
stead it  will  be  executed  at  the  next  5MINCT  mark. 

It  should  also  be  stated  that  the  equivalence  in  Routine  W  between  the  entry 

cycle  code  (ENTCYCL)  and  the  working  cycle  length  code  (CYCL)  which 

existed  in  the  first  generation  software  was  removed  in  the  second  generation 

version.     This  was  necessary  because  the  entry  cycle  length  code  may  be 

manipulated  twice  during  the  process  of  arriving  at  the  final  optimized  traffic 

signal  timing  plan  computed  by  the  second  generation  software.     The  first  is 

when  the  Time-of-Day  plan  is  read  in  from  the  RAD  to  obtain  the  fixed 

interval  durations.    The  second  is  after  the  LSTSQS  Routine  computes  the 
optimal  network  cycle  lengths.     If  this  process  were  carried  out  with  the 

entry  and  working  cycle  length  codes  equivalenced,   the  controller  do  loop 

routines  operating  at  Priority  Level  2  would  cause  the  Offset  Correction, 

Routine  U,   to  be  activated  when  it  was  not  desired  to  do  so. 


2.6.80.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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RESET  DISC 
CALL  FLAG; 
I.E.  DCF  =  0 


0 


Routine  W  (Continued) 
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KK  =  JJ*NPATB 


YES 


NN  =  0 

TSTN  =  X'7FFFFFFF' 

LSIM=  LNKSKJJ+1) 


11=  LNKSKJJ) 
TEST(NN)  =  0 


3 


=  SCTLN(II) 


TEST(NN)  =  TEST(NN)  + 
|HSBUF(NPTB4+NN+4*I) 
-  VOL(l)-KOCC*AHCC(l)| 


I 


11=  II +  1 


Routine  W  (Continued) 
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Routine  W  (Continued) 


32 


Routine  W  (Continued) 
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© 


ENABLE  LEVEL  2  INTERRUPT  (X'63') 
Ll=  LNKSKJJ),  LSIM=  LNKSKJJ+1) 


JL  =  SCTLN,  VHIST(JL)  =  LEFT  HALF  OF 
PTN(J),  AHIST(JL)  =  RIGHT  HALF  OF  PTN(J) 


JI  =  JI+  1 
Ll=  LI  +  1 


YES 


YES 


SET 

AUTOMATIC  CONTROL  DESIRED 

FLAGS  (AUTOF)  FOR  ALL 

CONTROLLERS  IN 

SECTION  (JJ)  WHICH 

DO  NOT  HAVE 

CONTROLLER  MALFUNCTION  (CMALF) 

FLAGS  SET 


SET  TRANSITION 
FLAGS  (TRNSF) 
FOR  ALL 

CONTROLLERS  IN 
SECTION  (JJ) 


6 


Routine  W  (Continued) 
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DISABLE  LEVEL  2 
INTERRUPT  (X'63') 

KK   =    BYTEO  OFPTN(JI) 
NSL  =    BYTE  1  OF  PTN(JI) 
ENOFF  (KK)   =    BYTE  2  OF  PTN(JI) 
CYCLC  (KK)   =    BYTE  3  OF  PTN(JI) 
SSSS  (KK)   =    PTNUI+1) 
GGGG  (KK)   =    PTNUI+2) 
GGGT  (KK)   =    PTNUI+3) 
TKM   =    BYTE  3  OF  PTN(JI+3) 
ENINT    (4*KK)   =    PTNUI+4) 
ENINT  (4*KK+1)  =    PTNUI+5) 
ENINT  (4*KK+2)   =    PTNUI+6) 
ENINT  (4*KK+3)   =    PTNUI+7) 
Jl   =   JI+8 


" 

GUABC  (KK)   =    PTN(JI) 
GLABC  (KK)   =    PTNUI+1) 
CTODF  (KK)   =    BYTE  3  OF 

PTN(JI) 
Jl   =   JI+  5 
JB=    4  *J| 

Jl  =  JI+  2 


Routine  W  (Concluded) 
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2.6.  81    SUBROUTINE  INTRP 

2.6.81.1  Purpose 

Subroutine  INTRP  enables  or  disables  the  Priority  Level  2  interrupt. 

2.6.81.2  Relationship 

This  subroutine  is  called  from  routines  K  and  W  by: 

CALL  INTRP  (I) 

where: 

I  is  the  mnemonic: 

I  =  0  to  disable  interrupt 
I  =  1  to  enable  interrupt 

2.6.81.3  Assumption 
None. 

2.6.81.4  Approach 

This  subroutine  is  called  from  Routines  K  and  W  in  Priority  Level  3.     It  is 
therefore  not  required  to  be  re-entrant. 

2.6.81.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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YES 


DISABLE 
PRIORITY 
LEVEL  2 
INTERRUPT 


ENABLE 
PRIORITY 
LEVEL  2 
INTERRUPT 


c 


RETURN 


> 


Subroutine  INTRP 
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2.6.82    SUBROUTINE  PATGENT 

2.6.82.1  Purpose 

Subroutine  PATGENT  computes  various  parameters  for  routine  optimi- 
zation. 

2.6.82.2  Relationship 

PATGENT  is  called  from  Routine  W  and  uses  Subroutine  SMINT. 

2.6.82.3  Assumption 
None. 

2.6.82.4  Approach 

Subroutine  PATGENT  uses  the  phase  definition  (starting  and  ending  inter- 
vals) and  the  signal  timing  pattern  interval  durations  and  offsets  for  each 
intersection  to  compute  the  time  from  the  start  of  the  first  interval  to  the 
start  of  A-phase,    the  duration  of  the  overlap,   the  offset  to  the  beginning  of 
each  phase,    and  the  lost  time  (all-red  time)  in  each  phase). 

2.6.82.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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C   ENTER  J 


!NTERVH=  INTHKINT) 
IP0SN  =  4*(INT-1)  +  1 
ASTRT  =  ABCSTRT  (IPOSN) 
AEND  =  ABCENT  (IPOSN) 
VARA=  VARINT  (IPOSN) 
VARA  =  VARA  +  1 


IP0SN  =  4*(INT-l)+2 
BSTRT  =  ABCSTRT  (IPOSN) 
BEND -ABCENT  (IPOSN) 
VARB=  VARINT  (IPOSN) 
VARB=  VARB  +  1 
CLC=  ENTCYCLC  (INT) 


I 


CY=  KCL(CLC  +  1) 
SUMTA=0 
SUMTB=  0 
SUMTC  =  0 
SUMT=  0 
SUMAC  =  0 
DEL=  0 


I 


IF0T=  0 
S2T=0 


I 


CALL  SMINT 


ATIME  =  DURATION  OF  A  PHASE 


CALL  SMINT 


BTIME=  DURATION  OF  B  PHASE 


Subroutine  PAT  GENT 
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IPOSN  =  16*  (INT-1)  +  VARA 
VARAT=  ENTRYINT  (IPOSN) 
SA=ATIME  -  VARAT 
IPOSN  =  16*  (INT-1)  +  VARB 
VARBT  =  ENTRYINT  (IPOSN) 
SB=  BTIME  -VARBT 
Bl  =  0 


IF  A  PHASE  DOES  NOT 
START  AT  INTERVAL  1 
THEN  Bl=  TIME  FROM 


AENDP1=  AEND  +  1 


AENDP1  J^ 

sJMO 

AENDP1  = 

INTERVH? 

AENDP1 -  INTERVH 

[YES 
L^ 

BENDP1=  BEND  +  1 


BENDP1  2N 

sNO 

BENDP1  = 

INTERVH? 

BENDP1 -  INTERVH 

yYES 
L« 

ASM1  =  ASTRT-  1 


Subroutine  PATGENT  (Continued) 
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BSM1  =  BSTRT-1 


BSM 


\     NO 

1  -  1 '   >                      to 

BSM1  =  BSM1  +  INTERVH 

-/YES 

Lm> 

K3PHT=K3PH(INT) 


K3PHT  =  0  ? 


TYPE  0  CONTROLLER 


AENDP1=^ 

vNO 

BSTRT?     ., 

CALL  bMiN  i 

^T^ES 

l^ 

THIS  CALL  TO  SMINT  SETS 
SUMTA  EQUAL  TO  LOST 
TIME  BETWEEN  A  AND  B 
PHASES 


BEND  PI  =^ 

vNO 

ASTRT? y 

CALL  bMlIM  1 

YYES 

Urn, 

SUMT=  CY-ATIME-BTIME 
ABCTIME=  CY-SUMT 
1000 


ST=ATIME* 

VARC  =  0 

SC=0 

IF0T=0 


ABCTIME 


"© 


SUMTB  =  LOST  TIME 
BETWEEN  BAND  A 
PHASES 


Subroutine  PATGENT  (Continued) 
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IP0SN  =  4*(INST-l)  +  3 
CSTRT  =  ABCSTRT(IPOSN) 
CEND  =  ABCEND  (1POSN) 
CENDP1=  CEND  +  1 


CENDP1  <;^ 
INTERVM?^ 

^NO  . 

— 

CENDP1  =  CENDP1 
-  INTERVH 

^ms 

CSM1=  CSTRT-1 


CSM1 


NO 

CSM1  = 

CSM1  + 

^  1?  . 

INTERVH 

YES 

CALL  SMINT 


CTIME  =  THE  TIME  OF  C  PHASE 


K3PHT  = 


\P°) 


K3PHT  =  2 


K3PHT  =  3 


■© 


€) 


Subroutine  PATGENT  (Continued) 
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SUMT=  CY-  ATIME 
-  BTIME-CTIME 
ABCTIME=  CY-  SUMT 

ST  =  ATIME*  1000 

ABC  TIME 
..-,-r     CTIME*1000 
^  '  ABC  TIME 


I 


IP0SN  =  4*(INT-l)  +  3 
VARC=  VARINT(IPOSN) 
VARC  =  VARC  +  1 
IPOSN  =  16*<INT-1)  +  VARC 
VARCT=  ENTRYINT(IPOSN) 
SC=  CTIME- VARCT 




IFOT=0 


VARB  >  \     YES 
VARC? 


<B 


AENDPl \ 
=  BSTRT?  > 

NO 

- 

... 

CALL  SMINT 

i 

^fYES 
m 

JSUMTA  =  LOST  TIME  FROM  C  TO  B 


BENDP1  \ 
CSTRT?^ 

.NO 

- 

CALL  SMINT 

^""yes 

La 

SUMTC  =  LOST  TIME  FROM  A  TO  C 


CENDP1 
= ASTRT? 


NO 


CALL  SMINT 


ISUMTB  =  LOST  TIME  FROM  B  TO  A 


Subroutine  PATGENT  (Continued) 
43 


AENDP1  = 
CSTRT? 


YES 


SUMTA=  LOST  TIME 
BETWEEN  C  TO  B  PHASES 


CENDP1=\ 
BSTRT?  y 

.NO 

-* 

CALL  SMINT 

N^YES 

SUMTC=  LOST  TIME 
BETWEEN  A  TO  C  PHASES 


BENDPl^ 
ASTRT?    > 

OMO 

— 

CALL  SMINT 

Ssl^ES 

SUMTB=  LOST  TIME 
BETWEEN  B  TO  A  PHASES 


Subroutine  PATGENT  (Continued) 
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STRTAC  =  ASTRT 
STACM1  =  ASM1 


YES /ASTRT 
CSTRT? 


YES 


SJRTAC=  CSTRT 
STACM1  =  CSM1 


J 


CALL  SMINT 


CALL  SMINT 


I 


DEL  =  TIME  FROM  C 
PHASE  START  TO  A 
PHASE  START 


DEL=  DEL  +  128 


I 


-H*- 


_  JDEL  =  TIME  FROM  A  PHASE 
START  TO  C  PHASE  START 


ENDAC= AEND 
ENACP1  =  AENDP1 


© 


Subroutine  PATGENT  (Continued) 
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ENDAC=  CEND 
ENACP1=  CENDP1 


ATIME=TIMEOF 
--  A/C  PHASE 


ENACP1=N 

.NO 

CALL  SMINT 

BSTRT?    > 

I^/ES 

SUMTA  =  LOST  TIME  FROM  THE 
A/C  PHASE  TO  THE  B  PHASE 


CALL  SMINT 


SUMTB=  LOST  TIME  FROM  THE 
B  PHASE  TO  THE  A/C  PHASE 


SUMT=  CY-ATIME 

-BTIME 

ABCTIME=  CY-SUMT 

SUMAC  =  ATIME-YARAT 

SC=  CTIME-VARAT 

ZJ     ATIME*10QO 

bl     ABC  TIME 


Subroutine  PATGENT  (Continued) 
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<=> 


STRTBC=  BSTRT 
STBCM1=  BSM1 


BSTRT  =  \YES 
CSTRT? 


YES 


CALL  SMINT 


STRTBC=  CSTRT 
STBCM1=  BSM1 


I 


CALL  SMINT 


I 


STRTBC  =  CSTRT 
STBCM1=  CSM1 


I 


DEL  =  DEL+  128 


DEL  =  TIME  FROM  C- 
PHASE  START  TO  B- 
PHASE  START 


s 


DEL=  TIME  FROM 
B-PHASE  START  TO 
C-PHASE  START. 


Subroutine  PATGENT  (Continued) 
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ENDBC  =  BEND 
ENBCP1  =  BENDP1 


YES 


ENDBC  =  BEND 
ENBCP1  =  CENDP1 


BTIME=  THE  TIME 
OF  THE  B/C  PHASE 


SUMTA  =  LOST  TIME  BETWEEN 
A  PHASE  TO  B/C  PHASE 


CALL  SMINT 


SUMTB=  LOST  TIME  BETWEEN 
B/C  PHASE  TO  A  PHASE 


SUMT=  CY-ATIME-BTIME 
ABCTIME=  CY-SUMT 
SUMAC  =  BTIME-VARBT 
SC=  CTIME-VARBT 
=T     ATIME*1Q0Q 
51         ABCTIME 


© 


Subroutine  PATGENT  (Continued) 
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© 


AOFF  =  ENTRYOFF(INT) 
AOFF  =  AOFF-BI 


CALL  SMINT 


I 


TBOFF=TIME  FROM 
A-PHASE  TO  B-PHASE 


BOFF=  BOFF  +  AOFF 


COFF=  0 


Subroutine  PATGENT  (Continued) 
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yxK3PHT=  0?\  N0 

jTyes 

CALL  SMINT 

— . 

TCOFF=  TIME  FROM 
A-PHASE  TO  C-PHASE 

* 

COFF=  TCOFF  +  AOFF 

s^                     ^v      N 

<    COFF  <  CY    >► 

0 

COFF=COFF-CY 

L^ES 

^ 

\ 

1 

;3PHT^ 

0  .OR. 

LFOT=  0? 


YES 


NO 


BOFF  =  AOFF 


CALL  SMINT 


TAOFF=  TIME  FROM 
B-PHASE  TO  A-PHASE 


Subroutine  PATGENT  (Continued) 
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BIAS(INT)  =  Bl 
DELT(INT)=  DEL 
SUMACBC(INT)  =  SUMAC 
IPOSN=  2*INT-1 
SPLIT(IPOSN)  =  ST 
IPOSN=  2*INT 


I 


SPLIT(IPOSN)=  S2T 
IP0SN  =  4*INT-3 
IO(IPOSN)=  AOFF 
SMTABC(IPOSN)  =  SUMTA 
SABC(IPOSN)=  SA 
IP0SN=4*INT-2 


I 


IO(IPOSN)=  BOFF 
SMTABC(IPOSN)=  SUMTB 
SABC(IPOSN)=  SB 
IP0SN  =  4*INT-1 
SMTABC(IPOSN)=  SUMTC 
SABC(IPOSN)  =  SC 


I 


IPOSN  =  4*INT 
IO(IPOSN)=  IFOT 
SMTABC(IPOSN)  =  SUMT 


^RETURN  J 


Subroutine  PATGENT  (Concluded) 
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2.6.83    SUBROUTINE  SMINT 

2.6.83.1  Purpose 

This  subroutine  sums  the  interval  duration  (ISM)  for  controller  IC  from 
interval  KL  to  interval  KH. 

2.6.83.2  Relationship 

Subroutine  SMINT  is  called  from  Subroutine   PATGENT  by: 
CALL  SMINT  (IC,    KL,    KH,    ISM) 

where 

IC  specifies  the  controller 

KL  specifies  the  first  interval 

KH  specifies  the  last  interval 

ISM    returns  the  sum  of  the  interval  durations  from  interval 
KL  to  interval  KH. 

2.6.83.3  Assumption 
None. 

2.6.83.4  Approach 


The  duration  of  the  intervals  are  summed.     Provisions  for  cycle  warp- 
around  (i.e.,    interval  KL  occurs  after  interval  KH)  are  provided. 


2.6.83.5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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Q  ENTER   J 


ISM=0 

!L=KL 

IH  =  KH 

'L      ^VYES 

0?        / 

IL=  IL  +  JH 

Yfto 
k 

IH    "Ns^YES 

0?       >"^ 

IH  -  Ih  +  jh 

YfJo 
M 

IL       >v    YES 


JH 


NO 


IL=  IL~  JH 


IH 

Jh 


YES 


IH=  IH  -  JH 


Subroutine  SMINT 
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1=  IL 

JC=  16*  (IC  -  1) 


Subroutine  SMINT  (Concluded) 
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2.6.84    SUBROUTINE  WCODE 

2.6.84.1  Purpose 

This  subroutine  temporarily  holds  arguments.     It  contains  no  executable 
code. 

2.6.84.2  Relationship 

This  is  a  storage  area.     This  subroutine  is  not  called. 

2.6.84.3  Assumption 
None. 

2.6.84.4  Approach 
N/A 

2.6.84.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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c 


ENTER 


3 


ALOTTED  STORAGE  FOR 

WCDDE 

WTEST 

IR3 


C 


RETURN 


j 


Subroutine  WCODE 
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2.6.85    SUBROUTINE  WSUB 

2.  6.85. 1    Purpose 

Subroutine  WSUB  reads  in  the  proper  pattern  (operator  chosen  or  time -of- 
day)  from  the  RAD. 

2.6.85.2  Relationship 

This  subroutine  is  called  from  Routine  W  by: 
CALL  WSUB 

2.6.85.3  Assumption 
None. 


2.6.85.4    Approach 

This  subroutine  is  called  from  Routine  W  and  is  therefore  not  required 
to  be  re-entrant. 

The  flag  WCODE  indicates  whether  a  pattern  or  a  history  selection  should 
be  read  from  the  RAD.     If  a  RAD  read  error  is  encountered,  the  flag  WTEST 
is  set  negative. 


2.6.85.5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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C     ENTER  J 


ERROR  FLAG 
WTEST=  0 


DOES  \  YES 
MCODE=  1? 


LOOK  UP  DISC  STARTING  LOCATIOf\ 
FOR  DESIRED  PATTERN  NUMBERUJ 


START  READING  IN 
DESIRED  PATTERN 
(JJ)  FROM  DISC 


YES 


DSPP=  DISKSP 
READ  IN  HISTORY  SELECTION 
FROM  RAD  INTO  HSBUF 
STARTING  POINT  =  DSPP 
#  OF  WORDS  =  HBUFLNG 


ERROR  FLAG 
WTEST=  -1 


Subroutine  WSUB 
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2.6.86    ROUTINE  X 

2.  6.86.  1    Purpose 

Semi-actuated  controllers  normally  remain  in  A-Phase  Green,   and  go 
through  a  B-Phase  vehicle  or  pedestrian  sequence  (or  a  combination  there- 
of) only  when  a  vehicle  or  a  pedestrian  demand  is  present  at  a  specified 
time  within  the  controller  cycle.     The  proposed  UTCS  computer  program 
will  transmit  a  yield  pulse  to  the  semi-actuated  controller  on  a  once-per- 
cycle  basis.     The  controller  will  use  this  pulse  to  initiate  one  of  the 
following  B-Phase  sequences  as  a  function  of  the  existing  demand  at  the 
intersection: 

1.  Pedestrian  sequence 

2.  Vehicle  sequence 

3.  Combination  of  the  above  two  sequences 

4.  No  B-Phase  sequence. 

For  (1),   (2),   or  (3)  above,  the  local  controller  will  initiate  and  time  the 
A-Phase  clearance  intervals  and  the  required  B-Phase  vehicle  or  pedestrian 
intervals.     The  total  time  required  for  the  sum  of  these  intervals  will  vary 
as  a  function  of  the  demand,   but  the  requested  sequence  will  always  be  com- 
pleted within  a  prescribed  maximum  time  duration  after  the  yield  pulse. 

The  function  of  the  semi-actuated  part  of  the  UTCS  computer  program  is  to 
generate  the  yield  pulse  which  will  meet  the  following  requirements: 

1.      Maintain  the  same  cycle  length  at  the  semi-actuated  con- 
troller as  that  which  is  in  effect  at  the  adjacent  controllers 
in  the  section. 
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2.      Transmit  the  yield  pulse  at  such  a  time  that  the  subse- 
quent A-Phase  clearance  intervals  and  the  longest 

B-Phase  intervals  will  cause  the  next  sequential  begin- 
ning of  A-Phase  Green  to  occur  at  the  desired  main 
street  offset  time. 

The  procedure  of  (2)  will  maintain  the  desired  main  street  vehicle  pro- 
gression speeds.     The  start  of  A-Phase  Green  may  occur  earlier  than 
the  offset  time  for  the  case  of  shorter  B-Phase  intervals,  but  it  will  never 
occur  later  than  this  time. 

The  specific  functions  performed  by  Routine  X,   Semi-actuated  Controller 
Routine  are  to: 

1.  Transfer  the  semi-actuated  controller  from  standby  to 
computer  control  smoothly. 

2.  Time  the  intervals  required  to  transmit  the  yield  pulse 
to  the  controller. 

3.  Perform  controller  state  error  detection  to  verify  that 
the  controller  is  returning  to  A-Phase  Green. 

4.  Maintain  the  required  offset  relationship  with  respect 
to  the  adjacent  controllers  in  the  section. 

5.  Provide  for  transition  from  one  cycle  length  and  off- 
set to  another. 

6.  Transition  the  semi-actuated  controllers  from  com- 
puter control  to  standby. 
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2.6.86.2    Relationship 

This  routine  is  part  of  the  controller  loop  Level  2  Routines.     The  routine 
is  executed  following  the  completion  of  Routine  V,   Critical  Intersection, 
which  is  the  last  routine  in  the  Do-Loop  for  traffic-adjusted  controllers. 
The  routine  is  executed  once  per  second  for  each  semi-actuated  controller 
in  the  system. 


2.6.86.3  Assumption 
None. 

2.6.86.4  Approach 

Routine  X  will  perform  the  following  series  of  computations.     (Refer 
to  the  Routine  X  Flow  Charts  for  a  detailed  functional  flow  of  these 
computations, 

1.  The  semiactuated  controller  will  be  picked  up  from 
standby  when  the  Automatic  Control  Desired  Flag 
(AUTOF)  is  set,  and  the  controller  is  in  A-Phase 
Green.     The  A -Phase  Green  will  be  extended  by  the 
amount  required  to  synchronize  the  controller  to 
the  computer-selected  cycle  length  and  offset,   Entry 
Cycle  Length  Code  (CYCLC)  and  Entry  Offset  (ENOFF), 
respectively. 

2.  The  computer  will  time  the  following  two  intervals: 

a)      The  first  interval,  GA,   will  extend  from  the  desired 
offset  time  to  the  beginning  of  the  first  A-Phase 
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clearance  interval,  at  which  time  the  "Yield  Pulse" 
will  be  set.     The  controller  will  then  initiate  the 
A-Phase  clearance  intervals  if  a  B-Phase  vehicle 
or/and  pedestrian  demand  is  present. 

b)      The  second  computer-timed  interval,   GB,  will 
extend  from  the  beginning  of  A-Phase  pedestrian 
clearance  to  the  desired  offset  time.     It  should  be 
noted  that  this  computer-timed  interval  corres- 
ponds to  the  summation  of  the  following  intervals 
at  the  semi -actuated  controller: 

•  A-Phase  clearance  intervals. 

•  Maximum  B-Phase  vehicle  or  pedestrian  passage 
intervals. 

•  B-Phase  clearance  intervals. 

3.  The  computer  program  will  perform  controller  state  error 
detection  by  checking  for  the  presence  of  A-Phase  Green 
during  the  1  second  time  period  prior  to  the  normal 
transmission  of  the  "Yield  Pulse.  "    If  this  return  is  not 
present,  the  program  will  take  the  following  steps: 

a)  Place  the  controller  into  standby  (HOLD  =  0). 

b)  Set  the  "Controller  Malfunction  Flag.  " 

c)  Call  the  CRT  Failure  Table  Add  Subroutine,   which 
is  part  of  Routine  K,   CRT  Control  and  Data  For- 
matting. 

4.  The  computer  will  compare  the  desired  offset,   Working 
Offset  (WRKOF),  with  the  measured  offset,  Tck,  at  the 
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beginning  of  the  first  computer -timed  interval.     If  an  ad- 
justment is  required  due  to  either  an  operator  offset  change 
or  an  offset  error,   the  program  will  adjust  the  offset  of  the 
first  interval  by  the  required  amount  subject  to  the  con- 
straints that  the  change  be  greater  than  -4  seconds  (i.e., 
4  seconds  earlier)  and  less  than  cycle  minus  4  seconds. 
It  should  be  noted  that  these  restrictions  also  apply  to  a 
transition,    which  is  described  below. 

5.      The  program  will  provide  a  smooth  transition  from  one 
selection  of  cycle  length  and  offset  to  another  at  the 
start  of  the  first  programmed  interval.     At  this  time, 
the  following  data  transfers  will  take  place: 


• 


• 


Entry  QA  to  GA 
Entry  to  GB 


•  "Entry  Cycle  Length  Code"  to  "Working  Cycle  Length 
Code." 

•  "Entry  Offset"  to  "Working  Offset.  " 

The  actual  transition  will  be  effected  by  either  decreasing 
(up  to  4  seconds)  or  increasing  the  first  computer-timed 
interval. 

6.      The  semi-actuated  controller  will  be  transferred  from 
computer  control  back  to  standby  at  the  end  of  the  first 
computer-programmed  interval,  GA,  when  the  "Standby 
Desired  Flag"  is  set.     This  technique  will  maintain  at 
least  as  much  A-Phase  Green  at  the  controller  as  that 
which  existed  under  computer  control.     The  initial 
A-Phase  Green  may  be  longer  than  this  duration  by  as 
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much  as  one  controller  standby  cycle  length  when 
there  are  vehicle  and/or  pedestrian  demands  present 
at  the  intersection.     The  A-Phase  Green  may  extend 
indefinitely  if  there  is  no  cross  street  demand. 

Routine  X  is  called  by  the  Level  2  Executive  Routine,   Routine  Al.     Routine 
Al  calls  Routine  X  if  there  is  at  least  one  semi-actuated  controller  in  the 
group  of  32  controllers  being  processed.     Upon  entry,   Routine  X  examines 
each  of  the  32  controllers,   with  array  SEMSK,  to  determine  if  it  is  a 
semi-actuated  type.     When  a  semi-actuated  controller  is  located,  array 
HOLD  is  examined  to  see  if  the  controller  is  in  standby.     A  semi-actuated 
controller  in  standby  is  brought  on  line  if  the  automatic  control  desired 
flag,  AUTOF,   and  corrected  A-Phase  Green  Return,   AOG,  are  set. 

After  a  semi-actuated  controller  has  been  brought  on  line,   AUTOF  is  reset 
and  the  countdown  register,   CNTDR,   is  counted  down -from  its  initial  value 
of  ENINT  (1+16*CNTRN)  which  contains  duration  entry  GB.     Entry  GB 
extends  from  the  beginning  of  A-Phase  pedestrian  clearance  up  to  the 
desired  offset  time.     This  interval  is  composed  of  the  A-Phase  clearance 
intervals,  maximum  B-Phase  vehicle  or  pedestrian  passage  intervals,   and 
B-Phase  clearance  intervals. 

After  CNTDR  has  counted  down  to  zero,  the  interval  number,  INTNM,    is 
set  to  1,   and  an  adjustment  is  made  for  the  required  offset  relationship 
with  respect  to  the  adjacent  controllers  in  the  section.     The  desired  offset, 
WRKOF,   is  compared  with  the  measured  offset,   KCL-KCTC,  to  form 
deviation  EPS.     If  an  adjustment  is  required  (EPS/0)  due  to  either  an 
operator  offset  change  or  an  offset  error,   the  program  will  adjust  the  first 
interval  by  the  required  amount  subject  to  the  constraints  that  the  change 
be  greater  than  -4  seconds  and  less  than  the  cycle  length  minus  4  seconds. 
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If  the  transition  flag,  TRNSF,   is  on,   a  smooth  transition  is  provided  from 
one  selection  of  cycle  length  and  offset  to  another  at  the  start  of  the  first 
interval,  by  transfer  of  appropriate  entry  data  to  working  storage.     These 
operations  are  performed  prior  to  the  offset  corrections  described  in  the 
previous  paragraph. 

After  the  offset  corrections  are  made,  the  countdown  register  is  set  to  entry 
GA  plus  the  final  offset  correction,    if  any.     Entry  GA  extends  from  the 
desired  offset  time  to  the  beginning  of  the  first  A-Phase  clearance  interval. 
After  CNTDR  has  counted  down,   INTNUM  is  set  to  2  and  AOG  is  examined. 
If  AOG  is  zero,  the  controller  is  assumed  to  have  malfunctioned.     At  this 
time  the  controller  malfunction  flag,   CMALF  is  set  and  the  controller  is 
dropped  into  standby  by  setting  HOLD  to  zero.     After  this  operation,   the 
ADD  FAILURE  UPDATING  SUBROUTINE  is  called  to  record  this  failure. 
If  AOG  is  a  "one"  and  the  standby  desired  flag,  STBYD,   is  set,  the  flag 
is  reset  and  the  controller  is  placed  in  standby.     If  AOG  is  set  and  STBYD 
is  not  set,   CNTDR  is  set  to  entry  GB  and  the  yield  pulse,   YIELD,   is  set. 

For  each  call  to  Routine  X,  the  above  operations  are  repeated  for  each  of 
the  eight  controllers  being  examined.  Upon  examining  all  32  controllers, 
the  routine  is  exited. 

Routine  X  has  been  modified  to  transfer  the  entry  cycle  length  code 
(ENTCYCL)  to  the  working  cycle  length  code  (CYCL)  immediately  prior  to 
transition.     This  transfer  is  necessary  because  ENTCYCL  is  not  equated 
to  CYCL  in  the  second  generation  software  for  reasons  described  in  Routine 
W. 

The  analyst  should  refer  to  Section  2.  25  of  Reference  2  to  trace  the  following 
discussion.     On  Figure  2.  25.  5-1  of  Page  2.  25.  7/2.  25-8  of  Reference  2,   code 
has  been  added  at  the  location  equivalent  to  that  of  the  lower  right  hand  pro- 
gram box  to  perform  the  cycle  length  transfer.     This  code  is  immediately 
after  the  transfer  of  the  entry  offset  to  the  working  offset. 

65 


The  cycle  length  code  transfer  also  occurs  in  the  right  most  program  box 
on  Figure  2.  25.  5-3  of  Page  2.  25-10  of  Reference  2.  The  transfer  occurs 
immediately  prior  to  resetting  the  transition  flag  (TRANSF)  to  zero. 


2.6.86.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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ENTER 
FROM  ROUTINE  Al 


]=  1 

IE11  =  16*(1C0N+1)  -15 
IE12  =  IE11+1 
ISHF=  ICON+1 


YES 


NO 


AUTOF  (l)=  0 


YES 


TRNSF  (l)=  0 
STBYD  =  0 
FSETCF(I)=  0 


CNTDR  (ISHF) 
CNTDR  (ISHF)-  IDECR 


YES 


INTNM(ISHF)=  1 


YES 


YES 


INTNM(ISHF)   =  2 

INTDU(IEIl)     =  ENINT(IEll) 

INTDUCIE12)   =  ENTNT(IE12) 

CNTDR(ISHF)  =  INTDUCIE12) 

HOLD(l)  =  1 

WRKOF(I)        =  ENOFF(I) 


©  © 


Routine  X 
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YES 


FSTCF(I)=  0 
INDEX  =  CYCLE(ISHF) 
CUROF(l)=  KCTC(INDEX) 
EPS  =  WRKOFU)-KCTC(INDEX) 


INTDU(IEI1)=  ENINT(IEIl) 
INTDU(IEI2)=  ENINTCIEI2) 
WRKOF(l)  =  ENOFF(J) 
TRNSF(I)=  0 


Routine  X  (Continued) 
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CMALF=  1 
HOLD  (I)  =  0 


SET  BIT  0    . 

IN  MAST  (KMAST) 


CONTB  (CLENG) 
=  ISHF 


CLENG = 
CLENG+1 


YES 


YES 


CNTDR  (ISHF) 
INTDU  (IE12) 
YIELD  (l)=  1 


fr«- 


© 


STBYD  (I)  =  0 
HOLD  (l)=  0 


Routine  X  (Continued) 
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NO 


YES 


EPS  = 

EPS-KCL  (INDEX) 


EPS  = 

EPS+KCL  (INDEX) 


CNTDR  (ISHF) 
INTDU  (IE12) 
+  EPS 


0  0 


CNTDR  (ISHF) 
INTDU  (IE12) 


Routine  X  (Concluded) 
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2.6.87    ROUTINE  Y 

2.  6.87.  1    Purpose 

Routine  Y  is  entered  at  the  start  of  real  time  operations.     After  the  program 
is  read  in  from  the  RAD,  this  routine  requests  the  CRT  System  Status  page 
and  initializes  the  displayed  data.     The  Startup  Routine  contains  data  from 
constants,   fixed  and  variable  arrays,   and  assigns  reserve  areas  for 
dynamic  storage.     The  constants  and  arrays  are  initialized  by  the  assem- 
bler and  the  dynamic  storage  area  is  initialized,  to  zero,  by  this  routine. 


2.6.87.2    Relationship 

The  routine  is  only  entered  when  system  operation  under  computer  control  is 
initiated.     The  startup  procedure  is  requested  by  the  operator  through  a 
monitor  key-in.     Upon  completion  it  exits  to  the  level  initialization  loop. 


2.6.87.3  Assumption 

The  shutdown  file,    which  may  be  used  at  startup,    has  not  been  disturbed  and 
remains  correct  from  the  previous  shutdown. 

2.6.87.4  Approach 

When  the  program  is  read  in  from  the  RAD,    at  the  operator's  request,    it 
initializes  all  dynamic  arrays  to  zero.     The  arrays  stored  at  shutdown  are  read 
in  from  the  RAD  if  sense  switch  4  is  in  the  down  position.     The  routine  then 
initializes  the  TOD  clocks  as  per  the  operator  inputted  TOD.     It  then  transfers 
directly  to  the  Level  1  executive  initialization  loop  of  Routine  A  which  sets  up 
the  priority  structure  of  the  program  and  begins  the  real  time  processing. 
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Routine  Y  was  modified  to  initialize  the  FORTRAN  library  upon  entering 
the  routine.     This  initialization  is  accomplished  through  a  call  the  system 
routine,   9INITIAL.     This  initialization  is  required  to  support  the  FORTRAN 
routines  of  the  second  generation  software. 

Routine  Y  was  also  modified  to  initialize  the  Priority  Level  5  software 
counter,    5MINCT.     The  assumption  used  in  initalizing  this  counter  is  that 
it  is  zero  at  midnight  of  the  previous  day.    The  counter  is  used  by  the 
system  executive  Routine  (RTNA)  to  schedule  the  Priority  Level  5  loop 
when  5  MINCT  goes  to  zero.     The  user  may  specify  the  scheduling  frequency 
of  the  Level  5  software  through  merely  changing  one  parameter  in  the  data 
base,   OPTIME.     OPTIME  is  the  number  of  minutes  desired  between  con- 
secutive calls  by  the  system  executive,   RTNA  to  the  Priority  Level  5 
executive,   MACRO.     The  surveillance  constants  in  RTN5,  which  are  time 
dependent,  are  computed  as  a  function  of  the  value  of  OPTIME. 


2.6.87.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 


72 


c 


ENTER  FROM  RBM 

I 


") 


ZERO  OUT  ALL 
DYNAMIC  ARRAYS 


YES 


READ  THE  SHUTDOWN 
PARAMETERS  FROM 
RAD  AND  RESTORE 


TRANSFER  THE  BPSCFF 
AND  BPSDFF  ARRAYS  TO 
BPSCFFLAND  BPSDFFL, 
RESPECTIVELY 


I 


CRTDC1  =  2,  CRTDC2  =  2 
CRTFL1  =  1,  CRTFL2  =  1 
SCHEDULE  ROUTINE  K  IN  MAST 
GET  RBM'S  TFME  AND  CONVERT 
TO  MINUTES  SINCE  MIDNIGHT^) 


FMTODC=  (K/15) 
FMC=K/15  -  (K/15) 
I.E.FMC=K(MOD  15) 


DETERMINE  THE 
DAY  OF  WEEK 
AND  STORE  IN 
DAYOWK 


(EXIT  TO  LEVLL  :  EXECUTIVE'S  'N 
INITIALIZATION  LOOP  J 


Routine  Y 
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2.6.88     ROUTINE  MACRO 

2.  6.88. 1    Purpose 

Routine  MACRO  serves  as  the  executive  for  the  TANSTP  MACRO  loop.     It 
calls  the  Level  5  Surveillance  Routine,  the  Predictor  Routine,  the  Network 
Optimization  Routine,  the  Subnetwork  Interfacing  Routine,   and  the  Transition 
Optimization  Routine. 

2.6.88.2  Relationship 

Routine  MACRO  is  called  from  the  Executive  Routine  A  upon  triggering  the 
hardware  interrupt  X'69'.     This  interrupt  is  at  Priority  Level  5  in  the  UTCS/ 
2GCF  system  and  is  triggered  every  5MINCT  seconds  in  Routine  A  after  the 
Priority  Levels  2  and  3  interrupts  have  been  triggered.     At  increments  of 
5MINCT  seconds,    the  Priority  Level  Five  loop  becomes  active.     Routine 
MACRO  calls  routines  FMODE,    INTFC,    IOCHK,    LSTSQS,    PREDICT,    RTSND, 
TRFLAG,    and  TRNPAR. 

2.6.88.3  Assumption 
None. 

2.6.88.4  Approach 

Upon  entry,   Routine  MACRO  sets  the  Level  5  Active  Flag  (LEVL5A)  to 
indicate  to  the  software  system  that  the  Level  5  loop  is  active.     A  check 
is  then  made  on  the  IPASS  flag  to  determine  if  a  system  shutdown  is  in 
progress.     If  so,  the  routine  is  exited.     If  the  IPASS  flag  is  not  set,  the 
reoptimization  required  flag,    IREOPT,   is  reset  and  a  call  is  made  to 
Subroutine  SHTDWN  to  determine  if  a  Level  5  shutdown  sequence  should 
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begin.     If  so,  the  IP  ASS  flag  is  set  for  subsequent  entries  to  the  MACRO 
Routine.     In  any  event,   the  subnetwork  interfacing  flag,   INTFLG,    is  reset 
and  a  call  is  made  to  Subrbutine  IOCHK.    This  subroutine  determines  if  any 
input  or  output  (I/O)  operations  are  in  progress.     If  so,  the  subroutine  is 
continually  recalled  until  the  I/O  has  been  completed. 

A  call  is  then  made  to  RTN5  to  compute  the  surveillance  parameters  of 
volume,   speed,  and  occupancy  and  to  write  the  Level  5  surveillance  tape. 
A  test  is  then  made  on  the  Level  5  first  pass  flag,  ICOUNT,  to  determine 
if  the  current  entry  to  MACRO  is  the  first.     If  so,   no  optimization  occurs 
due  to  lack  of  sufficient  surveillance  data  and  MACRO  is  exited. 

If  the  current  entry  is  other  than  the  first,  the  Predictor  Routine  is  called 
to  predict  the  values  of  volume  and  speed  during  the  next  optimization  time 
period.     If  the  Priority  Level  5  shutdown  procedure  is  to  be  initiated 
(IPASS^O),  the  MACRO  Routine  is  exited. 

If  shutdown  is  not  to  be  initiated,   a  call  is  made  to  Subroutine  FMODE  to 
determine  if  any  section  of  controllers  is  in  the  traffic  responsive  mode,     If 
not  (K=0),   MACRO  is  exited. 

If  any  section  is  in  the  traffic  responsive  mode,  the  Level  5  Start  Flags 
(LEV5SF)  are  tested  to  determine  which  sections  are  in  the  traffic  respon- 
sive mode.     The  first  generation  Routine  W  has  been  modified  to  set  and  re- 
set these  flags.     If  none  is  set,   program  control  is  returned  to  the  Monitor. 

If  any  of  the  LEV5SF  flags  is  set,    the  intersection  core  groups  that  should 
be  considered  for  combining  into  subnetworks  are  placed  in  the  LEV5S 
array  for  Routine  RTSND.     In  addition,    if  any  of  the  LEV5S  is  set  non  zero, 
a  test  is  made  to  determine  if  any  one  has  changed  since  the  last  5MINCT 
pass  through  the  program.     If  any  has  changed,    reoptimization  required  flag, 
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IREOPT,    is  set  to  indicate  that  different  sections  have  been  selected  for 
traffic  responsive  control  since  the  last  pass  and  hence  re  optimization 
of  the  network  should  occur. 

Routine  MACRO  then  checks  the  flag,  IRADRD,  to  determine  if  a  new  pattern 
has  been  read  into  the  entry  buffers  by  Routine  W.     (Routine  W  has  been 
modified  to  set  this  flag  under  these  conditions.  )    If  so  the  IREOPT  flag  is 
set  indicating  reoptimization  should  occur. 

Routine  RTSND  is  then  called  to  determine  subnetwork  configurations, 
intersection  splits,  and  cycle  lengths. 

Routine  MACRO  then  tests  the  flag  IREOPT  to  determine  if  reoptimization 
should  occur.     If  it  is  not  set,   program  control  is  returned  to  the  Monitor. 

If  reoptimization  is  to  occur  (IREOPT  =  1)  Routine  LSTSQS  is  called  to 
determine  the  offsets  and  cycle  lengths  which  minimize  delay  in  each  sub- 
network. 

A  test  is  then  made  on  the  subnetwork  interface  flag  (INTFLG)  to  determine 
if  all  sections  have  been  selected  for  traffic  responsive  control.     If  so,  the 
Routines  INTFC  and  TRNPAR  are  called  to  optimize  the  subnetwork  inter- 
faces and  to  optimize  the  transition,   respectively.     The  transition  flags  are 
then  set  by  Subroutine  TRFLG. 

In  any  event,  the  IRADRD  flag  is  reset,  the  ICOUNT  flag  is  set,  the 
LEVL5A  is  reset,  and  control  is  returned  to  the  Monitor. 


2.6.88.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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(ROUTINE 
MACRO  ENTER 
FROM  ROUTINE  A 


I 


SET  LEVEL  5  ACTIVE 

FLAG 

LEVEL5A=  1 


YES 


*© 


CALL  SHTDWN  SUBROUTINE 
TO  CHECK  FOR  SHUTDOWN  NOW 


CALL  SHTDWN  (SHTFLG) 


RESET  FLAG  FOR  NO 
INTERFACING 
INTFLG=  0 


Routine  MACRO 
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*»" 


CALL  IOCHK  SUBROUTINE 

TO  CHECK  FOR  INPUT/OUTPUT 


CALL  IOCHK  (^  0) 


YES 


CALL  RTN5  SUBROUTINE  TO 
COMPUTE  SURVEILLANCE 
DATA 


CALL  RTN5 


YES 


OVERLAY  LOAD  SEGMENT  3 


CALL  SEGLOAD  (3) 


I 


CALL  PREDICTOR 


CALL  PREDICT  (SHTFLG) 


■© 


Routine  MACRO  (Continued) 
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CALL  FMODE  TO  DETERMINE 
IF  ANY  SECTIONS  IN 
TRAFFIC  RESPONSIVE  MODE 


CALL  FMODE(K) 


ZERO  SECTION  COUNTER 
J=0 


I 


SET  FLAG  SO  THAT 
PREDICTED  SURVEILLANCE 
IS  USED 
IPRD=  1 


Routine  MACRO  (Continued) 
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TEST  LEVEL  5  START 
FLAGS  TO  DETERMINE 
SECTIONS  FOR  TRAFFIC 
RESPONSIVE 


UPDATE  SAVED  LEV5S 
LEV50(I)=  LEV5SU) 


Routine  MACRO  (Continued) 
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f 


■0 


YES 


SET  REOPTIMIZATION 
FLAG  IREOPT=  1 


YES 


SET  INTERFACE  FLAG 
INTFLG=  1 


OVERLAY  LOAD  SEGMENT  2 
CALL  SEGLOAD(2) 


Routine  MACRO  (Continued) 
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Routine  MACRO  (Continued) 
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0 


© 


CALL  SUBNETWORK 
INTERFACING 


CALL  INTFC 


I 


COMPUTE  TRANSITION 
PARAMETERS  ' 


CALL  TRNPAR 


I 


SET  TRANSITION  FLAGS 


CALL  TRFLG 


I 


RESET  RAD  READ  FLAG 
IRADED=  0 


3 


SET  FIRST  PASS  COMPLETE 

COMPLETE 

FLAG 

ICOUNT=  1 


I 


RESET  LEVEL  5  ACTIVE  FLAG 
LEVL5A=0 


I 


(RETURN  TO       } 
ROUTINE  A  J 


Routine  MACRO  (Concluded) 
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2.6.89     SUBROUTINE  FMODE 

2.6.89.1  Purpose 

Subroutine  FMODE  determines  if  any  section  of  controllers  is  in  the  Traffic 
Responsive  Mode  of  operation. 

2.6.89.2  Relationship 

This  subroutine  is  called  by: 

CALL  FMODE(K) 

whe  re : 

K  is  the  return  argument 

K=0  =  No  section  is  in  Traffic  Responsive  Mode 

K=l  =  At  least  one  section  is  in  the  Traffic  Responsive  Mode 

Subroutine  FMODE  is  called  by  MACRO  and  returns  to  that  routine. 

2.6.89.3  Assumption 
None. 

2.6.89.4  Approach 

The  mode  flags  (MODE)  for  all  sections  are  tested  for  a  value  of  3.     If  any 
flag  is  3,    K  is  set  1.     Otherwise  K  is  set  0. 

2.6.89.  5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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TRAFFIC  RESPONSIVE 
THE  DESIRED  MODE? 


Subroutine  FMODE 
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2.6.90    SUBROUTINE  OUTTBF 

2.6.90.1  Purpose 

Subroutine  OUTTBF  outputs  the  2GC  Tape  Buffer. 

2.6.90.2  Relationship 

Subroutine  OUTTBF  is  called  from  Routine  RTN5  by: 

CALL  OUTTBF  (TBUFF) 
where: 

TBUFF  specifies  the  tape  buffer. 

2.6.90.3  Assumption 
None. 

2.6.90.4  Approach 

OUTTBF  first  examines  sense  switch  2.     If  it  is  set,    OUTTBF  returns  to* 
the  calling  routine.     OUTTBF  then  addresses  the  proper  I/O  channel,    checks 
for  errors,    and  then  writes  the  array  TBUFF  to  the  tape.     OUTTBF  then 
returns  to  the  calling  routine. 

2.6.90.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 


86 


Q              ENTER     ^ 

ir 

EXIT  IF  SENSE 
SWITCH  2  SET 

SFT  UP  BUFFER 
ADDRESSING 

ir 

CHECK  FOR  ERRORS 

1 

f 

WRITE  ON  TAPE 

_3_ 

f 

r            RETURN     ^ 

Subroutine  OUTTBF 
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2.6.91    ROUTINE  RTN 5 

2.6.91.1  Purpose 

Routine  RTN5  is  the  second  generation  software  surveillance  data  processor. 
Its  purpose  is  to  compute  average  link  volumes,    average  link  speeds,    average 
link  occupancy,    free  flow  link  volume,    free  flow  link  speed,    free  flow  link 
occupancy,    and  link  surveillance  data.     These  data  are  stored  in  the  second 
generation  MOE  buffer,    TBUFF.     This  buffer  is  accessed  by  Subroutine 
LDMOE5  to  retrieve  input  data  for  the  level  5  optimization  routine,    LSTSQS. 
RTN5  is  executed  at  the  frequency  of  optimization  specified  by  the  user 
through  input  of  the  data  base  parameter,    OPTIME.     The  surveillance  data 
computed  by  RTN5  is  the  average  over  each  OPTIME  interval. 

2.6.91.2  Relationship 

Routine  RTN5  is  called  from  Routine  MACRO  by: 

CALL  RTN5 
The  Routine  calls  subroutines  OUTTBF  and  TODAY. 

2.6.91.3  Assumption 
None. 

2.6.91.4  Approach 

Routine  RTN5  first  loads  the  current  date  and  time  into  the  buffer,    TBUFF. 
Then  RTN5  loops  over  all  detectors  for  a  given  link  and  sums  each  element 
of  raw  detector  data  CC2P,    RAWSPP,    and  RAWOCC,    accumulated  by  Routine 
C  for  all  detectors  on  the  link.     Routine  RTN5  then  computes  the  average 
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link  volume,    average  link  speed,    average  link  occupancy,    and  loads  this  data 
into  the  buffer,    TBUFF.     Free  flow  data  is  obtained  in  much  the  same  manner 
from  the  detector  on  each  link  which  is  the  greatest  distance  from  the  stop 
bar.     RTN5  resets  the  raw  data  cccumulators  CC2P,    RAWSPP,    and 
RAWOCC. 

The  above  procedures  are  carried  out  for  each  link  in  the  UTCS  network. 
After  all  links  have  been  processed,    RTN5  calls  an  internal  subroutine, 
OUTAPE,    to  write  TBUFF  to  tape  on  tape  unit,    El. 

It  should  be  noted  that  if  a  detector  has  failed  on  a  link,    the  second  genera- 
tion link  fail  flag,    FAILP,    is  set  and  zeroes  are  loaded  into  TUBFF  for  the 
average  link  volume,    speed,    and  occupancy.     If  the  free  flow  detector  has 
not  failed  on  the  failed  link,    free  flow  data  will  be  calculated.     The  free 
flow  detector  is  the  upstream  detector  for  two  detector  links,   and  is  the 
detector  for  single  detector  links.     If  the  free  flow  detector  has  failed,   the 
free  flow  fail  flag,    DETM,    is  set  and  zeroes  are  loaded  into  TBUFF  for 
free  flow  volume,    occupancy,    and  speed. 


2.6.91.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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©■ 


CLEAR  TBUFF 
ARRAY-  ELEMENTS 
5  THROUGH  TBSIZE 
1=  NMVED  +  1 


-  .'/. 


±. 


IDET=0 
IRLNC=0 
IRLSP=0 
IRLNO=  0 
IOCCP=  0 
ISPEED  =  0 
ICOUNT=  0 


1=1-1 


YES 


CALL 

OUTTBF 

(TBUFF) 


-|/RETURN  ) 


LINK  =  DWL(I)  +  1 


Routine  RTN5 
90 


1  =  1  +  1 


IOC  =  (RAW0CP  +  8)/16 
IRLNC=  IRLNC  +  CC2P(I) 
IRLSP  =  IRLSP  +  RAWSPPd) 
IRLNO  =  IRLNO  +  IOC 


CLEAR 
CCZP(I) 
RAWSPP(I) 
RAWOCP(I) 


NDD=  NDD-  1 


YES 


Routine  RTN5  (Continued) 
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YES 

yOFML  (LINIOS 
>V       =0?    yS 

M 

NO 

IVOL  = 
(IRLNC*OPTIN) 

Nil 

IV0L  = 

VHIST  (LINK)/2 

ISP  = 

/URLSP*KV1)\*LL1 
\      IRLNC        / 

I0CC  = 

AHIST  (LINK)*100 
215 

i 

r 

_jr 

ISP  = 

(i$  *  (430) 

ISP  = 

VHIST  (LINK)  *215 

AHIST  (LINK)*  180 

1 


TBUFF(ISTX)  = 

=  LINK 

TBUFFUSTX  + 

1)  = 

IFAIL 

TBUFFUSTX  + 

2)  = 

IVOL 

TBUFF(ISTX  + 

3)  = 

ISP 

TBUFF  (ISTX  + 

4)  = 

I0CC 

TBUFF0STX  + 

5)  = 

IDF 

ISTX  =  ISTX  + 

10 

Routine  RTN5  (Concluded) 
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2.6.92    SUBROUTINE  TRFLG 

2.6.92.1  Purpose 

This  subroutine  interleaves  the  STB  YD  array  with  ones. 

2.6.92.2  Relationship 

Subroutine  TRFLG  is  called  from  Routine  MACRO  by: 

CALL  TRFLG 
This  subroutine  calls  subroutine  PTHLF. 

2.6.92.3  Assumption 
None. 

2.6.92.4  Approach 

All  bits  in  every  other  word  in  the  STBYD  array  are  set. 

2.6.92.  5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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Q  ENTER  J 


l  =  -l 

k 

1=1  +  2 

_3L- 

SET  ALL  BITS 
OF  STBYD  (1) 
TO  1. 

NO 

/     1  = 

=  7?      > 

YES 


RETURN 


Subroutine  TRFLG 
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2.6.93    ROUTINE  PREDICT 

2.6.93.1    Purpose 

The  objective  of  the  real-time  predictor  is  to  forecast  link  volumes  and 
speeds  for  the  5MINCT  period  ahead  of  the  current  time.     The  predictor  is 
presently  dimensioned  to  forecast  these  parameters  for  150  links. 


2.6.93.2    Relationship 

Routine  PREDICT  is  called  from  Routine  MACRO  by: 

CALL  PREDICT 

Routine  PREDICT  calls  subroutines  BUFFERIN,    BUFFEROUT,    LDMOES, 
and  PRPRED. 


2.6.93.3  Assumption 
None. 

2.6.93.4  Approach 

Routine  PREDICT  is  called  every  15  minutes  on  the  quarter  hour.     Upon 
entry  the  current  quarter  hour  time  is  extracted  from  location  HEAD21+8.     If 
the  predictor  was  called  at  the  previous  15 -minute  mark,  the  forecast 
errors  over  this  period  are  summarized  and  the  statistical  summary  files 
(files  15,    16,   and  17)  are  updated.     (See  Reference  3  for  discussion  of  these 
files.  )    If  a  link  has  failed  during  the  last  15  minutes,   its  error  statistics 
are  not  updated,   since  the  true  traffic  parameters  for  that  link  are  not 
available.    The  statistical  error  analysis  for  each  link  contain  the  following 
for  both  predicted  volumes  and  speeds: 
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•  Number  of  predictions 

•  Average  forecast  error 

•  Standard  deviation  of  the  forecast. errors 

•  Average  relative  forecast  error 

•  Standard  deviation  of  relative  forecast  error 

•  Average  absolute  relative  forecast  error  (AMI) 

•  Number  of  predictions  where  AMI  >  .  05 

•  Number  of  predictions  where  AMI  >  .  10 

•  Number  of  predictions  where  AMI  >  .  20 

The  present  form  of  the  traffic  volume  predictor  can  be  summarized  as 
follows: 

Let         r(t)       =        f(t)  -  m(t) 
'where    f(t)       =       observed  traffic  volume  at  time  t. 

m(t)     =        Fourier  series  value  for  traffic  volume  at  time  t. 

Let 

and 

The  predicted  value  at  time  t+1  is  then  computed  as: 
v(t+l)  =  m(t)  +  c(t)  +  d(t) 

The  speed  predictor  is  the  same  form  except  that: 
d(t)  =  .5h(t-l) 


c(t) 

= 

.9c(t-l)  +  .  lr(t-l), 

h(t) 

= 

r(t)  -  c(t), 

d(t) 

= 

.2h(t-l) 
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The  Fourier  series  values  m(t)  are  based  upon  historical  15-minute  traffic 
volume  and  average  speed  data  over  150  links  during  the  period  24  October 
1972  through  30  November  1972.  This  data  was  collected  on  the  UTCS  net- 
work during  the  time  period  0715  to  1900. 

After  the  statistical  files  have  been  updated,  forecasts  are  computed  for 
link  volumes  and  speeds.     The  model,   outlined  in  Reference  3  is  employed 
except  when: 

1)  The  actual  day  of  the  week  is  a  Saturday,   Sunday  or 
legal  holiday. 

2)  The  time  of  day  (for  the  prediction)  is  not  between 
0715  and  1900. 

3)  A  link  is  listed  in  the  LBAD  array,   indicating  that 
there  was  insufficient  data  available  to  estimate  the 
Fourier  coefficients  for  this  link. 

In  these  three  cases  the  actual  volumes  and  speeds  observed  over  the  last 
15  minutes  are  used  as  the  forecasts. 

The  Fourier  series  estimates   [m(t)]    of  the  traffic  parameters  are  obtained 
from  file  11.     This  file  contains  the  Fourier  series  values  for  each  link  for 
the  time  periods.     Corrections  are  then  applied  to  this  estimate,   except 
when  a  link  has  failed  during  the  last  fifteen  minutes.     In  this  case  the 
estimate  m(t)  is  used  as  the  forecast. 

If  the  predictor  was  executed  at  the  last  15-minute  mark  the  forecast  is 
computed  as: 
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v(t+l)        =     m  (t)  +  c  (t)  +  d  (t),  for  volumes 

s(t+l)        =     m  (t)  +  c  (t)  +  d  (t),   for  speeds 
s  s  s 


where 


c   (t)  =     .9c   (t-1)  +  .  lr  (t-1) 

V  V  V 

ca(t)  =     .9cQ(t-l)+  .lr  (t-1) 

fa  o  o 

r   (t-1)  =     f   (t-1)  -  m   (t-1) 
v  r  v 

r  (t-1)  =     f  (t-1)  -  m   (t-1) 
s  s  s 

h  (t-1)  =     r  (t-1)  -  c   (t-1) 

V  V  V 

h  (t-1)  =     r  (t-1)  -  c  (t-1) 
s  s  s 

d  (t)  =     .2  h  (t-1) 


v 


V 


d  (t)  =     .5h(t-l) 

f  (t)  =     observed  traffic  volume  at  time,  t 

v 

f  (t)  =     observed  traffic  speed  at  time,  t 

s 

m  (t)         =     Fourier  series  value  for  traffic  volume  at  time,  t 
m  (t)         =     Fourier  series  value  for  traffic  speed  at  time,  t. 

fa 

If  the  predictor  was  not  executed  at  the  last  15-minute  mark,  the  forecast 
is  computed  as: 

v(t+l)        =     mv(t)  +  cy(t) 
s(t+l)        =     rn(t)  +  c  (t) 

S  >3 

where  c   (t)  and  c   (t)  are  the  corresponding  corrections  that  were  applied  the 

v  fa 

last  time  the  predictor  was  called. 
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Routine  PREDICT  interfaces  with  the  real-time  second  generation  via  com- 
mon block  BLK3.     The  common  area  contains  the  arrays  PV(1)  and  PS(I), 
and  the  predicted  volumes  and  speeds  are  stored  there,    respectively. 


2.6.93.5    Listing 

See  Volume  III,    2GCF  Routine  Listings, 
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C  ENTER  J 

IPRD 

NN=  0 


I 


INPUT  CP 
ARRAY 


MIN=  60*HR  +  MIN 
TIME=  (MIN/5)-  84 


NN  =  0 

IT=  ITSAVE  +  1 
(UNIT  =  12 
IP=  NSTAT 


Routine  PREDICT 
100 


© 

Routine  PREDICT  (Continued) 
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KK  = 


0? 

s^YES 

KK  -  ± 

NO 

LL=  0? 


NO 


YES 


LL=  1 


RR(1)  =  JB(2)-PV(I) 
RR(2)  =  JB(5)-PS(I) 


© 


Routine  PREDICT  (Continued) 
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©■ 


STAT  (1, 

J, 

IP)  = 

=.STAT  (1 

J, 

IP) 

+  1 

STAT  (2, 

J, 

IP)  = 

STAT  (2, 

3, 

IP) 

+  RR(J) 

STAT  (3, 

J, 

IP)= 

STAT  (3, 

J, 

IP) 

+  RR(J)**2 

J  =  J  +  l 

RR(1)=  RR(1)/JB(2) 
RR(2)=  RR(2)/JB(5) 
J  =  1 


1 


STAT(4,J,  IP)=  STAT  (4,  J,  IP) 

+  RR(J) 

STAT  (5  J,  |P)=  STAT  (5,  J,  IP) 

+  RR(J)2 

RR(J)=  ABS(RRU)) 


I 


STAT  (6,  J,  IP)=  STAT  (6,  J,  IP) 
+  RR(J) 


Routine  PREDICT  (Continued) 
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©- 


YES 


STAT  (7,  J,  IP)  = 
STAT  (7,  J,  IP)  +  1 


YES 


STAT  (8,  J,  IP)  = 
STAT  (8,  J,  IP)  +  1 


YES 


STAT  (9,  J,  IP)  = 
STAT(9,  J,  IP)  +  1 


i 


1=1  +  1 


H© 


OUTPUT 

STAT 

ARRAY 


6 


Routine  PREDICT  (Continued) 
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CLEAR  PVAND 
PS  ARRAYS 
HOLDF=HLIDF 


YES 


YES 


IPRD=  1 
NN=  0 
IP=  NDATA 
LB=  1 


© 


CALL  LDM0E5 
TO  OBTAIN 
JB(2)&  JB(5) 


I 


PV(I)  =  JB(2) 
PS(I)  =  JB(5) 
1=1  +  1 


Routine  PREDICT  (Continued) 


105 


INPUT  HISTORICAL 
PREDICTION  CO- 
\EFFICIENTS  DATA 

T 


NN=  NN  +  13 
IP=  1 


L=0 

JJ=  DATA  (IT,  IP) 


Routine  PREDICT  (Continued) 
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€> 


LB=  LB  +  l 


-** 


CALL  LDM0E5 
TO  OBTAIN 
JB(2)&  JB(5) 


VP  =  JB(2) 
SP=  JB(5) 


YES 


YES 


YES 


CALL  LDM0E5 
TO  OBTAIN 
JB(2)&  JB(5) 


JJ=  DATAOT-1,  IP) 

SSS  =  JJ(2) 

RRR  =  JJ(1) 

R  =  JB(2)-RRR 

T=R-  CP(I) 

CP(I)  =  .9*CP(I)  +  .1*S 


I 


S  =  JB(5)-  SSS 

U=  S-  CS(I) 

CS(I)=  .9*CS(I)  =  .1*S 


© 


\Fv 


-»»/Ho  j 


"•■\H°) 


Routine  PREDICT  (Continued) 


107 


©- 
©- 


VP  =  VP  +  CP(I) 
SP  =  SP  +  CS(I) 


YES 


*© 


ITSAVE  \  NO 
=  IT-1 

? 


*® 


VP=  VP  +  .2*1 
SP=  SP  +  .5*U 


I 


PV(I)  =  VP 
PS(I)=  SP 


Routine  PREDICT  (Continued) 
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1  =  1 


CALL  LDMOES 
TO  OBTAIN 
JB(2)&  JB(5) 


PV(I)  =  JB(2) 
PS(I)  =  JB(5) 


ITSAVE=  IT 
NN=  0 


OUTPUT  CP 
ARRAY 


C  return") 


Routine  PREDICT  (Concluded) 
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2.6.94    SUBROUTINE  LDMOE5 

2.  6.94. 1    Purpose 

Subroutine  LDMOE5  retrieves  the  measures  of  effectiveness  data  computed 
by  RTN5  for  a  given  link.  Volume,  speed,  and  occupancy  are  returned  for 
the  specified  link. 

2.6.94.2  Relationship 

This  subroutine  is  called  by: 

CALL  LDMOE5  (LNK,  IVOL,  ISP,  IOCC) 

where: 

LNK  =  Input  argument  (first  generation  link  number) 
IVOL  =  Average  volume 
ISP  =  Average  speed 
IOCC  =  Average  occupancy 

This  subroutine  is  called  by  either  of  the  Level  5  Routines  RTSND,  PREDICT, 
or  LSTSQS.     The  subroutine  returns  to  the  calling  routine. 

2.6.94.3  Assumptions 
None. 
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2.  6.  94.4    Approach 

The  subroutine  searches  the  Level  5  buffer  (TBUFF,    loaded  by  RTN5)  to 
find  the  link  number  corresponding  to  the  value  of  LNK.     The  average 
volume,    speed,    and  occupancy  for  this  link  are  retrieved  and  returned  in 
the  output  arguments.     If  the  link  number  entered  is  255,    or  if  the  link 
fail  flag  (FAILP)  is  set,    or  if  the  link  is  not  found  in  the  buffer,    default 
values  for  average  volume,    speed,    and  occupancy  are  returned  (default 
volume  =  260,    default  speed  =  35,    default  occupancy  =  4%). 


2.6.94.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 


Ill 


1ST  =  IST  +  10 


YES 


IST=9 

NBK  =  (TBSIZE-4)/5 

1=0 

w 

1  =  1  +  1 

IV0L=  260 
ISP=  35 
I0CC  =  4 


C  RETURN^ 


Subroutine  LDMOE5 
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TBUFF  (IST  +  2) 
=  JLIM  (4) 


TBUFF  (JST  +  2) 
=  JLiM  (3) 


TBUFF  (IST  +  3) 
=  JLIM  (1) 


TBUFF  (IST  +  3) 
=  JLIM  (2) 


Subroutine  LDMOE5  (Concluded) 
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2.  6.  95.  ROUTINE  RTSND 

2. 6.95. 1  Purpose 

The  objective  of  Routine  RTSND  is  to  provide  intersection/ link  configurations 
to  be  optimized  by  the  network  optimization  Routine  LSTSQS.     In  addition, 
splits  at  major  intersections  are  computed  based  upon  current  and  predicted 
volumes.     Also,   cycle  lengths  are  established  for  each  subnetwork  to  be 
used  by  the  network  optimization. routine  to  compute  the  optimal  signal 
timing  patterns.     Checks  are  also  made  in  the  routine  to  determine  if 
traffic  conditions  have  changed  significantly  during  the  last  fifteen  minute 
period  to  justify  network  re -optimization. 

2.6.95.2  Relationship 

RTSND  is  called  from  Routine  MACRO  by: 

CALL  RTSND  (LEV5S,    NSUB) 

where: 

LEV5S  is  array  that  specifies  which  sections  are  in  traffic  responsive 
NSUB  a  return  parameter  that  specifies  the  number  of  core  groups 

RTSND  calls  subroutines  IRTN,    LDMOE5,    LPHASV,    SMSCYL  and  PSUB2. 

2.6.95.3  Assumption 
None. 

2.6.95.4  Approach 

The  methodology  of  the  subnetwork  determination  routine  is  based  upon  using 
a  predicted  volume  to  form  subnetworks.     The  procedure  used  involves  the 
following  steps: 
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a)      Webster's  equation  is  used  to  arrive  at  an  optimal  cycle 
length  for  each  major  intersection  based  upon  predicted 
volumes  on  each  major  intersection  approach.     The 
specific  formula  for  each  intersection  is: 


1.  5  K+5 


C    = 


1  -  Y 

where: 

C    =    Cycle  length 

K    =    Lost  time  for  all  phases  due  to  starting  delay 

Y    =    Sum  of  ratios  of  predicted  critical  lane  volume 
to  saturation  volume  on  each  approach. 

b)  Limit  the  cycle  lengths  computed  in  step  a)  to  a  minimum 
equalling  the  greatest  minimum  cycle  length  of  all  con- 
trollers in  the  core  group  and  a  maximum  of  120  seconds. 

c)  In  each  of  six  core  groups  of  major  intersections,   determine 
maximum  cycle  length  found  in  step  a). 

d)  Combine  neighboring  core  groups  having  cycle  lengths  which 
are  nearly  equal  (currently  set  as  10  seconds). 

The  UTCS  network  core  groupings  currently  used  are  listed  in  Table  1. 
These  groupings  were  selected  because  of  the  expected  differences  in  flow 
conditions  on  each  of  the  major  arterials  in  each  control  section. 
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Table  1.     UTCS  Network  Intersection  Core  Groups 


Section 

Core 

Controller 

Number 

Group 

Number 

Name 

1 

3 

15 
17 

M  Street 

2 

1 

120 

9 

10 

11 

Upper  Wisconsin 

2 

2 

121 

123 

12 

Lower  Wisconsin 

3 

4 

84 
86 
82 
97 
72 
59 
61 
88 
99 

K  and  L  Streets 

3 

5 

70 
65 
91 
66 
68 
94 
55 

Penn.  Avenue,   H  and  I  Stree 

4 

6 

53 

Constitution  Avenue  and  17th 

31 

Street 

50 

33 

40 

41 

42 
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Routine  RTSND  is  called  by  the  MACRO  loop  executive  routine.     Upon  entry 
the  flags  LEV5S  are  tested  to  determine  which  control  sections  have  been 
selected  by  the  operator  for  the  traffic  responsive  mode  of  control.     Based 
upon  which  sections  have  been  selected,   "scratch"  group  numbers  are  placed 
in  the  array  IGRUP  to  indicate  which  neighboring  groups  should  be  con- 
sidered for  combining  into  subnetworks. 

Subroutine  IRTN  is  then  called  for  each  of  the  "scratch"  groups  to  obtain  the 
intersection  numbers  in  each  group,  the  approach  link  numbers  to  each 
intersection,  and  the  minimum  cycle  length  of  all  the  intersections  in  each 
group. 

A  "Do-Loop"  is  then  entered  to  process  each  major  intersection  in  each 
core  group. 

The  predicted  volumes  on  each  approach  link  to  each  major  intersection  are 
obtained  from  the  array  PV.  Intersection  splits  are  then  computed  for  each 
major  intersection  via  the  following  equation: 

VPHASA 
SPLIT    =    


VPHASA  +  VPHASB 

where: 

VPHASA  is  A  phase  critical  approach  volume 
VPHASB  is  B  phase  critical  approach  volume 

The  critical  approach  for  each  phase  is  defined  as  the  approach  which  has  the 
maximum  volume  out  of  all  links  associated  with  the  phase. 

A  C -phase  split  is  computed  in  the  same  manner  for  three  phase  inter- 
sections.    Subroutine  SPLPK  is  called  to  pack  the  splits  into  the  data  base. 
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The  cycle  length  then  computed  for  the  intersection  via  Webster's  equation 
described  earlier  and  the  maximum  cycle  length  found  in  each  "scratch" 
group  is  saved. 

The  controller" Do- Loop"  is  repeated  for  each  "scratch"  group.     The  maxi- 
mum cycle  lengths  of  neighboring  groups  are  then  compared.     If  they  are 
within  a  tolerance  (currently  10  seconds)  the  groups  are  combined  into  a 
subnetwork.     The  groups  contained  in  each  subnetwork  are  stored  in  the 
array  ISB  for  use  by  the  network  optimization  routine.     No  more  than  three 
core  groups  may  be  combined  into  one  subnetwork  due  to  the  dimension 
limitations  of  the  network  optimization  package.     These  dimensional  limi- 
tations are  due  to  processing  time  being  a  quadratic  function  of  the  size  of 
the  subnetwork. 

After  the  core  groups  have  been  arranged  into  subnetworks,  the  ISB  array 
is  tested  to  determine  which  core  groups  have  been  combined.     The  sub- 
network interfacing  flags,  INTF,   which  are  all  set  upon  initial  entry  to 
Routine  RTSND  are  reset  for  those  group  interfaces  which  are  not  to  be 
optimized  because  the  groups  have  combined  into  one  subnetwork. 

Tests  are  then  made  to  determine  if  the  number  of  subnetworks  to  be  optimi- 
zed is  different  than  the  number  that  were  optimized  on  the  last  15-minute 
pass  through  the  MACRO  loop.     If  so,  the  flag,  IREOPT,  is  set  to  indicate 
to  the  MACRO  loop  executive  that  re -optimization  should  occur. 

In  addition,  the  maximum  cycle  lengths  for  each  subnetwork  determined  by 
Webster's  approach  described  earlier  are  compared  to  the  cycle  lengths 
currently  being  used  by  each  subnetwork.     If  any  is  different,  the  flag, 
IREOPT,   is  set  to  indicate  to  the  MACRO  loop  executive  that  re -optimization 
should  occur. 
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For  the  controllers  whose  corresponding  split  hold  flag  (HLDSPL)  is  set, 
the  current  rather  than  the  optimitized  split  values  are  retained. 

Before  exiting  to  the  MACRO  loop  executive  the  cycle  length  to  be  searched 
by  the  network  optimization  routine  is  established  by  storing  the  value 
computed  by  the  Webster's  equation  into  the  CYCL  array  for  each  sub- 
network. 


2.6.99.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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(   ENTER  J 

^r 

1SCT(1-124)  =  SPLITd-124) 
INTF(l-5)=  0 
ISECT(l-8)=0 
CYCC1-8,  1-5)=  0 
NCYCLE(l-8)=0 
ISB(l-8,  1-8)=  0 


0. 


i 


IGRUPU-8),  1-8)  =0 

L=0 

M=0 


3 


M=  M  +  1 


YES 


*© 


Routine  RTSND 


120 


Routine  RTSND  (Continued) 


121 


© 


IGRUP(K, 

L)  = 

1 

IGRUPdO 

1, 

L)  = 

2 

IGRUP(IO 

2, 

L)  = 

3 

ISECT(l) 

=  1 

ISECK2) 

=  1 

ISECK3) 

=  1 

I 


INTF(4)  =  1 
INTF(3)  =  1 
M=  M  +  l 


LEV5S(3)\I\I0 
=  3? 


*© 


Routine  RTSND  (Continued) 
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IGRUP(K,  L)=  1 
IGRUP(K  +  1,  L)=  2 
ISECT(1)=  1 
ISECT(2)=  1 
INTF(4)=  1 


IGRUP(K,  L)=4 
IGRUP(K  +1,  L)=  5 
INTF(1)=  1 
ISECT(4)=  1 
ISECT(5)=  1 


IGRUPCK,  L)=  6 
ISECT(6)=  1 


Routine  RTSND  (Continued) 
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J=  LEV5SCM) 
NN=  M  +  1 
N=  NN 


J  = 

AND 


=  1       \YES 

K=  2     > : *" 

INTF(3)=  1 

?     ^^ 

s/NO 

J=  1        > 

v  YES 

AND  K  =  3 

INTF(^;  =  l 

N.       7        S^ 

>^NO 
W— — 

AND 


-  3       \YES 

K  =  4    >>               ** 

INTF(5)=  1 

>j/NO 

Fo 


Routine  RTSND  (Continued) 
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M  =  M  +  1 


Routine  RTSND  (Continued) 


125 


©■ 


CYMAX(l-8)=  0 
ITKMIN(l-8)=  0 
ICOUNT=  0 
M=  1 


5 


N=  1 


YES 


ICOUNT  =  ICOUNT  +  1 
IGP=  IGRUP(N,  M) 


I 


CALL  IRTN 


I 


ITKMIN(ICOUNT)  =  ITKI 
INT  =  1 


3 


ISN  =  INTRS(INT) 


YES 


*© 


*© 


Routine  RTSND  (Continued) 
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CALL  LPHASV 


CALL  SMSCYL 


I 


SUMT  =  £L^pjMI2 
ILN=  1 


Routine  RTSND  (Continued) 
127 


Q 


VSUM=  0.0 
TL0SS=  0.0 


YES 


YES 


VA=  VPHASA 


© 


VPHASA  =  PV(LNKl) 


VPHASA  =  PV(LNK2) 


VPHASA  =  PV(LNK3) 


VPHASB=  PV(LNK4) 


VPHASC=  PVCLNK5) 


Routine  RTSND  (Continued) 
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K3PH  = 


:    2? 

vYES 

^~ 

VA  = 

^ 

AMAXKVA,  VPHASC) 

NO 

VB= VPHASB 


YES 


VB  = 

AMAXKVB,  VPHASC) 


VSUM=  VA  +  VB 


=  1? 

v  YES 

1  LUbb  -   1  Lu^j  t  3.u 

o° 

Routine  RTSND  (Continued) 
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Y=  VSUM*CONS 


Y> 


.95 

w    YES 

^- 

Y=  O.vd 

NO 

CYCLE     d»5(TL0SS+5) 


(l.-Y) 
CYCLE  =  0.85*CYCLE 
ICYSAVE=  CYCLE 


YES 


YES 


CYCLE  *         \    YES 
CYMAX(ICOUNT) 


CYMAX(ICOUNT)=  CYCLE 


X 


INT  =  INT  +1 


CYCLE  =  ITKMIN 
(ICOUNTJ  +  5 


CYCLE = 120 


Routine  RTSND  (Continued) 
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N=  N  +  l 


M=  M  +  1 


*© 


Routine  RTSND  (Continued) 
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NSUB=  1 
M=  0 
K  =  0 
IM  =  0 
L=  1 


K=  K  +  l 
IM  =  IM  +  1 

M=  M  +  1 


YES 


ISBdM,  NSUB) 
=  IGRUP(M,  L) 


CVCL(NSUB)  = 
CYMAX(K  +  1) 


YES 


YES 


YES 


YES 


CYMAX(K  +  1) 
=  CYMAX(K) 


-MM 


Routine  RTSND  (Continued) 
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CYCL(NSUB)=  CYMAX(K) 
ISBdM,  NSUB) 
=  IGRUPCM,  L) 
NSUB=  NSUB  +  1 
IM=  0 


Routine  RTSND  (Continued) 
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B(2,  1)    \YES 

INTF(4)=  0 

=  2?        / 

|N0 

m 

IS 


B(3,  1)    \  YES                 _ 

INT-F(3)=  0 

=  3?        y? 

l« 

SB(2;  2)x 

s^YES 

=  3?       . 

INTF(3;  -  u 

|NO 
U« 

ISB(2,  2)\YES 


5? 


NO 


INTF(1)=  0 


ISB(2,  3)   X.YES 

INTFCL)  -  0 

=  bV         / 

|NO 

La 

ISB  (2,  4)  ^ 

V   YES 

INTF(1)=  0 

=  5? 

[NO 
La 

INTFCL)  ^0 

v    YES 

INTF(2)=  0 

|   NO 
m 

Routine  RTSND  (Continued) 
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1  =  1 


I 


JCYCL=  CYCL(I) 
MCYCL =JCYCL 


MCYCL    \YES 

MCYCL  =  MCYCL  -  100 

>100.''    />                        •" 

La 

KrYr.   _  MCYCL 
K.LYLL-  !0*10 

KCYCL=MCYCL-KCYCL 
ICYCLE=  0 


KCYCL 

\YES 

ICYCLE  =  5 

>  2? 

^ 

plo 

La 

KCYCL 

VYES 

ICYCLE  -  10 

>7? 

Le — 

CYCCl,  1)  =  JCYCL  + 
ICYCLE  +  KCVCL 
NCVCLE(I)=  1 
1=1  +  1 


Routine  RTSND  (Continued) 
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N 


SUBO    v 
NSUB 

\.  YES 

IRE0PT=1 

?    s^ 

jfio 

©■ 


rQ 


NSUBO  =  NSUB 
M=  1 


3 


CYCLE  =  CYC(M,  1) 
C0NS1  =  CONS*CYCLE 
N=  1 


3 


IGP  =  'ISN(N,  M) 


YES 


CALL  IRTN 

(IGP,  INTRS,  ITKM) 


I 


SCT0T=0.0 
SCCNT  =  0.0 
INT  =  1 


Routine  RTSND  (Continued) 


*© 
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GAMIN 
GAMIN 
GBMIN 
GBMIN 
GCMIN 
GCMIN 


IGAMIN 

GAMIN/2 

IGBMIN 

GBMIN/2 

IGCMIN 

GCMIN/2 


I 


CALL  SMSCYL 

I 


GET  THE  LOST  TIME  AND  THE 

CYCLE  LENGTH  FOR  THE  INTERSECTION 


SUMT=  ISUMT 
SUMT=  SUMT/2 
GABC=CYCLE-SUMT 
ILN=  1 


Routine  RTSND  (Continued) 
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Routine  RTSND  (Continued) 


138 


YES 


YES 


YES 


YES 


GDA=  VPHASA*C0NS1 
GDB=  VPHASB*C0NS1 
GDC  =  VPHASC*C0NS1 


VPHASA=  PV(LNKl) 


VPHASA=  PV(LNK2) 


VPHASB=  PV(LNK3) 


VPHASB=  PV(LNK4) 


VPHASC=  PV(LNK5) 


Routine  RTSND  (Continued) 
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rtIN 

VYES 

* 

GA=  GAMIN 

0 

GABC 

<GBM 


:-ga^ 

MM?/ 

vYES 

** 

GA=  GABC-GBMIN 

NO 

GA 


SPLITA  =  GABC 

ISPLT=  SPLITA*1000  +  5 

ISPLT2=  0 


Routine  RTSND  (Continued) 
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£ 


bA     GDA  +  GDB=GDC 
GC-GDC 


GDA=  GDB=  GDC 


*GABC 
-*GABC 


GA<  GAMIN^ 

YES 

^- 

GA-  Gamin 

|N0 
m 

GC<GCMIN 


YES 


GC  =  GCMIN 


GABC-(GA  +  GC)\YES 
*  GBMIN? 


NO 


GA=   (GDAGD+AGDC)   3ABC-GBMIN 


tA<  GAMIN' 


YES 


GA=  GAMIN 


■GABC-(GA-GBMIN> 

s  YES 

GA=  GABC  - 

.       <GCMIN?       > 

^ 

(GCMIN  +  GBMIN) 

JNO 

l^ — . 

GC  =  GARC- 
(GA  +  GdMIN) 


Routine  RTSND  (Continued) 
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® 


SPL!TA=GA/GABC 
ISPLT=  SPLITA*1000  +  5 
SPLIT2=GC/GABC 
ISPLT2=  SPLIT2*1000  +  5 


© 


© 


GDAC  =  AMAXKGDA,  GDC) 


GB  = 


-GOB. 


GDAC  +  GDB 


*GABC 


GB<GBMIN* 

SYES 

GB=GBMIN 

\^NO 

GABC-GB^ 
<GAMIN?    . 

v  YES 

* 

GB=  GABC-GAMIN 

iiO 

GAC  =  GABC-GB 
SPLITA=  GAC/GABC 
ISPLT=  SPLITA*1000  +  5 
ISPLT2=0 


Routine  RTSND  (Continued) 
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GDBC  =  AMAX1  (GDB,  GDC) 


GA 


&DA- 


GDA  +  GDBC 


*GABC 


YES 


GBCMIN  = 
AMAXKGBMIN,  GCMIN) 


GABC 
GBCf 


-GA^ 

v   YES 

GA  =  GABC-GBCMIN 

NO 

SPLITA=  GA/GABC 
ISPLT=  SPLITA*I000  +  5 
ISPLT2=  0 


Routine  RTSND  (Continued) 
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JSN  =  2*(ISN-1)  +  1 


I 


JSHF=HLDSPL(ISN) 


YES 


ISCT(JSN)  =  ISPLT 
ISCT(JSN1)=  ISPLT2 
SPLIT(JSN)=  ISCAC 


I 


SCAC=  FLOAT(ISCAC)*XD 
SCAN=  FLOAT(ISPLT)*XD 


SCTOT=  SCTOT  + 


ABS(SCAC-SCAN) 


"SCAU" 


SCCNT=  SCCNT  +  1.0 


YES 


•© 


*© 


Routine  RTSND  (Continued) 
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<s> 


ISCC(JSN  +  1)=  SPLIT 
SCCC  =  FLOAT(ISCCC)*XD 
SCCN=  FL0AT(ISPLT2)*XD 

ABS(SCC-SCCN) 
5CCC 


SCTOT=  SCTOT  + 


SCCNT  =  SCCNT  +  1 


3 


INT=  INT  +  1 


N=  N  +  1 


*® 


*© 


^ 


Routine  RTSND  (Continued) 
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CYCLOdTN 
*  CYCL(I)    > 

YES 

IREOPT=  1 

^Sj^NO 

INTFO(I) ^V  YES 

IREOPT  =  1 

/  INTKI)     J>                        w 

L^° 

Routine  RTSND  (Continued) 
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KSPHLD(I)=  HLDSPL 


SPHLD 


)=  0  7s 

YES 

SPLIT(I)  =  ISCT(I) 

NO 

Routine  RTSND  (Concluded) 
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2.6.96    SUBROUTINE  IRTN 

2.6.96.  1    Purpose 

Subroutine  IRTN  accesses  the  array  INCG  for  a  given  group  number  and 
returns  the  instrumented  controller  numbers  in  the  group  in  array  INT. 
In  addition  all  controllers  in  the  group  are  examined  to  determine  the 
largest  minimum  cycle  length  of  all  the  controllers  in  the  group. 

2.6.96.2  Relationship 

Subroutine  IRTN  is  called  as  follows: 
CALL  IRTN  (IGR,  INT,   ITKM). 

where: 

IGR  =  Input  argument  (group  number) 

INT  =  Array  of  instrumented  controller  numbers 

ITKM  =  Largest  minimum  cycle  length  in  group. 

Subroutine  IRTN  is  called  by  Routine  RTSND  and  returns  to  that  routine. 

2.6.96.3  Assumption 


None, 
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2.6.96.4    Approach 

The  instrumented  controller  starting  index  for  the  group  is  determined  by- 
accessing  the  ICGSI  array  for  the  group.     The  INCG  array  is  then  trans- 
ferred into  the  output  array  INT  for  all  instrumented  controllers  in  the 
group. 

The  controller  starting  index  for  the  g^oup  is  then  determined  by  accessing 
the  GCNTSI  array.     The  minimum  cycle  length,    TKM,    for  each  controller 
in  the  group  is  accessed  through  a  call  to  Subroutine  GETBY  and  the  largest 
TKM  is  stored  in  the  output  argument,    ITKM. 


2.6.96.5    Listing 


See  Volume  III,    2GCF  Routine  Listings, 
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r  enter) 

_£_ 

CLEAR  INT 
ARRAY 


I 


ITKM  =  0 
1  =  0 


1=1  +  1 


I 


N=  IST  +  1 
INT(I)=  INCG(N) 


NO  v^NICIGP(IGR) 
? 


i 


1=1  +  1 
N=  IST  +  1 

ITM=  (TKM(ICON)  +  1) 
2 


Subroutine  IRTN 
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2.6.97     SUBROUTINE  LPHASV 

2. 6.97. 1  Purpose 

Subroutine  LPHASV  loads  first  generation  link  numbers  by  phase  for  a  given 
controller.     In  addition,  the  controller  type  (K3PH)  is  returned  along  with 
time  available  for  A  and  B  phases  (GAB),  time  available  for  B  and  C  phases 
(GBC),   and  time  available  for  A,   B,   and  C  phases  (GABC).     Also  returned 
are  the  minimum  phase  durations  for  A,   B,   and  C  phases,   respectively. 
This  data  is  used  by  Routine  RTSND  to  compute  splits  at  instrumented  con- 
trollers in  each  core  group. 

2.6.97.2  Relationship 

The  calling  sequence  for  this  subroutine  is  as  follows: 

CALL  LPHASV  (ISN,   LNK,   K3PH,    GAB,    GBC,    GABC,    GAM, 
GBM,   GCM) 

where: 

ISN  =  Input  argument  (first  generation  controller  number) 
LNK  =  Array  of  link  numbers  associated  with  the  controller 


where: 


LINK(l)  =  Phase  A 1  link  number 
LINK(2)  =  Phase  A2  link  number 
LINK(3)  =  Phase  Bl  link  number 
LINK(4)  =  Phase  B2  link  number 
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LINK(5)  =  Phase  CI  link  nutnber 
Link(6)  =  Phase  C2  link  number 
K3PH  =  Controller  type 

=  0  two  phase,   no  overlap 

=  1  three  phase,   no  overlap 

=  2  three  phase,   A  and  C  overlap 

=  3  three  phase,   B  and  C  overlap 
GAB  =  Time  available  for  A  and  B  phases  (sec,    scaled  Bl) 
GBC  =  Time  available  for  B  and  C  phases  (sec,    scaled  Bl) 
GABC  =  Time  available  for  A,    B,    and  C  phases  (sec,    scaled  Bl) 
GAM  =  Minimum  A  phase  time  (sec,    scaled  Bl) 
GBM  =  Minimum  B  phase  time  (sec,    scaled  Bl) 
GCM  =  Minimum  C  phase  time  (sec,    scaled  Bl) 

This  subroutine  is  called  by  Routine  RTSND  and  returns  to  that  routine. 

2.6.97.3  Assumption 
None. 

2.6.97.4  Approach 

The  subroutine  accesses  arrays  LSI,   CONTLN,   and  LDSET  to  determine  the 
link  numbers  and  their  phase  association  for  controller  ISN.     Arrays 
GABBCABC  and  GMINABC4  are  accessed  to  determine  the  allowable  phase 
times  and  minimum  phase  times. 

2.6.9  7.5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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(enter) 


+ 

CLEAR  LNK  ARRAY 

IDX  = 

4*(INT-1)  +  1 

GAB  = 

GGGT(IDX) 

GBC  = 

GGGTODX  +  1) 

GABC 

=  GGGKIDX+  2) 

GAM  = 

=  'GGGG(IDX) 

GBM  = 

=  GGGGCIDX  +  1) 

GCM: 

=  GGGGCIDX  +  2) 

K3PH    ^  2^ 

YES 

fc 

GCM  = 

W 

SSSSCIDX  +  2) 

\^0 

1ST  =  6*(INT-1) 

ISTT  =  IST+  NQLCINT) 


-►(return  ) 


Subroutine  LPHASV 
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2.6.98    SUBROUTINE  SMSCYL 

2.6.98.1  Purpose 

Subroutine  SMSCYL  fetches  the  lost  times  following  the  A,    B,    C  phases 
and  the  total  lost  time  per  cycle. 

2.6.98.2  Relationship 

Subroutine  SMSCYL  is  called  from  Routine  RTSND  by: 

CALL  SMSCYL  (INT,    SUMA,    SUMB,    SUMC,    SUMT,    CYCLE) 

where: 

INT  denotes  the  intersection 

SUMA  returns  the  lost  time  following  A  phase 

SUMB  returns  the  lost  time  following  B  phase 

SUMC  returns  the  lost  time  following  C  phase 

SUMT  returns  the  total  lost  time. 

CYCLE  returns  the  cycle  length. 

It  returns  to  the  calling  routine. 

2.6.98.3  Assumption 
None. 

2".  6.98.4    Approach 

Subroutine  SMSCYL  makes  repeated  calls  to  subroutine  GETBY  to  retrieve 
the  cycle  length,    the  lost  times  following  the  A,    B,    C  phases  and  the  total 
lost  time  per  cycle.     SMSCYL  then  returns  to  the  calling  routine. 

2.6.98.5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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(enter) 


CYCLE  =  KCL(ENTCYCdNT)) 
IP0SN  =  4*(INT-1) 


SUMA=  SSSSCIPOSN  +  1) 
SUMB=  SSSS(IPOSN  +  2) 
SUMC=  SSSSCIPOSN  +  3) 
SUMT=  SMTABCdPOSN  +  4) 


(^RETURN  J 


Subroutine  SMSCYL 
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2.6.99    ROUTINE  LSTSQS 

2.6.99.1  Purpose 

Routine  LSTSQS  determines  the  signal  settings  (cycle  length  and  offsets)  for 
intersections  such  that  vehicle  delay  is  minimized  in  each  subnetwork  deter- 
mined by  Routine  RTSND. 

2.6.99.2  Relationship 
LSTSQS  is  called  by  MACRO. 

LSTSQS  calls  LINKD,    LDMOE5,    LLGTMF,    CORRECT,    and  LOADAT. 

2.6.99.3  Assumption 
None. 

2.6.99.4  Approach 

The  methodology  used  by  the  optimization  algorithm  coded  in  Routine  LSTSQS 
is  based  upon  the  SIGOP  Traffic  Signal  Optimization  Program  described  in 
Reference  5  and  modified  by  Reference  6.     The  algorithm  assumes  that 
vehicle  delay  on  a  link  in  the  traffic  signal  network  is  a  quadratic  function 
with  respect  to  offset  difference  between  the  upstream  and  downstream 
nodes,  (controllers)  on  the  link.     The  assumption  is  an  approximation  to  the 
true  relationship  between  delay  and  offset  difference.     However,  the  approxi- 
mation is  very  accurate  in  the  neighborhood  of  the  minimum.     The  over- 
whelming advantage  of  using  this  approximation  is  that  the  system  of  linear 
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equations  which  result  when  applying  the  necessary  conditions  for  a  minimum 
from  the  calculus  can  be  solved  very  rapidly  in  relationship  to  other  approaches 
which  assume  piecewise  linear  delay  functions,   or  conduct  steepest  ascent 
searches  over  an  unknown  hypersurface  to  seek  a  minimum.     Because  of  the 
great  advantage  in  computational  speed,  the  quadratic  delay  approximation 
lends  itself  nicely  to  real-time,   on-line  network  signal  timing  optimization. 
In  fact,   it  is  the  only  technique  to  date  which  meets  the  speed  requirements 
necessary  for  real-time  optimization  of  a  signal  network. 

The  explicit  objective  function  assumed  in  arriving  at  the  LSTSQS  algorithm 
is: 

D=2S  A..  (R.    +  M      +  d    -  d.)2  (1) 

1        1  -I  o  o  J  J 

where: 

D  =  Assumed  network,  delay 

A. .  =  A  link  weighting  factor  for  the  link  connecting  node  i  to 

J 

j  composed  of  link  volumes  and  a  link  importance  factor 

R..  =  Link  idealized  offset  for  the  link  connecting  node  i  to  i 

(computed  as  the  link  length  divided  by  the  average  speed 
of  the  vehicles  traveling  the  link) 

d.  =  Offset  of  node  (intersection)  j 

d.  =  Offset  of  node  i 
l 

R. .,   d-,   d.  are  normalized  by  cycle  length  in  the  above 
expression. 

M..  are  integers  chosen  because  of  periodicity  of  the  offset 
such  that 
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-  1/2  s(R..  +  M..  +d.  -  d.)  <  1/2  (2) 

13  13        3         1 

The  LSTSQS  Routine  solves  for  the  d.  from  the  equations  which  arise  from 
applying  the  necessary  conditions  to  the  above  objective  function.     These 
optimal  offsets  are  determined  for  a  range  of  candidate  cycle  lengths  and 
the  objective  function  is  evaluated  for  each  set  of  offsets.     The  cycle  length 
and  associated  offsets  which  provide  for  the  least  delay  are  chosen  for 
implementation  and  packed  into  the  entry  offset  array. 

The  solutions  obtained  from  the  system  of  linear  equations  provide  for 
excellent  progressions  in  signalized  networks  wherein  a  cycle  length  can 
be  implemented  which  allows  the  delay  function  to  be  driven  near  zero. 
Through  observation  of  the  objective  function  of  Eq  (1),   it  can  be  shown 
that  these  conditions  can  be  achieved  by  employing  a  cycle  length  which  is 
twice  the  average  travel  times  on  all  links  in  the  network. 

To  choose  such  a  cycle  length  is  not  always  possible,   however.     In  networks 
where  block  lengths  are  very  short  (less  than  400  feet),  travel  times  are 
such  that  a  cycle  length  which  is  twice  the  average  travel  time  may  be  less 
than  40  seconds.     It  is  not  feasible  to  implement  such  a  small  cycle  length 
in  heavily  congested  networks. 

When  larger  cycle  lengths  are  imposed  in  the  objective  function  of  Eq  (1)  for 
networks  of  this  type,    the  assumption  that  a  quadratic  delay  relationship 
exists  tends  to  break  down.     This  breakdown  occurs  because  a  solution  is 
forced  which  is  not  near  the  region  of  convergence  as  described  earlier. 

As  a  result,  the  offset  errors  are  relatively  large.     That  is,  vehicles  re- 
leased from  one  node  (intersection)  arrive  at  the  downstream  node  well 
after  the  downstream  signal  changes  to  green.     This  offset  error  is  not 
accounted  for  in  establishing  the  offset  for  the  second  downstream  node. 
That  is,  the  objective  function  as  written  in  Eq  (1)  assumes  that  vehicles 
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are  released  at  the  upstream  nodes  at  the  start  of  green,   i.  e. ,  with  zero 
offset  error. 

This  shortcoming  is  easily  overcome  by  including  continuity  links  in  the 
network  as  shown  in  Figure  1  below. 


Figure  1.     Simple  Network  Demonstrating  Continuity  Links 

The  objective  function  for  the  simple  network  of  Figure  1  then 
becomes 

D  =  A21  (R21  +  M21  +  d2  -  dl)2  +  A12  (R12  +  M12  +  dj  -  d/ 

+  A32  <R32  +  M32  +  d3  "  d/  +  A32  <R32  +  M32  +  d3  "  d/ 
+  A31  (R31  +  M31  +  d3  -  dt)2  +  A13  (R13  +  M13  +  dt  -  d3)+ 

where  the  following  constraints  hold 


R13  +  M13=R12  +  M12+R23  +  M23 


(3) 


and 


R31+M31=R32  +  M32  +  R21  +  M21 
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When  the  constraints  of  Eq  (3)  are  substituted  into  Eq  (1),   and  the  necessary- 
conditions  for  optimization  are  applied  (i.e.,   taking  the  partial  derivatives 
with  respect  to  d.,   setting  to  zero,   and  solving  for  the  d.)  the  following 
equations  result. 

N-l 
d.    =      2        E..C.  (4) 

3  =  1 

dN  =  ° 


where 


N  is  the  number  of  nodes  (intersections)  in  the  subnetwork 

d.  are  the  optimal  offsets 

E-  •  is  the  inverse  of  the  matrix  B. . 

where 


Bjk  " "  <Akj  +  V  i '  K 
Bkk  ■  f  %  +  V 

C,    is  the  vector 
k 

-f<A3kT3k-AkJVk=1 N 

where 


Tij  ■  M«  +  Rij 
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Because  of  periodicity,   many  local  minima  exist  in  offset  space.     Each 
minima  depends  upon  the  choice  of  the  M...     Thus,   a  'gradient  type'  search 
through  neighboring  local  minima  is  performed  starting  from  a  randomly 
selected  initial  point  to  obtain  a  second  order  minimum.     The  starting 
points  are  not  preselected.     A  second  order  minimum  is  one  which  yields 
a  value  of  the  objective  function  which  is  smaller  than  any  of  the  neighboring 
minima.     Several  second  order  minima  may  exist.     Therefore,   a  series 
of  "Monte  Carlo"  games  is  played,   each  starting  from  a  randomly  selected 
starting  point.     The  values  of  the  objective  function  for  all  second  order 
minima  obtained  from  all  random  games  are  compared  to  determine  the 
global  minimum. 

The  technique  used  in  the  gradient  type  search  is  to  choose  the  subscripts 

of  the  following  expression  which  yield  the  minimum  value  for  the  expression: 


X  -    |  Y  |  (5) 

m,n       ■     m,n' 


where 


Xm,n  -T      Aln    [Aln  Vn  +  *]    +   Amn[Amn  »ma  +1] 

+  A  [A        v  +ll+A,     A  f&  -e         +  <f>        +  P       1 

mo    I     mo  T  mn  In     mo    I    mn        mn        mn      ymnj 


1    +  A        A         fp         +v        "  tf      1 1 
1  mo     mn  rmn       'mn        rnnlr 


+  A,     A  Tv        -  e         +B|+A        A         [fl         +  y        -  $ 

In     mn      'mn        mn       Hmn  mo     mn  rmn       'mn        mn 


and 


Y  =  A,     T      +  A        T         +  A        T 

m,  n         In     in         mn     mn         mo      mo 

where  subscripts  1  and  o  for  link  m,   n  are  as  indicated  below. 
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I       1      J W     m      j  w      n       V 


The  terms  in  the  expression  for  X  are  as  follows: 

^  m,  n 

Y  =  E,      -  E„  -  E       +  E   . 
T  mn         In  11         nn         nl 

B  =E         -E  -E+E 

w  mn         mn         mm  nn         nm 

Y  =E         -E  -E+E 
mn         mo         mm         oo         om 

6  =  E,      -  E„  -  E       +  E   . 

mn         lo         11         oo         ol 

e  =  E,       -  E„  -  E  +  E     . 

mn         lm         11         mm         ml 

</>  =  e        -E        -E       +E 

mn         no         nn         oo         on 

The  subscripts  selected  indicate  which  M. .  should  be  changed  by  unity  to 
achieve  a  minimum  smaller  than  the  current  minimum. 


A  second  order  minimum  is  one  at  which  no  choice  of  the  subscripts  m, 
n  in  Eq  (5)  yields  a  negative  value  for  the  expression. 

The  mathematical  basis  for  the  equations  are  given  in  Reference  6. 

The  LSTSQS  Routine  is.  called  by  Routine  MACRO  of  the  Priority  Level  5 
loop.     The  LSTSQS  Routine  performs  a  "Do  Loop"  on  all  subnetworks  to  be 
optimized.     Upon  entry  into  this  loop  Subroutine  LINKD  is  called  to  obtain 
the  following  data: 

a)  The  first  generation  controller  numbers  in  each  subnetwork. 

b)  The  first  generation  link  numbers  connecting  the  controllers. 
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c)  The  first  and  second  controller  upstream  numbers  associated 
with  each  link. 

d)  The  first  and  second  controller  numbers  associated  with  each 
link. 

e)  The  second  generation  link  numbers  associated  with  each  link. 

Predicted  speeds  and  volumes  are  then  obtained  for  each  of  the  major  links 
from  the  PS  and  PV  array  respectively.     Subroutine  LLGTHF  is  called  for 
each  link  to  obtain  link  lengths  and  link  importance  factors.     The  idealized 
offsets  (the  R..  of  Equation  (1)  are  computed  for  each  link  and  the  weighting 
factors  (the  A.,  of  Equation  (1)  are  then  computed  for  each  link.     The 

J 

idealized  offsets  and  weighting  factors  are  then  stored  into  the  FIDEAL 
array  and  AMN  array  for  each  link. 

A  subscript  matching  loop  is  then  entered  to  automatically  insert  the  con- 
tinuity links  into  the  subnetwork  to  be  optimized. 

The  E  matrix  and  vectors  X  and  Y  of  Eq  (4)  and  (5)  are  then  computed  for 

« 

the  subnetwork  under  consideration.    A  series  of  random  games  is  then 
played  using  the  expression  of  Eq  (5)  to  determine  the  M..  and  associated 
offsets  which  yield  the  minimum  value  for  the  objective  function  of  Eq  (1). 
The  starting  points  for  the  random  games  are  not  preselected. 

The  above  process  is  repeated  for  each  cycle  length  in  the  cycle  length 
array,   CYC,   established  by  Routine  RTSND  for  each  subnetwork.     The 
cycle  length  and  offsets  which  yield  the  minimum  value  of  the  objective 
function  are  saved  in  the  arrays  CL  and  AE  respectively.     LSTSQS  then 
computes  the  variable  interval  durations  for  each  controller  which  main- 
tains the  optimal  cycle  length  and  fixed  intervals  at  each  controller. 
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After  all  subnetworks  have  been  optimized,  Routine  LSTSQS  calls  LOADAT 
to  insure  that  the  relationships  formed  when  an  operator  modified  the  splits 
and/or  offset  is  maintained.     Routine  LSTSQS  then  returns  to  the  MACRO 
loop  execution. 


2.6.99.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 


164 


ENTER   ROUTINE 
LSTSQS  FROM 
1ACRO 


I 


DEFINE  FUNCTION 
LINF(l)=(l*(l-l))/2 


I 


INDEX  FOR 
SUBNETWORKS 
1111=  1,  IIJJ 


SUBNETWORK  MAPPING  SUBROUTINE,  LINKD 

RETURNS  LIST  OF  1ST  GENERATION  CONTROLLER 
NUMBERS  IN  SUBNETWORK  (INTS),  1ST  GENERATION 
DATA  LINK  NUMBERS  ASSOCIATED  WITH  EACH  SUB- 
NETWORK LINK  (LINK),  LIST  OF  2ND  UPSTREAM  SUB- 
NETWORK NODE  NUMBERS  ASSOCIATED  WITH  EACH 
LINK  IN  THE  SUBNETWORK  (LNDEX),  LIST  OF  1ST 
UPSTREAM  SUBNETWORK  NODE  NUMBERS  ASSOCIATED 
WITH  EACH  LINK  IN  THE  SUBNETWORK  (INDEX),  LIST 
OF  1ST  DOWNSTREAM  SUBNETWORK  NODE  NUMBERS 
ASSOCIATED  WITH  EACH  LINK  IN  THE  SUBNETWORK 
(JNDEX),  LIST  OF  2ND  DOWNSTREAM  SUBNETWORK 
NODE  NUMBERS  ASSOCIATED  WITH  EACH  LINK  IN 
THE  SUBNETWORK  (ONDEX),  LIST  OF  2ND  GENERA- 
TION LINK  NUMBERS  ASSOCIATED  WITH  EACH  LINK 
IN  THE  SUBNETWORK  (LINK2) 


Routine  LSTSQS 
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NO 


SET  NUMBER  OF 

INTERSECTION 

NN 

MDATAQ,  1111)-=  I 


NO 


YES 


INDEX  FOR  1ST  GENERATION  LINK  NUMBERS 
1=1,175   • 


SET  NUMBER  OF  LINKS 
LINKS       =  I 
MDATA(2,  1111)=  I 


1=1  +  1 


NO 


Routine  LSTSQS  (Continued) 
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1 


IR1\IDM=  827 


INITIALIZE  LINK  INDEX 
KMIN=  1 
IX  =  0 


INDEX  FOR  1ST  GENERATION 
LINK  NUMBERS 
IS=  1,  LNKS 


YES 


KD 


INCREMENT  COUNTER  INDICATING  NUMBER 
OF  ACTUAL  LINKS  WITH  UPSTREAM  AND 
DOWNSTREAM  NODES 
IX  =  IX  +  1 


SET  ALL  NODES  FOR  THIS  LINK 
LNDEX  (IX)  =  LNDEX  (IS)  2ND  UPSTREAM 
INDEX  (IX)  =  INDEX  (IS)  1ST  UPSTREAM 
JNDEX  (IX)  =  JNDEX  (IS)  1ST  DOWNSTREAM 
ONDEX  (IX)  =  ONDEX  (IS)  2ND  DOWNSTREAM 


Routine  LSTSQS  (Continued) 
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1 


LINK  DATA  SUBROUTINE  LDM0E5 


RETURNS  VOLUME  (I  VOL),  SPEED 
(I  SPEED)  AND  OCCUPANCY  (I  OCC) 


I 


INTERSECTION  STATUS 
SUBROUTINE  LLGTHF 


RETURNS  LENGTH  (LGTH), 
PHASE  (IFZE),  AND  WEIGHTING 
FACTOR  (LIMF) 


I 


CONVERT  LENGTH  TO  A  REAL  NUMBER 
FLNGT  =  LGTH 


IVOL       =  PV  (LINK(IS)) 
ISPEED=  PV(LINKdS)) 


CONVERT  SPEED,  VOLUME,  AND  WEIGHTING 
FACTOR  TO  REAL  NUMBERS 

VOL  =  IVOL      FLIM=LIMF      SPEED  =  ISPEED 


I 


COMPUTE  WEIGHTING  FACTOR  FOR  LINK  IX 
AMN  (IX)  =  XD*VOL*FLIM 
COMPUTE  IDEALIZED  OFFSET  FOR  LINK  IX 
FIDEAL(IX)  =  FLNGT/SPEED 


I 


IFAZE(IX)=IFZE 


I 


Routine  LSTSQS  (Continued) 
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INCREMENT  LINK 

INDEX 

IS=  IS  +  1 


SET  MAXIMUM  NUMBER  OF  ACTUAL  LINKS 
KMAX=  IX 


I 


INITIALIZE  INDICATOR  FOR  NUMBER  OF 

ACTUAL  LINKS  PLUS  CONTINUITY  LINKS 

IXX  =  IX 

RESET  COUNTER  FOR  NUMBER  OF 

ACTUAL  LINKS 

IX  =  0 

©■ 


INDEX  FOR  1ST  GENERATION  LINK 
NUMBERS  FOR  OUTER  LOOP 
IS=  1,  LNKS 


I 


INITIALIZE  INDEX  INDICATING  NUMBER  OF 
LINKS  WITH  EXISTING  UPSTREAM  AND 
DOWNSTREAM  NODES  IN  INNERLOOP 
IIX=0 


IS  1ST  UPSTREAM  NODE 
OR  1ST  DOWNSTREAM  NODE 
FOR  THIS  LINK  EQUAL  TO  ZERO  ? 
INDEX(IS)=0 
JNDEX(IS)=  0 


YES 


<D 


NO 


INCREMENT  COUNTER  INDICATING 
NUMBER  OF  ACTUAL  LINKS  IN  OUTER  LOOP 
IX  =  IX  +  1 


I 


Routine  LSTSQS  (Continued) 
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2 


INDEX  FOR  1ST  GENERATION  LINK  NUMBERS 
FOR  INNER  LOOP 
ISS  =  1,  LNKS 


INCREMENT  COUNTER  INDICATING  NUMBER 
OF  ACTUAL  LINKS  IN  INNER  LOOP 
IIX=  IIX  +  1 


L=  LNDEX(IX) 
INDEX(IIX) 


o 


1=  INDEX(IX) 
J  =  JNDEX(IIX) 


<D 


J  =  JNDEX  (IX) 
0=  ONDEX  (IIX) 


<D 


INCREMENT  COUNTER  INDICATING 
NUMBER  OF  ACTUAL  LINKS  PLUS 
CONTINUITY  LINKS 
IXX  =  IXX  +  1 


■© 


Routine  LSTSQS  (Continued) 
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£ 


COMPUTE  WEIGHTING  FACTOR  FOR  THIS 

CONTINUITY  LINK 

AMN(IXX)  =  (AMN(IX)  +  AMN(IXX))/2 


COMPUTE  IDEALIZED  OFFSET  FOR  THIS  CONTINUITY 
LINK  FIDEAL  IXX)  =  FIDEAL  (IX)  +  FIDEAL  (MX) 


STORE  UPSTREAM  AND  DOWNSTREAM  NODES  FOR 

THIS  CONTINUITY  LINK 

INDEX  (IXX)  =  LNDEX  (IX)   JNDEX  (IXX)  =  JNDEX  (IX) 


STORE  CONTINUITY  LINK  NUMBER  IXX  ASSOCIATED 
WITH  ACTUAL  LINK  NUMBERS  IX  AND  MX 
MO  (IIX)=  IXX      LN  (IX)   =  IXX 


STORE  WEIGHTING  FACTOR  FOR  CONTINUITY  LINK 
NUMBER  IXX  ASSOCIATED  WITH  ACTUAL  LINK 
NUMBERS  IX  AND  IIX 
ALN  (IX)  =  AMN(IXX) 
AMO  (IIX)=  AMN(IXX) 


Routine  LSTSQS  (Continued) 
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1 


SET  MAXIMUM  NUMBER  OF  ACTUAL 
LINKS  PLUS  CONTINUITY  LINKS 
KMAXS=IXX      KMAX  =  KMAXS 
JCYCLE=  NCYCLE(IIII) 


I 


GET  THE  CANDIDATE  CYCLE  LENGTHS 

CURV  (JJJJ)  AND  SET  THE  INITIAL  MINIMUM 

LEAST  SQUARES  RESIDUES  ADL  (JJJJ)  FOR 

SUBNETWORK  -  Mil 

CURV  (JJJJ)  =  CYC  (llll,  JJJJ) 

ADL  (JJJJ)  =  lOOOOOO 


I 


ASSEMBLE  A  MATRIX  OF  WEIGHTING 
FACTORS  IN  THE  ARRAY  E. 


I 


INVERT  THE  MATRIX  BY  FORMING  A  TRIANGULAR 
MATRIX,  INVERTING  THE  LOWER  TRIANGLE,  AND 
CALCULATING  THE  ELEMENTS  OF  THE  INVERSE 
MATRIX.   THE  INVERTED  MATRIX  IS  STORED  IN 
THE  ARRAY  E. 


I 


FILL  LAST  MATRIX  ROW  WITH  ZEROS 


I 


COMPUTE  ARRAY  G  FROM  THE  INVERSE  OF  LINK 
WEIGHTING  FACTOR  MATRIX  AND  LINK  WEIGHT- 
ING FACTORS 

G(K)  =  G(K)  +  E  (LLJ)  -  E(LLL) 


Routine  LSTSQS  (Continued) 
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£ 


SET  WEIGHTING  TERMS  FOR  EACH  CONTINUITY  LINK  ASSOCIATED 
WITH  LINK  K  WHERE  K  =  KMIN,  KMAX 

ALPHA  (K)  =  G(K)  FOR  LINK  LJ 

BETA(K)     =  G(K)  FOR  LINK  IJ 

GAMMA  (K)=  G(K)  FOR  LINK  10 

DELTA  (K)  =  G(K)  FOR  LINK  LO 

EPSIL(K)   =  G(K)  FOR  LINK  LI 

PHI  (K)        =  G(K)  FOR  LINK  JO 


I 


COMPUTE  DECISION  VECTOR  X  FOR  EACH  ACTUAL  LINK  K  WHERE 

K  =  KMIN,  KMAX 

X(K)=(ALN  (K)* 

(ALN  (K)  * 

ALPHA  (K)  + 

1 

+  AMN  (K) 

*(AMN(K) 

*  BETA  (K)  + 

1 

+  AMO(K) 

*  (AMO  (K) 

*  GAMMA  (K) 

+  1 

+  ALN  (K) 

*AMO(K) 

*  (DELTA  (K) 

-EPSIL(K)- 

PH  1(10  + BETA  (K)) 

+  ALN(K)- 

k  AMN(K)*  (ALPHA  (K)  ■ 

■  EPSIL(K) 

+  BETA  (K)) 

+  AMO(K) 

*AMN(K)" 

*(BETA(K)  + 

GAMMA  (K) 

-  PHKK)))*  .5 

I 


SET  FLAT  TO  INDICATE  OPTIMAL  CYCLE  LENGTH  NOT  YET  FOUND 
IOP=  0 


©■ 


I 


INDEX  FOR  CANDIDATE  SIGNAL  CYCLE  LENGTHS 
JJJ  =  1,  NCYCLE(IIII) 


COMPUTE  LINK  K  IDEALIZED  OFFSETS  IN 
FRACTION  OF  CYCLE  LENGTH  WHERE  K  =  KMIN,  KMAX 
SCL=  1./CURVUJJ) 
Q(l<)=  SCL(  FIDEAL(K) 


SET  CRITERION  THRESHOLD  FOR  SECOND  ORDER  MINIMUM  (EPSI) 
MODIFY  SEED  FOR  RANDOM  NUMBER  GENERATOR  (IRNDM),  AND 
INITIALIZE  MINIMUM  DELAY  (DLYMN) 

EPSI  =.00001       IRNDM=  2*  IRNDM  +  1      DLYMN  =  I.OE  +  35 


Routine  LSTSQS  (Continued) 
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INITIALIZE  FOR  RANDOM  START 

FOR  SEARCH  FOR  SECOND  ORDER  MINIMUM 


I 


INITIAL  VALUE  OF  OFFSETS  FROM  RANDOM 
NUMBER  GENERATOR  AT  START  OF  GAME 
FOR  EACH  INTERSECTION 
TD(I)  =  RNDMX  *  2.384185791015625  E-7 


I 


INITIALIZE  ARRAY  OF  RESIDUES  OF  THE 
OPTIMAL  OFFSET  DIFFERENCE  AND  IDEAL 
OFFSET  FOR  EACH  LINK 
T1  =  QQ(KK)  +  TD(I)-TD(J) 
QQ(KK)=  Tl  WHERE  -.5  <  Tl    <  .5 


I 


ARRAY  OF  WEIGHTED  OFFSET  DIFFERENCES 
FOR  EACH  LINK  PAIR 
C(l)=  C(D-  AMN(K)*T1 
C(J)=  C(J)  +  AMN(K)*T1 


©■ 


I 


IN  A  GIVEN  FIXED  REGION 

INITIAL  OPTIMAL  SIGNAL  OFFSETS 

FOR  EACH  INTERSECTION 

D(l)=  D(l)  +  E(LiJ)*C(J) 

TD(I)=  TD(I)*D(I) 


ARRAY  OF  RESIDUES  OF  THE 
OPTIMAL  OFFSET  DIFFERENCE  AND 
IDEAL  OFFSET  FOR  EACH  LINK 
IN  A  GIVEN  REGION  (LOCAL) 
QQ(KK)=  QQ(KK)  +  D(I)  -  D(J) 


I 


SEARCH  FOR  NEIGHBORING  OPTIMA 


Routine  LSTSQS  (Continued) 
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2. 


COMPUTE  DECISION  VECTOR  Y  FOR  EACH 
OPTIMIZABLE  LINK  KK  IN  THE  SUBNETWORK 
Y(KK)  =  ALN(KK)*  QQ(ILN) 
+  AMN(KK)*QQ(KK) 
+  AMO(KK)*QQ(IMO) 


COMPUTE  CHANGE  IN  DELAY  DUE 
TO  SHIFT  TO  NEIGHBORING  REGION 

/IENSIONAL  SPACE 
X(KK)-  IY(KK)I 


IN  N2  Dl 


CHDLY^ 


YES 


SET  INDEX  FOR  GRADIENT  CHANGE 
WHICH  PROVIDES  MAXIMUM 
REDUCTION  AND  SET  DELAY 
yiAXK  =  K 

;hmax  =  chdly 


YES 


*© 


>     CHMAXL=CHMAX 


YES 


BITERS  =  ITERATION  NUMBER  MINUS  50 
BITERS  =  ITERS  -50 


Routine  LSTSQS  (Continued) 
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0 


CITERS  =  ITERS  -100 
DITERS=  ITERS  -150 


YES 


YES 


YES 


FOR  THE  LINK  WHICH  RESULTED 

IN  THE  GREATEST  CHANGE  IN  DELAY 

DUE  TO  GRADIENT  CHANGE 


© 


COMPUTER  FINAL  DELAY 

OVER  ALL  LINKS 

DELAY  =  DELAY  =  AMN(K)  *  QQ(K)  *  QQ(K) 


® 


Routine  LSTSQS  (Continued) 
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© 


COMPUTE  NEW  RESIDUE  OF  THE 

OPTIMAL  OFFSET  DIFFERENCE  AND 

IDEALOFFSET 

QQ(KK)  =  QQ(KK)  +  EN  WHERE 

EN=  ±1.0 


COMPUTE  WEIGHTING  FACTOR  FOR  EACH 
CONTINUITYLINK  ASSOCIATED  WITH 
LINK  K  LMNO=  1,  3 


AMNP=  AMN(IO*EN 
I  =  INDEX(K) 

J  =  JNDEX(K) 


YES 


AMNP=  AMO(K)*EN 
I  =  INDEX(K) 

J  =  ONDEX(K) 


AMNP  =  ALN(K)*EN 
l=LNDEX(K)      J  =  JNDEX(K) 


FOR  ALL  INTERSECTIONS  COMPUTE  OFFSETS 
L=  1,  NN 


J  >  L 


YES 


LLJ=  LINF(J)  +  L 


LLJ=  LINRD  +  J 


>  L 


YES 


LLI  =  LINF(1)  +  L 


LLI=  LINF(L)+  1 


S  =  AMNP*(E(LLJ)  -  E(LLD) 

DD(L)=  DD(L)  +  S 
D(L)    =  DD(L) 
TD(L)=  TD(L)+  S 


Routine  LSTSQS  (Continued) 
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MODULO  THE  INTERSECTION  OFFSETS  (D(D) 
0  <  D(l)<  1 


I 


Tl  =  DIFFERENCE  BETWEEN  DELAY  FOR  THIS  GAME  AND  MINIMUM  DELAY 

Tl  =  (DELAY  -DLYMN)/DLYMN  OR 

TI  =  DELAY  -DLYMN  IF  DLYMN  <  0.00001 


Routine  LSTSQS  (Continued) 
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Routine  LSTSQS  (Continued) 
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SET  THE  CYCLE  LENGTH 

WHICH  GAVE  THE  SMALLEST  DELAY 

CL  =  CURV  (J) 


I 


SET  CYCLE  LENGTH  FOR 
ONE  MORE  SET  OF  RANDOM 
GAMES  TO  RESTORE  THE 
OPTIMAL  OFFSETS  FOR  THIS 
CYCLE  LENGTH 
CURV(l)  =  CL 


I 


SET  FLAG  DENOTING  OPTIMAL 
CYCLE  LENGTH  HAS  BEEN  FOUND 
IOP  =  1 


RESET  IOP  FLAG  FOR 
NEXT  SUBNETWORK 
IOP  =  0 


I 


IS  =  INDEX  FOR  INTERSECTIONS 


INTERSECTION  ENTRY  DATA,  MACA 


NUMBER  OF  INTERSECTION  PHASES  (IINT(4)), 
A-PHASE  SPLIT  (IINT(6)),  AND  C-PHASE 
SPLIT  (IINT(7)) 


I 


A-PHASE  OFFSET  10(2)  AND 

B-PHASE  OFFSET  10(3)  FOR  OUTPUT  ARRAY 

10(2)  =  AE  (IS)*CL  +  0.5 

10(1)  =  CL 

10(3)  =  1111 


i 


Routine  LSTSQS  (Continued) 
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LOAD  OFFSETS 
AND  VARIABLE 
INTERVALS 


OPTIMIZE  OPERATOR 
ENTERED  OFFSETS 


(RETURN  TO^\ 
MACRO       J 


Routine  LSTSQS  (Concluded) 
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2.6.100    SUBROUTINE  CORRECT 

2.6.100.1  Purpose 

Subroutine  CORRECT  computes  a  correction  factor  for  phase  connection 
other  than  an  A-phase  to  A-phase  link. 

2.6.100.2  Relationship 

CORRECT  is  called  from  Routine  LSTSQS  by: 

CALL  CORRECT  (ITB,    ICORR) 

where: 

ITB(l)  denotes  downstream  intersection 
ITB(2)  denotes  upstream  intersection 
ITB(3)  denotes  head/tail  phase  code 
ITB(4)  denotes  cycle  length 
ICORR  returns  the  correction  factor 

2.6.100.3  Assumption 
None. 

2.6.100.4  Approach 

Subroutine  CORRECT  determines  the  A  and  B  phase  durations  of  two 
adjacent  nodes.     The  idealized  offset  is  then  adjusted  to  reflect  the  phase 
connection. 

2.6.100.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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f    ENTER   ) 

IPH(2)  = 

MOD  (ITB(3),  3) 

IPH(1)=  ITB(3)/3 


IDX  =  IT6CI) 
JDX  =  4*(!DX-1)  +  1 
ICL  =  2*ITB(4) 
-SMTABCUDX  +  3) 
INX=  2((IDX-1)  +  1 


T 


AA  =  ISP*ICL 
APH  =  AA/IOOO, 
+  SMTABC  (JDX) 
FAC(I)  =  APH 


Subroutine  CORRECT 
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1 


BA=  SPLIT  (INX  +  1) 

BB=1CL-(AA+BA) 

1000 
BPH  =  BB  +  SMTABC  (JDK  +  1) 
FAC(I)  =  FAC(I)  +  BPH 


1 


FAC(I)  =  0 


YES 


FAC(I)=  IDLT 


FAC(I)  = 
FAC(I)+  APH 


Subroutine  CORRECT  (Continued) 
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ICORR  = 

FAC(2)  -  FAC(l) 


(return} 


Subroutine  CORRECT  (Concluded) 
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2.  6.  101    SUBROUTINE  LINKD 

2.  6. 101.  1    Purpose 

Subroutine  LINKD  compiles  a  list  of  first  generation  controller  numbers 
included  in  each  subnetwork,   and  the  subnetwork  node  numbers  associated 
with  each  link  in  the  subnetwork. 

2.6.101.2    Relationship 

This  subroutine  is  called  by: 

CALL  LINKD  (ISB,   INT,    LNK,    LNDX,  INDX,   JNDX,   KNDX,    LNK2) 

where: 

ISB  =  Input  argument  containing  list  of  group  numbers  in  subnetwork 

INT  =  Output  argument  giving  list  of  first  generation  controller 
numbers  in  subnetwork 

LNK  =  Output  argument  giving  first  generation  data  link  numbers 
associated  with  each  subnetwork  link 

LNDX  =  Output  argument  giving  list  of  second  upstream  node 
numbers  associated  with  each  link  in  the  subnetwork 

INDX  =  Output  argument  giving  list  of  first  upstream  subnetwork 
node  numbers  associated  with  each  link  in  the  subnetwork 

JNDX  =  Output  argument  giving  list  of  first  downstream  subnetwork 
node  numbers  associated  with  each  link  in  the  subnetwork 
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KNDX  =  Output  argument  giving  list  of  second  downstream  subnetwork 
node  numbers  associated  with  each  link  in  the  subnetwork 

LNK2  =  Output  argument  giving  list  of  second  generation  link  numbers 
associated  with  each  link  in  the  subnetwork 

Subroutine  LINKD  is  called  by  LSTSQS  and  returns  to  that  routine. 


None. 


2.6.101.4    Approach 

LINKD  accesses  the  arrays  ISB,  CGRPN,   NCIGRP,   LMNO,   GCNTSI,  GCONT, 
1GLNK,   N2GCL,  and  2GCLSI,  to  build  the  INT,    LNK,   LNDX,  INDX,   JNDX, 
KNDX,  and  LNK2  arrays  which  contain  subnetwork  node  and  link  numbers 
for  each  link  in  the  subnetwork. 


2.6.101.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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C  ENTER    J 

INT(l-45)=  0 
LNK(1-175)=  0 
LNDX(1-175)=0 
INDX(1-175)=0 
JNDX(1-175)=  0 
KNDX(1-175)  =  0 
LNK2(1-175)  =  0 
ICT=  0 


I 


1=  1 


YES 


<D 


ISBN(I)=  ISB(I) 
ICT=  1 


NO 


Subroutine  LINKD 
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1  =  1  +  1 

11=  ISBN(I) 
IGN=  NCIGRP(II) 
1ST  =  GCNTSKI!) 

J  =  l 


11=  IST+J 
IND  =  GCONT(II) 

JST  =  JST  +  1 
INT(JST)=  IND 
J  =  J  +  1 


PLACE  NODES  OF 
SUBNETWORK  INTO 
TABLE  INT 


Subroutine  LBMKD  (Continued) 
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£ 


1=0 
NLK=0 


1=1  +  1 
11=  ISBN(I) 
LNKT=  N2GCL 
IST=  GC2LSI 


7 


11  =  ISB  +  J 
LN1  =  G1LNK(II) 
NLK=  NLK  +  1 
LNK(NLK)=  LNK1 
LNK2(NLK)  =  II 
JJ  =  4*(ll-l)  +  1 


I 


LNDX(NLK)=  LMNO(JJ) 
INDX(NLK)=  LMNOUJ  +  1) 
JNDX(NLK)=  LMNOUJ+  2) 
KNDX(NLK)=  LMNOUJ+3) 
J  =  J  +  1 


Subroutine  LINKD  (Continued) 
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£ 


i=  i 


CALL  SRCH 
(LNDX,  INT,  JST) 


CALL  SRCH 
(INDX,  INT,  JST) 


CALL  SRCH 
(JNDX,  INT,  JST) 


CALL  SRCH 
(KNDX,  INT,  JST) 


1=1  +  1 


NO 


Subroutine  LINKD  (Concluded) 
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2.6.102    SUBROUTINE  LLGTHF 

2.6.102.1  Purpose 

Subroutine  LLGTHF  retrieves  the  link  length,    phasing  code,    and  link  im- 
portance factor  for  a  specified  link  from  the  LINKL,    LINKHT,    IALPH, 
and  JALPH  arrays. 

2.6.102.2  Relationship 

The  calling  sequence  for  this  subroutine  is  as  follows: 
CALL  LLGTHF  (L2,    LGT,  IFZ,    LIM,   MSP) 

where: 

L2  =  Input  argument  (second  generation  link  number) 

LGT  =  Link  length  in  feet 

IFZ  =  Link  phasing  code 

LIM  =  Link  importance  factor 

MSP  =  Speed  Multiplier  Factor. 

This  subroutine  is  called  by  Routine  LSTSQS  and  returns  to  that  routine. 

2.6.102.3  Assumption 
None. 
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2.6.102.4    Approach 

The  elements  of  second  generation  data  base  arrays  LENGTH,    LFZE,   and 
IMP  associated  with  Link  (L2)  are  accessed  and  loaded  .into  the  return  argu- 
ments. 


2.6.102.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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f  ENTER  J 

g 

HALF 

(LINKL,  L2,  LGT) 


1 


GETBY 

(LINKHT,  L2f  IFZ) 


I 


GETBY 

(IALPH,  L2,  LIM) 


I 


GETBY 

(JALPH,  L2,  MM) 


I 


MSP  = 
MULT/100 


f  return) 


Subroutine  LLGTHF 
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2.  6.  103    SUBROUTINE  LOADAT 

2.6.103.1  Purpose 

This  subroutine  loads  offsets  and  variable  intervals,    and  provides  for  the 
maintenance  of  operator  selected  splits  and  offsets. 

2.6.103.2  Relationship 

LOADAT  is  called  from  LSTSQS  by: 

CALL  LOADAT  (INT,   JTB) 

where: 

INT  denotes  the  intersection 

JTB(l)  denptes  cycle  length 

JTB(2)  denotes  offset 

JTB(3)  denotes  number  of  subnetworks. 

2.6.103.3  Assumption 
None. 

2.6.103.4  Approach 

Subroutine  LOADAT  stores  the  working  offset  computed  by  LSTSQS.     It  then 
computes  the  phase  durations  based  on  the  split  and  the  controller  type 
(K3PM),    and  uses  these  durations  to  determine  the  durations  of  the  variable 
intervals.     LOADAT  stores  these  variable  interval  durations  and  returns  to 
the  calling  routine. 

2.6.103.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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YES 


ICL  =  2*JTB(1) 

IAOF=  2*JTB(2) 

ICOF=0 

IN2=  2*(INT-  1)  +  1 

IN4=4*(INT-  1)  +  1 

IN16  =  16*(INT-  1)  +  1 

IBS  =  BIAS  (INT) 


-© 


>  128 

v  YES 

«^ 

IBS  =  128  -  IBS 

7           S 

^^ 

Yno 

Subroutine  LOADAT 
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© 


ITC  =  8 

ICL=KCL(1TC) 


3 


ITC  =  ITC  -1 
ENTCYC=  ITC  (INT) 


WRKOF  =  JJ  (INT) 
K3P=K3PHI  (INT) 
IASP=  SPLIT  (IN2) 
IFX  =  SMTABC  (IN4  +  3) 
JCL=  ICL-  IFX 


7 


IAPH  =  yCLj^LSOO) 

IAV  =  VRINT0N4) 
JAV=  IN16  +  IAV 
IAFX=  SABC(IN4) 


*■  JJ  =  JJ  +  ICL 


JJ  =  JJ  -  ICL 


Subroutine  LOADAT  (Continued) 
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IAPV  =  IAPH  -  IAFX 
ENJNT  =  IAPV  (JAV) 
IFNA=  SMTABC  (IN4) 
IAD=  IAPH  +  IFNA 
IBOF=  IAOF  +  IAD 


Subroutine  LOADAT  (Continued) 
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© 


IBPH  =  JCL-  IAPH 
IBFX  =  SABC  (IN4  +  1) 


B  PHASE  DURATION 


K3P  =   0 

TYPE  0  CONTROLLER 

(TWO  NON-OVERLAPPING 
PHASES) 


YES 


IBFX=  SUMABAC  (INT) 


IBPV=  IBPH  -  IBFX 
IBV=  VRINTUN4  +  1) 
JBV=  IN16=  IBV 


ENINT=  IBPV(JBV) 


NO 


I0CIN4)  =  IAOF 
I0(IN4  +  1)=  IBOF 
I0(IN4  +  2)=  ICOF 
IEA(INT)  =  JTB(3) 


WRKOF(INT)=  IBOF 


K3P  = 


1? 

v  YES 

JCL  =  JCL  -  ICPH 

'no 

RETURN 


Subroutine  LOADAT  (Continued) 
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ICSP=  SPLIT  (IN2  +  1) 
irDu     (JCL*lCSP  +  500) 
,UKM_  1000 

IBPH  =  JCL-  IAPH  -  ICPH 
ICV=  VRINT0N4  +  2) 


I 


K3P  =    1 

TYPE  I    CONTROLLER 
(THREE  NON-OVERLAPPING 
PHASES  ) 


JCV=  IN16=  ICV 
ICFX=  SABC(IN4  + 
ICPV  =  ICPH  -  ICFX 


2) 


IBV  =  VRINT(IN4  +  1) 
JBV=  IN16  +  IBV 
IBFX=  SABC(IN4  +  1) 
IBPV=  IBPH  -  IBFX 
ENINT=  IBPV(JBV) 


I 


IFNS=  SUMTABC  (IN4  +  1) 
ICOF=  IBOF  +  IBPH+  IFNB 


YES 


ICOF=  ICOF-  ICL 


Subroutine  LOADAT  (Continued) 
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IDLT  =  DELT(INT) 


K3P  =    2 

TYPE  2  CONTROLLER 

(A  PHASE/CPHASE  OVERLAP) 


IDLT 


>128s 

VYES 

?     ^/ 

IDLT  -  l^o  -  iul  i 

Oo 

ICOF<0 


YES 


ICOF=  ICOF  +  ICL 


ICOF 


*  ICL 

?      y 

vYES 

IC0F=  ICOF-  ICL 

Y^iO 

IDLT*  0 


YES 


i/do) 


NO 


IBOF=  ICOF  +  IAD 


Subroutine  LOADAT  (Continued) 
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YES 


IDLT  =  128  -  IDLT 
IBOF=  IBOF-  IDLT 


YES 


K3P  =   3 

TYPE  3  CONTROLLER 

(B  PHASE/  C  PHASE  OVERLAP) 


ICOF=  IBOF  +  IDLT 


IBOF=  IBOF  +  ICL 


Subroutine  LOADAT  (Continued) 
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<s> 


©• 


IHOLD=  1 

ICON  =  JTB(IHOLD) 

ITEM  =  HLDOFF(ICON) 


YES 


JJTEM  =  ((ICON  -  1)*4)  -  1 
ITEM-  CYCLC  (ICON) 
IHCYC=  0.5*KCL  (ITEM) 
MIN=  10  000 
KTEM  =  1 


3 


JSUM=  0 
LTEM=  1 


3 


ISURCON 

=  SURCON  (KTEM  +  LTEM) 


YES 


WreturnJ 


*© 


Subroutine  LOADAT  (Continued) 
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© 


NOFF  =  ENOFF(ISURCN) 
NOFF=  NOFF-(KTEM-l) 


© 


y.    NOFF    \ 

YES 

v.       >  IHCYC      J 

NOFF- NOFF-nolu i tMJ 

\v     ?      y^ 

yNO 

YES 

\  <  ihcyc  y 

NOFF = NOFF + 
KCL(ITEM) 

T'nO 

ISUM=  ISUM + NOFF 
LTEM  =  LTEM  +  1 

<M_TEM  >4      J> 

NO 

i 

YES 


ISUM 
*  MIN 
? 


YES 


NO 


MIN  =  ISUM 
MO=  KTEM-  1 


1 


KTEM=  KTEM  +  1 


HI  H- 


NO    /  KTEM  \    YES 

>KCL(ITEM). 
? 


■© 


Subroutine  LOADAT  (Continued) 


204 


Subroutine  LOADAT  (Concluded) 
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2.6.104    SUBROUTINE  SRCH 

2.6.104.1  Purpose 

This  subroutine  scans  a  table  specified  by  the  second  input  argument  to  find 
an  element  equal  to  the  first  input  argument. 

2.6.104.2  Relationship 

Subroutine  SRCH  is  called  from  Subroutine  LINKD  by: 

CALL  SRCH  (J,    II,    N) 

where: 

J  =  the  element  to  be  searched  for 
II  =  the  table  to  be  searched 

N  =  returns  to  the  location  in  the  table  where  J  can  be  found. 
If  J  is  not  found,   N  is  set  to  zero. 

2.6.104.3  Assumption 
None. 

2.6.104.4  Approach 

The  specified  table  is  scanned  until  the  element  desired  is  found.     When 
found,    SRCH  returns  to  LINKD  with  the  third  parameter  equal  to  the  ele- 
ments location  within  the  table.     If  the  element  is  not  found  N  is  set  to  zero. 

2.6.104.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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Q  ENTER   J 

,r 

K=  0,1=0 

^ 

w 

Subroutine  SRCH 
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2.6.105    ROUTINE  INTFC 

2.  6.  105.  1     Purpose 

The  purpose  of  Routine  INTFC  is  to  optimally  interface  each  adjacent  sub- 
network established  by  Routine  RTSND.     The  interface  is  optimized  by 
insuring  that  vehicle  delay  is  minimized  on  the  interface  links  connecting 
two  subnetworks  at  the  midpoint  of  the  ensuing  5MINCT  time  interval. 

2.6.105.2  Relationship 

Routine  INTFC  is  called  from  Routine   MACRO  by: 

CALL  INTFC 
INTFC  calls  subroutines  GETINT,    LDMOE5,    LLGTHF;    MACA,    and  PERI05, 

2.6.105.3  Assumption 


None. 


2.  6.  105. 4    Approach 

The  methodology  used  is  to  determine  an  interface  parameter  which,   when 
added  to  the  newly  computed  entry  offsets  of  one  of  the  subnetworks  on  an 
interface  boundary  will  minimize  delay  on  the  interface  links. 

A  quadratic  delay  function,   similar  to  the  form  in  Routine  LSTSQS  is 
assumed  for  the  links  connecting  intersections  between  two  subnetworks: 


\fr  a..  [~io..  +o.  -  o.  +q1 
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where 

A. .  =  link  weighting  factor 

IO..  =  ideal  offset  on  link  ii 
i]  j 

O.  =  entry  offsets  of  the  interface  intersections  in  subnetwork  1 

O.  =  entry  offsets  of  the  interface  intersections  in  subnetwork  2 
3  J 

Q  =  interface  parameter 

D  =  total  delay  over  all  links. 

However,  this  delay  function  has  only  one  independent  variable,  Q.     O.  and 
O.  have  been  determined  by  the  network  optimization  routine.     The  function 
is  minimized  by  setting  the  partial  derivatives  with  respect  to  Q  equal  to 
zero  yielding: 

W  A.. A.. 
Q  =    1  3       1J     1J  (1) 

1  J      ij 

where  the  summation  is  over  all  links,  A.  is  the   link  weighting  factor,  and 

A..  =  IO..  +  O.  -  O.  (2) 

13  13         i         3 

Routine  INTFC  program  coding  is  described  below.     When  the  interface  flags 
INTF  indicate  computations  are  to  be  performed  for  the  interface,   the  entry 
cycle  length  and  time  in  the  cycle  are  determined  for  the  two  contiguous  sub- 
networks.    The  times  in  the  cycle  time  counters  of  Routine  N  are  extrapolated 
to  the  midpoint  of  the  ensuing  15 -minute  control  period. 
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ECYC  =  ECYC  +  450. 
ECYC2  =  ECYC2  +  450. 

where  ECYC  is  the  current  time  on  cycle  time  counter  of  Routine  N  for  the  entry 
cycle  length  in  the  cycle  of  subnetwork  1  and  ECYC 2  is  the  current  time  on 
the  cycle  time  counter  of  Routine  N  for  the  entry  cycle  length  of  subnetwork 
2.     These  times  are  then  normalized  by  their  respective  cycle  lengths. 

Data  are  obtained  for  each  of  the  interfacing  intersections.     The  entry  off- 
sets for  the  intersections  in  subnetwork  2  (XOFS2)  and  in  subnetwork  1 
(XOFS)  are  obtained.     A  link  weighting  factor  A(K)  and  an  idealized  offset 
ZO(K)  for  the  link  approaching  each  intersection  in  subnetwork  1  are  ob- 
tained.   Then  the  A.,  of  Equation  (2)  are  computed  for  the  links  approaching 
subnetwork  1  by  the  following  computations. 

XOFSP  =  ZO(K)  -  XOFS 
XOFSP  =  XOFSP  -  ECYC 
DEL  =  XOFSP  +  XOFS2 

Then  the  weighting  factors  A(K)  and  the  weighted  offset  differences  are 
added  to  the  sums  for  the  interface. 

SUMADL  =  SUMADL  +  A(K)  *  DEL 
SUMA  =  SUMA  +  A(K) 

The  A.,  for  a  link  departing  subnetwork  1  is  then  computed  as 

XOFSP  =  ZO(K)  +  XOFS 
XOFSP  =  XOFSP  -  ECYC 
DEL  = XOFSP  -  XOFS2 
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and  A(K)  and  A(K)  *  DEL  are  added  to  SUMA  and  SUMADL,   respectively. 

When  all  the  A. .  for  all  the  interfacing  links  have  been  processed  the  inter- 
face  parameter  is  computed. 

Q  =  SUMDEL/SUMA 

Using  the  index  for  INTF  flag  the  core  group(s)  in  subnetwork  2  are  ob- 
tained.    Then  the  intersection  numbers  in  subnetwork  2  are  obtained  and  Q 
is  added  to  the  entry  offset  of  each  intersection. 


2.6.105.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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c 


ENTER  SUBROUTINE 
INTFC  FROM  MACRO 


} 


1=  INDEX  FOR  CORE 
GROUP  INTERFACE 


OBTAIN  CYCLE  LENGTH  AND  TIME 
IN  CYCLE  LENGTH  FOR  TWO 
INTERFACING  SUBNETWORKS 


ENTRY  CYCLE  SUBROUTINE,  ENTCYCL 


RETURNS  TIME  IN  CYCLE  AND  CYCLE 
LENGTH  FOR  INTERSECTION  INT  IN 
JT(I)  AND  JT(2) 


I 


CYCLE  LENGTH  (CYC)  AND 
TIME  IN  CYCLE  (ECYC) 
CYC  =  JT(2) 
ECYC  =  JT(1) 


J. 


EXTRAPOLATE  TIME  TO  MIDPOINT  OF 
SIGNAL  SETTING  DURATION 
ECYC=  ECYC  +450. 


MODULO  SUBROUTINE,  PERI05 


MODULO  EXTRAPOLATED  TIME 
0  s  ECYC  <  CYC 


OBTAIN  CYCLE  LENGTH  AND  TIME 
IN  CYCLE  LENGTH  FOR  THE  SECOND 
OF  THE  INTERFACING  SUBNETWORKS 


<D 


Routine  INTFC 
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2 


ENTRY  CYCLE  SUBROUTINE, 
ENTCYCL 


TIME  IN  CYCLE  AND  CYCLE 
LENGTH  FOR  INTERSECTION  INT 


EXTRAPOLATE  TIME  TO  MIDPOINT 
OF  SIGNAL  SETTING  DURATION 
ECYC2=  ECYC2  +  450. 


MODULO  SUBROUTINE,  PERI05 


MODULO  EXTRAPOLATED  TIME 
0  ^  ECYC2<  CYC  2 


I 


ZERO  SUM  OF  WEIGHTED  OFFSET 
DIFFERENCES  (SUMADL)  AND 
SUM  OF  WEIGHTING  FACTORS  (SUMA) 
SUMADL  =  0. 
'SUMA=  0. 


INDX  =  POINTER  FOR  INTERSECTION  ARRAY 
INDX=  NGRPKD-  1 


& 


J  =  INDEX  FOR  INTERSECTIONS 
AT  INTERFACE 


I 


OBTAIN  PATTERN  OFFSETS  FOR 

INTERFACING  INTERSECTIONS  IN  SECOND  SUBNETWORK 


I 


INDX=  INDX  +  1 


(!) 


Routine  INTFC  (Continued) 
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2 


GET  INTERSECTION  NUMBER 
INT  =  INTN20NDX) 


INTERSECTION  ENTRY  DATA,  MACA 


RETURNS  A-PHASE  OFFSET  (IRT(D)  AND 
B-PHASE  OFFSET  (IRT(2)) 


I 


A-PHASE  OFFSET 
XOFS2=  IRT(l) 


IS  INTERFACE 
LINKAGE  TO  INTER- 
SECTION B-PHASE?, 

I.E.,  IS  LPHSEJ 
IDX)  =  1J 


NO 


YES 


SUBTRACT  EXTRAPOLATED  TIME 
X0FS2  =  X0FS2-  ECYC2 


X0FS2=  IRT(2) 


MODULO  SUBROUTINE, 
PERI05 


MODULO  OFFSET 
0  <:XOFS2<CYC2 


COMPUTE  OFFSET  RELATIVE  TO  FIRST 
SUBNETWORK  CYCLE  LENGTH 


Routine  INTFC  (Continued) 
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YES 


NO    4 


OBTAIN  PATTERN  OFFSETS  FOR 
FIRST  SUBNETWORK 


I 


INTERSECTION  NUMBER 
INT=  INTN(INDX) 


INTERSECTION  ENTRY  DATA, 
MACA 


RETURN  IRT(1)=  A-PHASE  OFFSET, 
IRT(2)=  B-PHASE  OFFSET 


p  X0FS2=X0FS2-CYC 


A-PHASE  OFFSET 
XOFS=  IRT(l) 


YES 


XOFS=  IRT(2) 


INDP=  INDEX 
INDP=  INDX*2 


I 


K=  INDEX  FOR  LINKS, 
ENTERING  AND  DEPARTING 


Routine  INTFC  (Continued) 
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© 


INDL  =  LINK  INDEX 
INDL=  INDP  +  K 


LINK  NUMBER 
ILNK=  LINKN(INDL) 


LINK  DATA  SUBROUTINE, 
LDMOE 


RETURNS  LENGTH,  VOLUME,  SPEED, 
AND  WEIGHTING  FACTOR 


A(K)=  LINK  WEIGHTING  FACTOR 

ZO(K)=  IDEALIZED  OFFSET  =  DISTANCE/SPEED 

A(K)=  IRT(6)*IRT(2)*.001 

ZO(K)=  IRT(1)/IRT(5) 


/IS  LINK  A  ^S. 
EPART1NG  LINK?,  I.E>s 

>  YES 

ki. 

SET  PATTERN  OFFSET 
NEGATIVE 

ISK=2?             S> 

w 

XOFS=-XOFS 

NO 

^ 

^ 

SUBTRACT  PATTERN  OFFSET 
FROM  IDEALIZED  OFFSET 
AND  THEN  SUBTRACT  TIME 
INTO  THE  CYCLE 
XOFSP=ZO(K)-XOFS 
XOFSP  =  XOFSP-ECYC 


Routine  INTFC  (Continued) 
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£ 


COMPUTE  SUMS  FOR 
INTERFACE  PARAMETERS 


YES 


SET  SECOND  SUBNETWORK 
PATTERN  OFFSET  NEGATIVE 
XOFS2=-XOFS2 


ADD  FIRST  AND  SECOND 
SUBNETWORK  OFFSETS 
DEL  =  XOFSP=XOFS2 


SUM  WEIGHTED  OFFSET  DIFFERENCES 
AND  WEIGHTING  FACTORS 
SUMADL=  SUMADL+A(K)*DEL 
SUMA=  SUMA  +  AOO 


INTERFACE  PARAMETER 
Q=  SUMADL/SUMA 


DETERMINE  CORE  GROUP  NUMBER 
FOR  SECOND  SUBNETWORK 


ICORE(2)=  0 


6 


Routine  INTFC  (Continued) 

2.17 


YES 


CORE  GROUP=  2 
ICORE(l)=  2 


YES 


ICORE(2)=  3 


YES 


IC0RE(1)  =  3 


YES 


ICORE(l)=  5 


YES 


YES 


IC0RE(1)=4 


CORE  GROUP  FOR 
INTERFACE  5 
ICORE(l)=  6 


ICORE(2)=  4 


OBTAIN  INTERSECTION 
NUMBERS  IN  SECOND 
SUBNETWORK 


<D 


GET  INTERSECTIONS  SUBROUTINE, 
GETINT 


RETURNS  CORE  GROUP  INTERSECTION 
NUMBERS  IN  INX  (1...24) 


Routine  INTFC  (Continued) 
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Routine  INTFC  (Continued) 
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Routine  INTFC  (Concluded) 
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2.6.106    SUBROUTINE  GETINT 

2.6.106.1  Purpose 

Subroutine  GETINT  gets  the  numbers  of  all  the  intersections  in  a  particular 
core  group. 

2.6.106.2  Relationship 

This  subroutine  is  called  by: 

CALL  GETINT  (NGROUP,    J) 

where: 

NGROUP  =  Group  number 

J  =  User's  buffer  for  the  intersection  numbers 

Subroutine  GETINT  is  called  by  INTFC  and  TRNPAR. 

2.6.106.3  Assumption 
None. 

2.6.106.4  Approach 

The  major  intersections  for  the  particular  group  are  obtained  from  the 
SUBMAJ  array  and  loaded  into  the  J  array.     Next,    the  minor  intersections 
are  obtained  from  the  SUBMIN  array  and  loaded  immediately  following  the 
majors.     When  the  number  of  majors  or  minors  is  less  than  32,    a  zero  is 
stored  in  the  first  unused  location  to  signify  the  end  of  the  data. 

2.6.106.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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Q  ENTER   J 


CLEAR  ARRAY 
J 


Subroutine  GETINT 
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2.6.107    SUBROUTINE  MAC  A 

2.  6.  107.  1    Purpose 

Subroutine  MAC  A  extracts  from  the  data  base  the  A-,   B-,  and  C -Phase 
entry  offsets,  the  intersection  phase,  the  entry  cycle  length,  and  the  inter- 
section splits  for  a  particular  intersection. 

2.6.107.2  Relationship 

This  subroutine  is  called  by: 
CALL  MAC  A(  INT,  IE  NT) 

where: 

INT  =  Intersection  number 

IENT(1,  . .  . ,  3)  =  A-,  B-,   and  C -Phase  entry  offsets,   respectively 

IENT(4)  =  Intersection  phase 

IENT(5)  =  Entry  cycle  length 

IENT(6)  =  A-Phase  intersection  split 

IENT(7)  =  C -Phase  intersection  split 

Subroutine  MACA  is  called  by  LSTSQS. 

2.6.107.3  Assumption 
None. 
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2.6.107.4    Approach 

Given  an  intersection  number  (INT)  MACA  extracts  the  A-,   B,   and  C -Phase 
entry  offsets,  the  intersection  phase,  the  entry  cycle  length,  and  the  inter- 
section splits  from  the  data  base. 


2.6.107.5    Listing 

See  Volume  III,    2GCF  Routine  Listings 
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1  =  1  +  1 

JDX-  IDX  +  I 
IOUT(l)  =  IQ(JDX) 
1 


I0UT(5)  = 

KCKENTCYL(INT)) 

I0UT(4)=K3PHI(FI\IT) 
IDX  =  2*(II\IT-1) 


3 


1=1  +  1 
I0UTCI  +  S) 
=  SPLIT(JDX) 


Subroutine  MACA 
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2.6.108    SUBROUTINE  PERI05 

2.6. 108. 1  Purpose 

The  purpose  of  Subroutine  PERI05  is  to  modify  (add  or  subtract)  a  given  num- 
ber by  the  cycle  length,    keeping  the  result  within  a  specified  range. 

2.6.108.2  Relationship 

This  routine  is  called  by: 

CALL  PERI05  (X,    LLC,    ULC,    CY) 

where: 

X  =  Number  to  be  tested  and  modified 

LLC  =  Lower  limit  for  resultant  value  of  X 

ULC  =  Upper  limit  for  resultant  value  of  X 

CY  =  Value  by  which  X  is  modified  to  put  it  in  the  desired  range. 

Subroutine  PERI05  is  called  by  level  5  routines  LSTSQS,    TRNPAR,   and 
INTFC. 

2.6.108.3  Assumption 
None. 
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2.6.108.4  Approach 

The  value  of  X  is  compared  to  LLC.     If  X  is  smaller,    CY  is  added  to  X  until 
X  is  larger  than  LLC.     Then  X  is  compared  to  ULC.     If  X  is  larger,    -CY  is 
added  to  X  until  X  is  smaller.     Thus,    the  result  is  necessarily  within  the 
desired  range  only  when  CY  <.  ULC  -  LLC,    which  is  always  the  case. 

2.6.108.5  Listing      . 

See  Volume  III,    2GCF  Routine  Listings. 
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(  ENTER    ^ 


^ 

^ 

v  NO 

hi 

X  -  LLC^ 

X  =  X  +  CY 

>0? 

F 

Subroutine  PER  105 
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2.6.109    SUBROUTINE  PUTEOFF 

2.6.109.1  Purpose 

The  purpose  of  the  Subroutine  PUTEOFF  is  to  put  the  entry  offset  value 
into  the  ENTRYOFF  array. 

2.6.109.2  Relationship 

This  subroutine  is  called  by: 

CALL  PUTEOFF  (INT,    LOC) 

where: 

INT  =  Intersection  number 

LOC  =  Location  of  the  entry  offset  value 

Subroutine  PUTEOFF  is  called  by  INTFC  and  TRNPAR. 

2.6.109.3  Assumption 
None. 

2.6.109.4  Approach  e 

The  entry  offset  value  is  obtained  from  LOC.     After  being  scaled  Bl,    it  is 
placed  in  the  ENTRYOFF  array  for  the  specified  intersection  number. 

2.6.109.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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f  ENTER    J 

IP0SN  =  4*<INT-1) 

N=3 

1=1 


X^K3PHT     >> 
V       =0?        > 

JNO 

.YES 

N  = 

2 

10  (IPOSN  +  1) 
=  2*E0FFS(I) 
1=1  +  1 

NO 

C         >N? 

Lyes 


CY=KCL(ENTCYC  +  1) 
EOFFA=  2*E0FFS(1) 
+  BIAS  (INT) 


Subroutine  PUTEOFF 
230 


YES 


E0FFA  =  2*E0FFS(2) 


ENOFF(INT)=  EOFFA 


(return) 


Subroutine  PUTEOFF  (Concluded) 
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2.6.110    ROUTINE  TRNPAR 

2.6.110.1  Purpose 

Routine  TRNPAR  computes  a  transition  parameter  which  minimizes  the 
sum  total  of  all  intersection  offset  changes  when  implementing  a  new  signal 
pattern.     This,   in  turn,   minimizes  the  sum  total  of  the  time  durations  of  the 
transition  from  the  old  to  the  new  patterns. 

2.6.110.2  Relationship 

Routine  TRNPAR  is  called  from  Routine  MACRO  by: 

CALL  TRNPAR 
Routine  TRNPAR  calls  subroutines  MACA,    PERI05,    and  PUTEOFF. 

2.6.110.3  Assumption 
None. 

2.6.110.4  Approach 

The  principle  of  the  method  is  to  determine  a  parameter,   which,   when 
added  to  all  the  entry  offsets  will  minimize  the  sum  of  the  squares  of  the 
distance  the  phase  offsets  must  change  in  the  offset  space  to  arrive  at  the 
new  pattern.     This  parameter  is  the  average  value  of  the  differences  between 
the  working  offsets  and  the  entry  offsets.     In  the  subroutine  program,   data 
to  compute  the  transition  parameter  are  obtained  by  a  "DO  LOOP"  on  the 
core  groups.     For  each  "DO  LOOP"  iteration  the  major  and  minor  inter- 
sections in  the  core  group  are  obtained  and  in  a  "DO  LOOP"  on  the  inter- 
sections the  following  are  computed. 
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WRKOFP  =  WRKT  -  WRKOF 
ENTOFP  =  ENTT  -  ENTOF 
DELOFF  =  WRKOFP  -  ENTOFP 

Where  WRKOFP  is  the  delta  between  the  anticipated  time  in  the  working 
cycle  when  transition  will  occur  (10-second  lead  time)  and  the  working  off- 
set,  and  ENTOFP  is  a  similar  term  for  the  entry  intersection  pattern. 

The  transition  parameter  is  computed 
XMIN  =  SUMDEL/NOINT 

where  SUMDEL  is  the  sum  of  DELOFF  for  all  intersections  and  NOINT  is 
the  total  number  of  intersections. 

The  transition  parameter  is  added  to  the  entry  offset  of  all  the  intersections 
in  the  network.     A  "DO  LOOP"  on  the  intersection  numbers  within  a  "DO 
LOOP"  on  the  core  groups  is  used  in  the  procedure. 

2.6.  110.  5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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c 


ENTER   ROUTINE 
TRNPAR  FROM  MACRO 


I 


D 


©■ 


ZERO  NUMBER  OF  INTERSECTIONS  (NOINT) 
AND  SUM  OF  DIFFERENCES  BETWEEN 
WORKING  OFFSETS  AND  ENTRY  OFFSETS 
NOINT=  0 
SUNDEL=  0. 


I 


1=  INDEX  FOR  CORE 
GROUPS 


GET  INTERSECTIONS  SUBROUTINE, 
GETINT 

RETURNS  INTERSECTION  NUMBERS 
FOR  CORE  GROUP  (I)  IN  INTNO 
(1...24) 


ENTRY  CYCLE  SUBROUTINE 
ENTCYCL 

RETURNS  TIME  IN  THE  CYCLE 
AND  ENTRY  CYCLE  LENGTH  FOR 
CORE  GROUP  (1) 


ADD  10  SECOND  LEAD  TIME  TO 
TIME  IN  THE  CYCLE 
ENTT  =  JT(I)  +  0. 


MODULO  SUBROUTINE, 
PERI05 


MODULO  LEAD  TIME  IN  CYCLE 
0  >  ENTT  <  ECYCLE 


J  =  INDEX  FOR  INTERSECTIONS 
IN  CORE  GROUP 


s 


Routine  TRNPAR 
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9 


GET  WORKING  OFFSETS  SUBROUTINE 
GETWORK 


RETURNS  A-PHASE  WORKING  OFFSET 
FOR  INTERSECTION  INTNO(J) 


GET  CYCLE  SUBROUTINE,  GETCYCLE 

RETURNS  WORKING  CYCLE  LENGTH  (WCYCLE)  AND 
TIME  IN  CYCLE  (WRKT)  FOR  I  NTERSECTION  INTNO(J) 


I 


ADD  10  SECOND  LEAD  TIME 
TO  TIME  IN  THE  CYCLE 
WRKT  =  WRKT  +  10 


I 


MODULO  SUBROUTINE,  PERI05 


MODULO  LEAD  TIME  IN  CYCLE 
0  .  <  WRKT  <WCYCLE 


DELTA  BETWEEN  LEAD  TIME 

IN  THE  CYCLE  AND  WORKING  OFFSET 

WRKOFP  =  WRKT  -  WRKOF 


MODULE  SUBROUTINE,  PERI05 


MODULO  DELTA  WORKING  OFFSET 
0  s  WRKOFP  <WCYCLE 


GET  ENTRY  OFFSET  SUBROUTINE, 
GETEOFF 


RETURN  A-PHASE  ENTRY  OFFSET 
FOR  INTERSECTION  INTNO(J) 


DELTA  BETWEEN  ENTRY  LEAD  TIME 
IN  THE  CYCLE  AND  ENTRY  OFFSET 
ENTOFP=  ENTT-  ENTOF 


Routine  TRNPAR  (Continued) 
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£ 


DELOFF =  DIFFERENCE'BETWEEN  WORKING 
OFFSETS  AND  ENTRY  OFFSETS 
DELOFF=  WRKOFP  -  ENTOFP 


& 


SUM  DELOFF  AND  NUMBER  OF 

INTERSECTIONS 

SUMDEL=  SUMDEL  +  DELOFF 

NOINT=  NOINT  +  1 


© 


<D 


COMPUTE  TRANSITION  PARAMETER 
XMIN=  SUMDEL/NOINT 


J. 


1=  INDEX  FOR  CORE  GROUPS 


GET  INTERSECTIONS  SUBROUTINE, 
GETINT 

RETURNS  INTERSECTION  NUMBERS 
FOR  CORE  GROUPC )  IN  INTNO  (1...24) 


© 


Routine  TRNPAR  (Continued) 
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ENTRY  CYCLE  SUBROUTINE, 
ENTCYCL 


RETURNS  TIME  IN  THE  CYCLE 
AND  ENTRY  CYCLE  LENGTH 


ENTRY  CYCLE  LENGTH 
ECYCLE  =  JT(2) 


I 


J=  INDEX  FOR  INTERSECTIONS 
IN  CORE  GROUP  (I) 


I 


GET  ENTRY  OFFSET  SUBROUTINE, 
GETEOFF 


RETURNS  A-PHASE  ENTRY  OFFSET 
FOR  INTERSECTION  INTNO(J) 


ADD  TRANSITION  PARAMETER  TO  ENTRY  OFFSET 
EOFF=  IRTN 
rEOFF=  EOFF+XMIN 


0* 


I 


MODULO  SUBROUTINE, 
PERI05 


MODULO  NEW  ENTRY  OFFSET 
0  <  EOFF<ECYCLE 


I 


PUT  ENTRY  OFFSET  SUBROUTINE, 
PUTEOFF 


STORE  NEW  ENTRY  OFFSET  FOR 
INTERSECTION  INTNO(J)  IN  ENTRY 
OFFSET  ARRAY 


AVE  ALL  INTER- 
YES  y    SECTIONS  IN  THE  CORE  \NO_ 
^GROUP  BEEN  PROCE! 


Routine  TRNPAR  (Continued) 
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Routine  TRNPAR  (Concluded) 
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2.6.111    SUBROUTINE  GETBI 

2.6.111.1  Purpose 

Subroutine  GETBI  unpacks  a  bit  from  a  bit  string, 

2.6.111.2  Relationship 

This  subroutine  is  called  by: 
CALL  GETBI  (X,   I,   Y,   J) 

where: 

X,  I,   Y,  and  J  are  the  mnemonics: 

X  =  Bit  string  (Array  name) 

I    =  Bit  index  for  X 

Y  =  Output  word 

J   =  Interleaved/ non-interleaved  flag 

J  =  1  for  interleaved 

J  =  0  for  non-interleaved 

2.  6.  llli  3    Assumption 
None. 


239 


2. 6. 111. 4    Approach 

This  subroutine  is  called  at  Priority  Levels  1,    3,   and  4  of  the  FORTRAN 
software  and  is  therefore  required  to  be  triply  re-entrant.     The  working 
registers  are  stored  immediately  upon  entering  the  routine  and  are  restored 
immediately  prior  to  exiting.     The  storage  array  is  accessed  on  a  "last  in, 
first  out"  basis  and  is  large  enough  to  accommodate  three  sets  of  registers. 
The  output  is  right  justified  in  the  output  word,    Y. 

It  should  be  noted  that  the  index  I  is  a  FORTRAN  index  and  as  such,   must  be 
decremented  by  one  in  the  routine  to  obtain  the  proper  assembly  language 
index  to  access  the  bit  string. 

Provision  is  made  fo.   accessing  bit  arrays  which  are  word  interleaved. 
Some  arrays  are  interleaved  for  the  Sigma  5  system  because  when  inputting 
or  outputting  a  word  of  data  via  the  multiplexer,    two  consecutive  words  of 
buffer  storage  are  required.     The  first  word  is  the  command;  the  second 
word  contains  the  data. 


2.6.111.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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C    ENTER  J 


SAVE  REGISTERS 


I 


IWORD=-^  =  WORD 
INDEX  FOR  ARRAY  X 


IBIT=(I-1)-32*J^- 
=  BIT  INDEX  * 


Y=BIT  IBIT 
OF  WORD  IWORD 
OF  ARRAY  X 


I 


RESTORE  REGISTERS 


~T~ 

f  return) 


Subroutine  GETBI 


241 


2.6.  112    SUBROUTINE  GETBY  OR  UPBYT 

2.6.  112.  1    Purpose 

Subroutine  GETBY  unpacks  a  specified  byte  from  a  byte  array  and  places  it 
into  the  rightmost  byte  of  a  specified  word. 

2.6.112.2  Relationship 

This  subroutine  is  called  by: 

CALL  GETBY  (X,    I,    Y) 

where: 

X,    I  and  Y  are  the  mnemonics: 
X  =  Byte  array 
I  =  Byte  index  for  X 
Y  =  Output  word. 

2.6.112.3  Assumption 
None. 
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2.6.112.4    Approach 

This  routine  is  called  at  four  priority  levels  of  the  FORTRAN  software 
and  is  therefore  required 'to  be  quadruply  re-entrant.     The  working  registers 
are  stored  immediately  upon  entering  the  routine  and  are  restored  im- 
mediately prior  to  exiting.     The  storage  array  is  accessed  on  a  "last  in, 
first  out"  basis  and  is  large  enough  to  accommodate  four  sets  of  registers. 
The  output  is  right  justified  in  the  output  word,   Y. 

It  should  be  noted  that  the  byte  index  I  is  a  FORTRAN  index  and  as  such, 
must  be  decremented  by  one  in  the  routine  to  obtain  the  proper  assembly 
language  index  to  access  the  array. 

It  should  also  be  noted  that  the  Priority  Level  2  routines  call  the  routine 
by  the  mnemonic  UPBYT.     Therefore  UPBYT  is  defined  as  an  entry  point 
to  this  routine  by  giving  it  the  same  memory  address  as  GETBY.     This 
capability  exists  in  any  assembly  language. 


2.6.112.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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(enter) 
save  registers 


LOAD  l-TH 

BYTE  FROM  THE  BYTE 

ARRAY  X  INTO  Y 


i 


RESTORE 
REGISTERS 


,  1  % 

(return) 

Subroutine  GETBY 
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2.6.113    SUBROUTINE  HALF 

2.6. 113. 1  Purpose 

Subroutine  HALF  unpacks  a  specified  halfword  from  a  halfword  array  in  a 
32 -bit  computer  configuration. 

2.6.113.2  Relationship 

This  subroutine  is  called  by: 
CALL  HALF  (X,   I,   Y) 

where: 

X,   I,   and  Y  are  the  mnemonics: 

X  =  Halfword  array 

I    =  Halfword  index  for  X 

Y  =  Output  word 

2.6.113.3  Assumption 
None. 
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2.6.113.4    Approach 

This  subroutine  is  called  at  four  priority  levels  of  the  FORTRAN  software 
and  is  therefore  required  to  be  quadruply  re-entrant.     The  working  registers 
are  stored  immediately  upon  entering  the  subroutine  and  are  restored  im- 
mediately prior  to  exiting.     The  storage  array  is  accessed  on  a  "last  in, 
first  out"  basis  and  is  large  enough  to  accommodate  four  sets  of  registers. 
The  output  is  right  justified  in  the  output  word,    Y. 

It  should  be  noted  that  the  index  I  is  a  FORTRAN  index  and  as  such,   must  be 
decremented  by  one  in  the  routine  to  obtain  the  proper  assembly  language 
index  to  access  the  halfword  array. 

If  the  FORTRAN  software  is  implemented  on  a  16 -bit  computer,   this  routine 
should  simply  access  a  word  array,    each  word  having  16  bits. 


2.6.113.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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C  ENTER   J 


SAVE  REGISTERS 


I 


DECREMENT 


i 


Y=HALFWORD  I 
OF  ARRAY  X 


i 


RESTORE 
REGISTERS 

TUT 

(^RETURN  J 

Subroutine  HALF 
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2.6.114    SUBROUTINE  PKB1T 

2.6.114.1  Purpose 

Subroutine  PKBIT  packs  thirty-two  flag  words  of  data  into  32  bits  of  a 
specified  bit  array. 

2.6.114.2  Relationship 

This  subroutine  is  called  by: 
CALL  PKBIT  (X,   I,    Y) 

where: 

X,   I,   and  Y  are  the  mnemonics: 

X  =  Bit  array 

I    =  Word  index  for  X,   indicating  into  which  word  packing 
is  to  begin 

Y  =  Input  array  dimensioned  32  words. 

2.6.114.3  Assumption 
None. 
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2.6.114.4    Approach 

This  subroutine  is  called  only  at  Priority  Level  2   of  the    FORTRAN  soft- 
ware and  is  not  required  to  be  re-entrant. 

It  should  be  noted  that  the  index  I  is  a  FORTRAN  index  and  as  such,    must 
be  decremented  by  one  in  the  subroutine  to  obtain  the  proper  assembly  lan- 
guage index  to  access  the  bit  array. 

In  the  current  Sigma  5  configuration,    all  bit  arrays  accessed  by  the  Priority 
Level  2  routines  are  word  interleaved  and  thus  Subroutine  PKBIT  assumes 
word  interleaving.     It  should  be  noted  that  when  using  a  16 -bit  machine,    pro- 
vision must  be  made  in  the  routine  to  access  two  consecutive  words  of  the 
bit  array  X.     In  a  32-bit  machine,    only  one  word  must  be  accessed. 


2.6.114.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 


249 


C  ENTER   J 

IWORD  =  2*CI-1) 
=  WORD  INDEX 
FORX 


I 


X  IS  AN 

INTERLEAVED 

ARRAY 


IBIT=  0 


STORE  WORD  YCIBIT  +  1) 
IN  BIT  IBITOFWORD 
X  (IWORD) 


I 


IBIT=  IBIT+1 


Subroutine  PKBIT 
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2.6.115    SUBROUTINE  PKBYT 

2.6.115.1  Purpose 

Subroutine  PKBYT  extracts  the  rightmost  byte  of  a  word  and  stores  it  in  the 
Jth  byte  of  another  word. 

2.6.115.2  Relationship 

This  subroutine  is  called  by: 
CALL  PKBYT  (I,    J,    K) 
The  subroutine  is  called  by  all  Priority  Level  2  routines. 

2.6.115.3  Assumption 
None. 

2.6.115.4  Approach 

Indirect  addressing  is  used  to  access  word  K.     The  rightmost  byte  of  K  is 
then  stored  (indirect  addressing)  into  the  Jth  byte  of  word  I. 

2.6.115.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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(enter  ^ 


LOAD  R3  WITH  K 
LOAD  R4  WITH  J 
R4=  R4-1 


STORE  BYTE 
FROM  R3  INTO 


(return) 


Subroutine  PKBYT 
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2.6.116    SUBROUTINE  PKHW 

2.6.116.1    Purpose 

Subroutine  PKHW  extracts  the  rightmost  half  word  of  a  word  and  stores  it 
in  the  Jth  half  word  of  another  word. 


2.6.116.2    Relationship 

This  subroutine  is  called  by: 

CALL  PKHW  (I,    J,    K) 

The  subroutine  is  called  at  the  Priority  Level  2  routines  to  pack  data  which 
is  in  half  word  format. 


2.6.116.3  Assumption 
None. 

2.6.116.4  Approach 

Indirect  addressing  is  used  to  access  word  K.     The  rightmost  half  word  of 
K  is  then  stored  (indirect  addressing)  into  the  Jth  half  word  of  word  I. 

2.6.116.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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f  enter) 


LOAD  R3  WITH  K 
LOAD  R4WITH  J 
R4=  R4-1 


STORE  HALFWORD 
FROM  R3  INTO  I 

"XT 

f  RETURN  J 


Subroutine  PKHW 
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2.6.117    SUBROU  TINE  PTBI 

2.6. 117. 1  Purpose 

Subroutine  PTBI  packs  a  bit  into  a  bit  string. 

2.6.117.2  Relationship 

This  subroutine  is  called  by: 
CALL  PTBI  (X,   I,   Y,   J) 

where: 

X,   I,   Y,   and  J  are  the  mnemonics: 

X  =  Bit  string 

I    =  Bit  index  for  X 

Y  =  Input  word 

J   =  Interleaved/non-interleaved  flag 

J  =  1  for  interleaved 

J  =  0  for  non-interleaved 

2.6.117.3  Assumption 
None. 
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2.6.117.4    Approach 

This  subroutine  is  called  at   Priority  Levels  1,    3,    and  4  of  the  FORTRAN 
software  and  is  therefore  required  to  be  triply  re-entrant.     The  working 
registers  are  stored  immediately  upon  entering  the  subroutine  and  are  re- 
stored immediately  prior  to  exiting.     The  storage  array  is  accessed  on  a 
"last  in,    first  out"  basis  and  is  large  enough  to  accommodate  three  sets 
of  registers. 

It  should  be  noted  that  the  index  I  is  a  FORTRAN  index  and  as  such,   must 
be  decremented  by  one  in  the  routine  to  obtain  the  proper  assembly 
language  index  to  access  the  bit  string. 

Provision  is  made  for  accessing  bit  arrays  which  are  word  interleaved. 
Some  arrays  are  interleaved  for  the  Sigma  5  system  because  when  inputting 
or  outputting  a  word  of  data  via  the 'multiplexer,   two -consecutive  words  of 
buffer  storage  are  required.     The  first  word  is  the  command;  the  second 
word  contains  the  data. 


2.6.117.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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f  ENTER  ) 
SAVE  REGISTERS 


I 


IWORD=-^  =  WORD 
INDEX  FOR  ARRAY  X 


IBIT=(l-l)-32* 
=  BIT  INDEX 


m 


BIT  IBIT  OF  WORD 
IWORD  OF  ARRAY  X 
=  Y 


I 


RESTORE  REGISTERS 


(^RETURN  J 


Subroutine  PTBI 
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2.  6.  118    SUBROUTINE  PTBY  OR  PKBYT 

2.6. 118. 1  Purpose 

Subroutine  PTBY  (or  PKBYT)  packs  input  data  into  a  specified  byte  of  a  byte 
array. 

2.6.118.2  Relationship 

This  subroutine  is  called  by: 
CALL  PTBY  (X,   I,   Y) 

where: 

X,  I,   and  Y  are  the  mnemonics: 

X  =  Byte  array 

I    =  Byte  index  for  X 

Y  =  Input  word 

2.6.118.3  Assumption 
None. 

2.6.  118.4    Approach 

This  subroutine  is  accessed  at'  four  priority  levels  of  the  FORTRAN  software 
and  is  therefore  required  to  be  quadruply  re-entrant.     The  working  registers  are 
stored  immediately  upon  entering  the  subroutine  and  are  restored  immediately 
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prior  to  exiting.     The  storage  array  is  accessed  on  a  "last  in,   first  out" 
basis  and  is  large  enough  to  accommodate  four  sets  of  registers. 

It  should  be  noted  that  the  byte  index  I  is  a  FORTRAN  index  and  as  such, 
must  be  decremented  by  one  in  the  routine,  to  obtain  the  proper  assembly 
language  index  to  access  the  array. 

It  should  also  be  noted  that  the  Priority  Level  2  routines  call  the  routine 
by  the  mnemonic  PKBYT.     Therefore  PKBYT  is  defined  as  an  entry  point 
to  this  routine  by  giving  it  the  same  memory  address  as  PTBY.     This 
capability  exists  in  any  assembly  language. 


2.6.118.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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(enter) 
save  registers 

i 

1=1-1 

i 


STORE  Y  INTO  THE 
l-TH  BYTE  OF  THE 
BYTE  ARRAY  X 


I 


RESTORE  REGISTERS 


(returnj 


Subroutine  PTBY 
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2.6.119    SUBROUTINE  PTHLF 

2.6.119.1  Purpose 

Subroutine  PTHLF  packs  a  halfword  into  a  halfword  array  in  a  32-bit  com- 
puter configuration. 

2.6.119.2  Relationship 

This  subroutine  is  called  by: 
CALL  PTHLF  (X,  I,   Y) 

where: 

X,   I,   and  Y  are  the  mnemonics: 

X  =  Halfword  array 
I  =  Halfword  array 
Y  =  Input  word 

2.6.119.3  Assumption 
None. 

2.6.119.4  Approach 

This  subroutine  is  called  at  four  priority  levels  of  the  FORTRAN  software 
and  is  therefore  required  to  be  quadruply  re-entrant.     The  working  registers 
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are  stored  immediately  upon  entering  the  subroutine  and  are  restored  im- 
mediately prior  to  exiting.     The  storage  array  is  accessed  on  a  "last  in, 
first  out"  basis  and  is  large  enough  to  accommodate  four  sets  of  registers. 

It  should  be  noted  that  the  index  I  is  a  FORTRAN  index  and  as  such,   must 
be  decremented  by  one  in  the  routine  to  obtain  the  proper  assembly  language 
index  to  access  the  halfword  array. 

If  the  FORTRAN  software  is  implemented  on  a  16-bit  computer,  this  routine 
should  simply  access  a  word  array,   each  word  having  16  bits. 


2.6.119.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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C  ENTER   J 


SAVE  REGISTERS 


I 


HALFWORD  1-1 
OF  ARRAY  X 
=  Y 


1 


RESTORE  REGISTERS 


(return) 


Subroutine  PTHLF 
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2.6.120    SUBROUTINE  UNBIT 

2.6.120.1  Purpose 

Subroutine  UNBIT  unpacks  32  bits  of  data  from  a  bit  array  and  places  then 
into  32  words  of  working  storage. 

2.6.120.2  Relationship 

This  subroutine  is  called  by: 
CALL  UNBIT  (X,   I,    Y) 

where: 

X,   I,   and  Y  are  the  mnemonics: 

X  =  Bit  array 

I    =  Word  index  for  X,   indicating  from  which  word  unpacking 
is  to  begin 

Y  =  Output  array  dimensioned  32  words. 

2.6.120.3  Assumption 
None. 

2. 6. 120.4  Approach 

This  subroutine  is  called  only  at   Priority   Level  2   of  the   FORTRAN  soft- 
ware and  is  therefore  not  required  to  be  re-entrant. 
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It  should  be  noted  that  the  index  I  is  a  FORTRAN  index  and  as  such,   must 
be  decremented  by  one  in  the  routine  to  obtain  the  proper  assembly- 
language  index  to  access  the  bit  array. 

In  the  current  Sigma  5  configuration,  all'bit  arrays  accessed  by  the  Pri- 
ority Level  2  routines  are  word  interleaved  and  thus  Subroutine  UNBIT 
assumes  word  interleaving.     It  should  be  noted  that  when  using  a  16-bit 
machine,  provision  must  be  made  in  the  routine  to  access  two  consecutive 
words  of  the  bit  array  X.     In  a  32 -bit  machine,   only  one  word  must  be 
accessed. 


2.6.120.5    Listing 


See  Volume  III,    2GCF  Routine  Listings. 
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C  enter) 

IWORD  =  2*(l-l) 

=  WORD  INDEX 
FORX 


I 


WORD  =  X(IWORD) 
=  WORD  TO  BE 
UNPACKED 


IBIT=  0 


Y(IB1T  +  1)=  BIT 
IBITOFWORD 


IB1T  =  IBIT+1 


X  IS  AN 
INTERLEAVED  ARRAY 


Subroutine  UNBIT 
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2.6.121    SUBROUTINE  UPBYT 

2.6.121.1  Purpose 

Subroutine  UPBYT  extracts  the  Jth  byte  of  a  word  I  and   stores   it   in 
word  K. 

2.6.121.2  Re  lationship 

This  subroutine  is  called  by: 
CALL  UPBYT  (I,    J,    K) 
The  subroutine  is  called  at  all  Priority  Level  2  subroutines. 

2.6.121.3  Assumption 
None. 

2.6. 121.4  Approach 

Indirect  addressing  is  used  to  access  the  Jth  byte  of  word  I.     This  byte  is 
then  stored  (indirect  addressing)  in  word  K. 

2.6.121.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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f  ENTER   ") 


LOAD  R4WITH  J 

R4=  R4-1 

LOAD  R3  WITH  BYTE 


I 


STORE  WORD 
FROM  R3  INTO  K 


(^RETURN  J 


Subroutine  UPBYT 
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2.6.122    SUBROUTINE  UPHW 

2.6.122.1  Purpose 

Subroutine  UPHW  extracts  the  Jth  halfword  of  a  word  I  and   stores   it   in 
word  K. 

2.6.122.2  Relationship 

This  subroutine  is  called  by: 

CALL  UPHW  (I,    J,    K) 

UPHW  is  called  at  the  Priority  Level  2  subroutines  to  unpack  data  which 
is  in  halfword  format. 

2.6.122.3  Assumption 
None. 

2.6.122.4  Approach 

Indirect  addressing  is  used  to  access  the  Jth  halfword  of  word  I.     This 
halfword  is  then  stored  (indirect  addressing)  in  word  K. 

2.6.  122.  5    Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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f ENTER  J 


LOAD  R4  WITH  J 

R4=R4-1 

LOAD  R3  WITH  HALFWORD 


I 


STORE  WORD 
FROM  R3  INTOK 


(return) 


Subroutine  UPHW 
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2.6.123    SUBROUTINE  FMTNUM 

2.6.123.1  Purpose 

To  temporarily  hold  arguments.     FMTNUM  contains  no  executable  code. 

2.6.123.2  Relationship 

Subroutine  FMTNUM  is  not  called.     It  contains  no  executable  code. 

2.6.123.3  Assumption 
None. 

2.6.123.4  Approach 
N/A 

2.6.123.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 


271 


(  ENTER   J 

DATA  IN 
INTRNUM 

(ret 

urn} 

Subroutine  FMTNUM 
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2.6.124    SUBROUTINE  INTRNUM 

2.6.124.1  Purpose 

Temporarily  holds  arguments.     Contains  no  executable  code, 

2.6.124.2  Relationship 

Data  storage  area.     Contains  no  executable  code. 

2.6.124.3  Assumption 
None.  ' 

2.6.124.4  Approach 
N/A 

2.6.124.5  Listing 

See  Volume  III,    2GCF  Routine  Listings. 
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r  ENTER    J 


ACCESS  DATA 


(return) 


Subroutine  INTRNUM 
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SECTION  3 
GLOSSARY 
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Definition  of  Variable  Names  Used  in  Subroutine  LINKD 


Variable 

Name 

Format 

GCNTSI 

byte 

GCONT 

byte 

G1LNK 

byte 

INT 

word 

INDEX 

word 

ISB 

word 

JNDEX 

word 

LINK 


LINK2 


LNDEX 


word 


word 


word 


NCIGRP 

byte 

N2GCL 

byte 

ONDEX 

word 

Definition 

Node  starting  index  for  subnetwork. 

Array  of  node  numbers  for  all  groups. 

List  of  first  generation  link  numbers 
associated  with  each  core  group. 

Output  argument  containing  list  of  group 
numbers  in  subnetwork. 

Output  argument  giving  list  of  first  upstream 
node  numbers  associated  with  each  link  in 
the  subnetwork. 

List  of  core  groups  within  a  subnetwork. 

Output  argument  giving  list  of  first  down- 
stream subnetwork  node  numbers  asso- 
ciated with  each  link  in  the  subnetwork. 

Output  argument  giving  first  generation 
data  link  number  associated  with  each 
subnetwork  link. 

Output  argument  giving  list  of  second 
generation  link  numbers  associated  with 
each  link  in  the  subnetwork. 

Output  argument  giving  list  of  second  up- 
stream node  numbers  associated  with  each 
link  in  the  subnetwork. 

Number  of  controllers  within  each  core 
group. 

Number  of  links  associated  with  each  core 
group. 

Output  argument  giving  list  of  second  down- 
stream subnetwork  node  numbers  asso- 
ciated with  each  link  in  the  subnetwork. 
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Definition  of  Variable  Names  Used  in  Subroutine  LPHASV 


Variable 

Name 

Format 

GAP 

byte 

GAPC 

byte 

GAM 

byte 

GBC 

byte 

GBM 

byte 

GCM 

byte 

GGGG 

4  bytes/ 
controllers 

GGGT 


4  bytes/ 
controllers 


INT 

word 

K3P 

word 

K3PHI 

byte 

LKPTR 

byte 

NQL 

byte 

ssss 

4  bytes/ 
controller 

Definition 

Time  available  for  A  and  B  phases. 

Time  available  for  A,   B,   and  C  phases. 

Minimum  A  phase  time. 

Time  available  for  B  and  C  phases. 

Minimum  B  phase  time. 

Minimum  C  phase  time. 

Byte  1:    Minimum  allowable  durations  of  A 
phase  green. 

Byte  2:    Minimum  allowable  duration  of  B 
phase  green. 

Byte  3:    Minimum  allowable  durations  of  C 
phase  green. 

Byte  4:    Minimum  allowable  minor  adjustable 
interval  duration. 

Byte  1:     Time  available  for  A  phase. 

Byte  2:    Time  available  for  B  phase. 

Byte  3:    Time  available  for  C  phase. 

Byte  4:    Minimum  transitional  cycle  length 
duration. 

Intersection  number. 

Type  of  CIC  control. 

Type  of  CIC  control. 

Link  phase  table. 

Number  of  multi-detector  links. 

Byte  1:    Summation  of  all  fixed  interval  in 
A  phase. 

Byte  2:    Summation  of  all  fixed  interval  in 
B  phase. 

Byte  3:    Summation  of  all  fixed  interval  in 
C  phase. 

Byte  4:    Summation  of  all  fixed  interval 
and  the  minimum  duration  of  all  adjustable 
intervals. 
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Definition  of  Variable  Names  Used  in  Subroutine  SMSCYL 


Variable 

Name 

Format 

CYCLE 

word 

ENTCYC 

byte 

ICYC 

word 

INT 

word 

IPOSN 

word 

KCL 

word 

SMTABC 

4  bytes/ 
controller 

SSSS 

4  bytes/ 
controller 

SUMA 

word 

SUMB 

word 

SUMC 

word 

SUMT 

word 

Definition 

Cycle  length  (in  sec^nds  )  for  controller  INT. 
Cycle  length  code  for  each  controller. 
Cycle  length  code  for  controllers  INT. 
Controller  number. 

Index  to  SSSS  array  for  controller  INT. 
Cycle  length  corresponding  to  ENTCYC  code. 
Same  as  SSSS. 

Byte  1:    Summation  of  all  fixed  intervals  in 
A  Phase. 

Byte  2:    Summation  of  all  fixed  intervals  in 
B  Phase. 

Byte  3:    Summation  of  all  fixed  intervals  in 
C  Phase. 

Byte  4:    Summation  of  all  fixed  intervals  and 
the  minimum  duration  of  all  ad- 
justable intervals  which  are  not 
part  of  A,   B,  or  C  phases. 

Summation  of  all  fixed  intervals  in  A-phase. 

Summation  of  all  fixed  intervals  in  B-Phase. 

Summation  of  all  fixed  intervals  in  C -Phase. 

Summation  of  all  fixed  intervals  and  the 
minimum  duration  of  all  adjustable  intervals 
in  a  controller  cycle  which  are  not  part  of 
A,  B,  or  C  Phase  Green. 
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APPENDIX  A 
SOFTWARE  MAINTENANCE  PROCEDURE 


Modifications  to  the  2GCF  software  is  accomplished  in  three  steps: 

•  Update  the  source  tape(s) 

•  Compile  or  assemble  the  new  source 

•  Load  and  link  the  new  object  files. 

These  procedures  are  rather  straightforward,  and  should  be  easily  performed 
by  personnel  familiar  with  SIGMA  5  operator  procedures. 

1.     Updating  Source  Tapes 

The  2GCF  source  language  resides  on  three  magnetic  tapes,  labeled  "MAIN 
FORTRAN  ROUTINES*',    "OVERLAY  FORTRAN  ROUTINES",  and  "ASSEMBLY 
ROUTINES".     The  order  and  contents  of  each  tape  are  shown  in  Tables  2, 
3  and  4,   respectively. 

To  change  a  source  program,  the  system  routine  UPDATE  has  been  included, 
in  the  2GCF  system  library.     A  user's  guide  to  the  procedures  for  using 
UPDATE  appears  in  Appendix  B. 

The  output  from  the  UPDATE  procedure  will  be  the  modified  compiler/ 
assembler-ready  code,   which  is  processed  by  the  next  step  in  the  Software 
Maintenance  Procedure.     UPDATE  provides  listings  of  the  modified  source 
programs,  which  eliminates  the  need  for  using  the  lengthy  and  time  con- 
suming compiler/assembler  listing  options. 
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Table  2.     Main  FORTRAN  Tape  Contents 

Deck  No.  Subroutine  Length  (Lines) 

1  RTNC  71 

2  RTNE  240 

3  RTNF  83 

4  RTNA1  172 

5  UPKMC  23 

6  UPBL2  75 

7  PKBL2  73 

8  BLOCK  DATA  29 

9  RTNO  46 

10  UNPO  17 

11  PCKO  17 

12  RTNP  111 

13  UNPP1  18 

14  UNPP2  20 

15  PCKP1  18 

16  RTNR  58 

17  RTNS  180 

18  UNPS  39 

19  RTNT  222 

20  PCKT  17 

21  UNPT  26 

22  MBYTS  18 

23  RTNU  97 

24  RTNV  267 

25  PSUB2  21 

26  UPKV  57 

27  GREEN  45 

28  GDSUB  45 

29  RTNX  189 

30  DECON  44 

31  HFMT  78 

32  MHEAD  33 

33  STMSG  71 

34  RTNJ  163 

35  BCDTH  24 

36  SDL  26 

37  CLS  29 

38  TCSLP  28 

39  CDL  23 

40  RTNK  823 

41  KFMT  139 

42  SRCFT  40 

43  CRTOUT  75 
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Table  2.     Main  FORTRAN  Tape  Contents  (Concluded) 

Deck  No.  Subroutine  Length  (Lines) 

44  KFTU  201 

45  CIC  42 

46  E1SRT  57 

47  SBYTC  22 

48  CNVT  117 

49  SCODE  56 

50  DCMAL  26 

51  CRLF  22 

52  DISNM  36 

53  DISIN  52 

54  SETCR  38 

55  RTNM  711 

56  MACRO  86 

57  PERI05  15 

58  FMODE  23 

59  RTN5  128 

60  LDMOE5  34 

61  TRFLG  20 

62  LLGTHF  13 

63  RTNA5  71 
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Table  3.     Overlay  FORTRAN  Routines  Tape  Contents 


Deck  No.  Subroutine  Length  (Lines) 

1  RTSND  38  6 

2  IRTN  33 

3  LPHASV  44 

4  SMSCYL  30 

5  LSTSQS  427 

6  LINKD  75 

7  SRCH  12 

8  CORRECT  49 

9  LOAD  AT  176 

10  INTFC  168 

11  TRNPAR  103 

12  MACA  29 

13  PUTEOFF  46 

14  GETINT  20 

15  PREDICT  186 

16  RTNH  222 

17  STAT  434 

18  LKDAT  383 

19  OUT2G  99 

20  PATGENT  354 

21  SMINT  27 

22  RTNW  248 

23  RTNQ  40 

24  RTNQ1  47 

25  MOVER  91 

26  RTNI  441 

27  PRPRED  177 
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Table  4.     Assembly  Language  Tape  Contents 


Deck  No, 

1 
2 
3 

4 

5 


6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 


Subroutine                 Length  (L 

RTNA 

797 

RTNB 

135 

CALOP,  CLOSE 

91 

ERROR,  FRAD,   OPEN 

UPBYT,   PKBYT,   UPHW, 

8 

PKHW 

GET  BY,   PTBY,   GETBI, 

320 

PTBI,   HALF,   PTHLF, 

UNBIT,   PKBIT 

MSUB 

839 

FMTNUM,  INTRNUM 

9 

MSUB2 

5 

WSUB 

92 

WCODE,  WTEST,   IR3 

10 

CHECK 

18 

OUTAP 

62 

PRNTF 

85 

CAL15,  STLP 

10 

CALRL,  STRLP 

10 

TSTAP,  TYPE  OUT 

40 

TSTLP 

61 

TIMEH 

58 

INTRP 

26 

TODAY 

27 

TRIG 

13 

RTNY 

2258 

BUFFERIN/OUT 

40 

OUTTBF 

75 

DEBUGLV5 

28 
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2.     Compiling /Assembly  UPDATE  Output 

The  output  from  the  UPDATE  procedure  is  ready  for  input  to  the  compiler 
or  assembler.     The  compile /assemble  control  card  sequence  is: 

!  ASSIGN  (M:SI,  7TAE1 ),  BIN,  PACK,  (RECL,  80) 
!  ASSIGN  (M:BQ,Dn,  filename) 
!  FORTRAN  SI,  BO,  NS,  S 

where  Dn  is  the  RAD  area  which  contains  the  binary  output  file,   and 
filename  is  the  binary  output  file  name. 

The  appropriate  control  cards  for  a  complete  recompilation/ reassembly, 
are  shown  in  Figure  2. 

The  binary  object  files  produced  by  the  compile/assemble  procedure  con- 
tain the  object  decks  for  groups  of  source  programs  as  described  in  Table 
5. 


3.     Load  and  Link 

After  performing  all  UPDATE  and  compile/assemble  operations,  the  new 
binary  files  must  be  loaded  and  linked  to  form  the  modified  executable  file 
KING.     The  control  card  sequence  to  perform  this  task  is: 

1JOB 

!  ATT 

! PAUSE  SYC 

IOLOAD  (FORE,  2A00),  (FILE,  FP,  KING),  (TEMP,  1000,  1000), 

(UDCB,  5),  (TASKS,  91); 

:(MAP,  ALL,  ALPHA) 

:ROOT  (FILE,  Dl,  SUB  1,  EOD),  FILE,  Dl,  SUB2,  EOD),  (FILE,  D4,  COC, 

EOD),  (ENTRY,  RTNY) 

:SEG  (LINK,  1,  ONTO,  0),  (FILE,  D4,  LST,  EOD) 

:SET  (LINK,  2,  ONTO,  0),  (FILE,  D4,  RTN,  EOD) 
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SEG(LINK,  3,  ONTO,  0),  (FILE,  D4,  PRE,  EOD) 
SEG  (LINK,  4,  ONTO,  0),  (FILE,  D4,  INT,  EOD) 
SEG  (LINK,  5,  ONTO,  0),  (FILE,  D4,  RTH,  EOD) 
SEG  (LINK,  6,  ONTO,  0),  (FILE,  D4,  RTW,  EOD) 
SEG  (LINK,  7,  ONTO,  0),  (FILE,  D4,  RTQ,  EOD) 
SEG  (LINK,  8,  ONTO,  0),  (FILE,  D4,  RTI,  EOD) 
ASSIGN(M:DO,  LPA02),  (VFC) 
ASSIGN  F:ll,  D4,COEFF) 
ASSIGN(F:12,  D4,  STAT  1) 
ASSIGN  (F:13,D4,STAT2) 
ASSIGN  (F;14,D4,  STAT3) 
ASSIGN  (F:15,  D4,STAT4) 
ASSIGN  (F:16,  D4,STAT5) 
ASSIGN  (F:17,  D4,STAT6) 
ASSIGN  (F:18,D4,COR) 
I  FIN 
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JOB 

ATTEND 

PAUSE  SYC 

PAUSE  MOUNT  ASSY  LANG 

ASSIGN  (F;ll,  7TAE0),  BIN,  PACK,  (RECL,80) 

ASSIGN  (F:12,  7TAE1),  BIN,  PACK,  (RECL,  8  0) 

UPDATE 

ASSEMBLY  LANGUAGE  (followed  by  appropriate  UPDATE  commands) 
EOD 

ASSIGN  (M:SI,  7TAE1),  BIN,  PACK,  (RECL,  80) 
ASSIGN  (M:B0,  D1,SUB1) 
AP  SI,B0,AC(D1),  BA 
REWIND  7TAE1 

PAUSE  MOUNT  OVERLAY  TAPE 
UPDATE 

OVERLAYS  (followed  bv  appropriate  UPDATE  commands) 
EOD 

ASSIGN  (M:SI,  7TAE1),  BIN,  PACK,  (RECL,  80) 
ASSIGN  (M:B0,  D4,RTN) 
FORTRAN  SI,  B9,  NS,  S,  BC(4) 
ASSIGN  (M:SI,  7TAE 1),  BIN,  PACK,  (RECL,  80) 
ASSIGN  (M:B0,  D4,  LST) 
FORTRAN  SI,  BO,  NS,  S,  BC(5) 
ASSIGN  (M:SI,  7TAE1),  BIN,  PACK,  (RECL,  80)- 
ASSIGN  (M:B0,  D4,INT) 
FORTRAN  SI,  BO,  NS,S,  BC(5) 
ASSIGN  (M:SI,  7TAE 1),  BIN,  PACK,  (RECL,  80) 
ASSIGN  (M:B0,  D4,PRE) 
FORTRAN  SI,  B0,  NS,  S,  BC(  1) 
ASSIGN  (M:SI,  7TAE1),  BIN,  PACK,  (RECL,  80) 
ASSIGN  (M:B0,  D4,  RTH) 
FORTRAN  SI,  B0,  NS,  S,  BC(4) 
ASSIGN  M:SI,  7TAE1),  BIN,  PACK,  (RECL,  80) 
ASSIGN  (M:B0,  D4,RTW) 
FORTRAN  SI,  B0,  NS,  S,  BC(3) 
ASSIGN  (M:SI,  7TAE1),  BIN,  PACK,  (RECL,  8  0) 
ASSIGN  (M:B0,  D4,RTQ) 
FORTRAN  SI,  B0,  NS,  S,  BC(3) 
ASSIGN  (M:SI,  7TAE1),  BIN,  PACK,  (RECL,  80) 
ASSIGN  (M:B0,D4,RTI) 
FORTRAN  SI,  B0,  NS,  S,  BC(2) 
REWIND,  7T AE 1 

PAUSE  MOUNT  MAIN  FORTRAN  TAPE 
UPDATE 

MAIN  FORTRAN  (followed  by  appropriate  UPDATE  commands) 
EOD 

ASSIGN  (M:SI,  7TAE1),  BIN,  PACK,  (RECL,  8  0) 
ASSIGN  (M:B0,  Dl,  SUB2) 
FORTRAN  SI,  BO,  NS,  S,  BC 

Figure  2.     Control  Card  Format  for  Complete  Recompilation/Reassembly 
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!  REWIND  7TAE1 

I  OLOAD  (FORE,  2A00),  (FILE,  FP,  KING),  (TEMP,  1000,  1000),  (UDCB,  5) 

TASKS,  9); 

:(MAP,  ALL  ALPHA) 

:ROOT  (FILE,  Dl,  SUB1,  EOD),  (FILE,  Dl,  SUB2,  EOD),  (FILE,  D4,  COC, 

EOD),  (ENTRY,  RTNY) 

SEG  (LINK,  l,ONTO,  0),  (FILE,  D4,  LST, EOD) 

SEG  (LINK,  2,  ONTO,  0),  (FILE,  D4,  RTN,  EOD) 

SEG  (LINK,  3,  ONTO,  0),  (FILE,  D4,  PRE,  EOD) 

SEG  (LINK,  4,  ONTO,  0),  (FILE,  D4,  INT,  EOD) 

SEG  (LINK,  5,  ONTO,  0),  (FILE,  D4,  RTH,  EOD) 

SEG  (LINK,  6,  ONTO,  0),  (FILE,  D4,  RTW,  EOD) 

SEG  (LINK,  7,  ONTO,  0),  (FILE,  D4,  RTQ,  EOD) 

SEG  (LINK,  8,  ONTO,  0),  (FILE,  D4,  RTI,  EOD) 

ASSIGN  (M:D0,  LPA02),  (VFC) 

ASSIGN  (F:ll,D4,COEFF) 

ASSIGN  (F:12,  D4,STAT1) 

ASSIGN  (F:13,D4,STAT2) 

ASSIGN  (F:14,D4,STAT3) 

ASSIGN  (F:15,  D4,  STAT4) 

ASSIGN  (F:16,D4,STAT5) 

ASSIGN  (F:17,  D4,  STAT 6) 

ASSIGN  (F:18,D4,COR) 
1FIN 


Figure  2.     Control  Card  Format  for  Complete 
Re  compilation/  Reassembly 
(Concluded) 
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Table  5.     Binary  Object  File  Definitions 


File  /RAD  Area 

RTN/D4 

LST/D4 

INT/D4 

PRE/D4 

RTH/D4 

RTW/D4 

RTQ/D4 

RTI/D4 

SUB1/D4 

SUB2/D1 


Subroutines  Contained  in  File 

RTSND,   IRTN,    LPHASV,  SMSCYL 

LSTSQS,   LINKD,  SRCH,  CORRECT,   LOADAT 

INTFC,  TRNPAR,   MAC  A,   PUTEOFF,   GETINT 

PREDICT 

RTNH,  STAT,  LKDAT,  OUT2G 

PAT  GENT,  SMINT,  RTNW 

RTNQ,  RTNQ1,  MOVER 

RTNI,  PRPRED 

All  routines  in  Table  4 

All  routines  in  Table  2 
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APPENDIX  B 
SOURCE  TAPE  UPDATING 


1.     General  Specifications 

•  First  data  card  is  a  title  card  (printer  at  top  of  index) 

•  Command  Card  Format: 

I,   M,   J,  K,   L, 

I    =  Command  Code 

0  =  create 

1  =  delete/ replace 

2  =  insert 

3  =  add  to  end  of  file 

4  =  compile  deck  list 

5  =  list  deck  list 

6-9  =  same  as  0-3  only  code  is  converted  from  026  to  029 

M  =  Deck  Number 

J    =  Card  Number 
K  =  Delete  Count 
L   =  Insert  Count 

For  type  4  and  5  cards,  the  decks  to  store  or  list  are  specified  on  the 
following  card  (up  to  16).    A  special  option  is  available  for  consecutive 
decks,  as  a  pair  of  the  form  (M,   -N)  on  the  deck  list  card  causes  output 
of  decks  M  through  N. 
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2.     Special  Options 

Type  2  Commands 

•  If  J  is  negative,  the  insert  is  before  the  card  specified. 

•  If  K  is  negative,  the  insert  count  (L)  specifies  the  number  of 
decks  to  insert. 

Type  0  and  3  Commands 

•  M  not  zero  causes  input  to  be  read  from  channel  specified  by 
M.  (This  file  or  tape  must  be  specified  by  an  J  ASSIGN  card; 
e.g.,    !  ASSIGN  (F  =  1,   Dl,   DATA) 

Type  4  Commands 

•  M  not  zero  identifies  file  to  which  decks  are  written  (Default  =  1) 
(do  not  use  11,    12)  these  must  be  specified  by  I  ASSIGN  cards. 

Type  4  and  5  Commands 

•  If  J  is  not  zero,   output  to  new  tape  is  suppressed.     (May  need 
1  ASSIGN  card  (F  =  12,  CO) 

Type  5  Commands 

•  If  type  5  card  is  absent,   all  changed  or  added  decks  will  be 
listed. 

•  If  M  is  not  zero,  tapes  will  not  be  rewound  following  an  update 
run.  (This  option  may  be  used  to  store  several  update  files  on 
the  same  tape. ) 

•  If  the  first  entry  on  the  deck  list  card  is  zero,  all  changed  decks 
will  be  listed  in  addition  to  those  specified. 
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3.  Constraints 

•  Type  4  and  5  cards  when  used  should  be  the  first  command  cards. 

•  A  maximum  of  200  decks  may  be  stored  on  a  single  file. 

•  A  maximum  of  10  input/ compile  files  may  be  specified. 

•  If  an  entire  deck  is  deleted  so  that  all  following  decks  are 
renumbered,  the  change  deck  listing  option  will  not  list 
these  decks. 

4.  Procedures 

Create  a  New  Tape  --  Command  Card:    0,   M, 

M  =  0  -  Input  cards  from  card  reader. 
M  /  0  -  Input  from  channel  M. 

Replace  /Delete  Specified  Cards  --  Command  Card:     1,   M,  J,  K,   L, 

M  =  Deck  number 
J  =  Card  number 

K  =  Number  of  cards  to  delete  starting  at  M,   J  (K      0) 
L  =  Number  of  new  cards  to  be  inserted. 

Insert  New  Cards /Decks  --  Command  Card:    2,   M,   J,   K,   L, 

M  =  Deck  number 
J  =  Card  number 
K  =  0  -  Insert  cards 
K  \0  -  Insert  decks 
L  *  Card/deck  count  to  insert 
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NOTE:    J  <  0  causes  the  insert  to  be  before  card  M,    -J, 

Add  Decks  to  the  End  of  the  Old  File  —  Command  Card:    3,   M, 

M  =  0  -  Input  decks  via  card  reader 
M  ^  0  -  Input  from  Channel  M 

Write  Specified  Decks  to  Specified  Files  --  Command  Card:    4,   M,   J,  K, 
Following  Card:    D  ,   D  ,    ---,   D  ,     n  <:  16 

M  =  Output  file  specification 
(default  =  1) 

J    f  0  -  No  new  tape  will  be  produced 

K  /  0  -  File  output  will  be  binary  if  M  <  10,   otherwise  BCD. 

D.  =  Deck  numbers  to  be  written  to  file  M.     If  D.  <  0,    decks  D.    , 
i  l  l-l 

through  -D.  will  be  written. 
NOTE:    More  than  one  type  4  card  may  be  used  if  needed. 
Print  Specified  Decks  on  Line  Printer  --  Command  Card:    5,  M,  J, 
Following  Card:    D       D  ,    ---,D  ,    n  <  16 

J.  d*  II 

M  ^  0  Inhibit  rewind  of  old/new  tapes  permits  multiple  files  on 
the  same  tape. 

J  ^  0  Inhibit  production  of  a  new  tape.     (List  only  option) 

D.  Decks  to  be  printed.     If  D.  <  0,    decks  D-_-.  through  -D. 

will  be  printed. 

NOTE:    If  no  type  5  card  is  used,   only  decks  changed  will  be  listed.     If  a 

type  5  card  is  included,   D    must  be  zero  to  cause  changed  decks  to 
be  printed. 
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026  Card  Input  --  To  use  026  code  inputs,   increment  command  codes  by  6 

(except  codes  4  and  5);  i.e.,   Code  6  =  Code  0  (Create) 

Code  7  =  Code  1  (Replace /delete) 
Code  8  =Code  2  (Insert) 
Code  9  =  Code  3  (Add-File) 


5.     Constraints /Special  Features 

A  maximum  of  200  decks/file  on  the  magnetic  tape. 

A  maximum  of  10  input /output  files  per  run. 

The  first  card  following  the  update  run  card  (I  UPDATE)  must  be  a  title  card. 

If  the  first  card  in  a  program  has  one  of  the  following  forms: 

C###*  ###    DECK  NAME 
*  ###*  ###    DECK  NAME 
(#  =  space) 

Then  the  deck  name  will  be  printed  on  the  output  summary  list.     'DECK  NAME' 
may  be  up  to  64  characters.     The  'END'  card  of  a  program  deck  should  have 
at  least  7  spaces  following  the  word  'END'. 


*U.S.  GOVERNMENT  PRINTING  OFFICE  :  1977  0-731-586/106 
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